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Nuclear Magnetic Resonance in Indium Antimonide 
I. The Effect of Impurities; 
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Services Electronics Research Laboratory, Baldock, Herts. 
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ABSTRACT 


A study has been made of the effect of impurities on the !5In nuclear 
magnetic resonance in indium antimonide. In n-type (Te-doped) and 
p-type (Zn-doped) samples, the amplitude and width of the resonance are 
very sensitive to the presence of donors or acceptors, a concentration of 
101° em~$ reducing the maximum of the derivative of the absorption curve 
by a factor of 10. This is ascribed to broadening arising from the interaction 
of the nuclear quadrupole moment with the Coulomb field of the ionized 
impurity atoms. The effect is slightly less pronounced in p-type than in 
n-type samples, which can be explained by the fact that holes are more 
effective than electrons in screening the ionic charges, because of their greater 
effective mass. Assuming the field gradient at a distance r from an ion to 
be given by 0?V/0r?=2fe/r>, B must be ~350 to fit the experiments, which 
implies an ‘ anti-shielding factor’ of ~1000. It is suggested that this large 
value is associated with the smallness of the forbidden gap in InSb. The 
replacement of a small fraction of In atoms by Ga in an InSb—GaSb mixed 
crystal has a qualitatively similar effect on the 1°In resonance, but the 
concentration of Ga atoms must be about 30 times that of Te or Zn for the 
same reduction in resonance amplitude. Since the Ga atoms are electrically 
neutral, their effect is attributed to strain set up in the InSb lattice by the 
disparity between the sizes of the Ga and In atoms. 

A search for a possible Knight shift in n- and p-type InSb containing 
101° carriers cm~* revealed no deviation in either case from the 11°In resonance 
frequency in intrinsic material, within an experimental accuracy of | in 104. 
The negative result in n-type material is attributed to the small effective mass 
of the electrons, which reduces the Pauli susceptibility, while the p-type 
result provides possible evidence for the p-character of hole wave functions. 


§ 1. INTRODUCTION 


In recent years nuclear magnetic resonance has become widely used as a 
tool for solid state investigations. In particular, an extensive review 
article by Bloembergen (1955) has drawn attention to its potentialities for 
studying defects in crystals. The possibility of using nuclear resonance for 
this purpose arises because the precise frequency at which a particular 
nucleus resonates depends upon the magnetic (and in some cases electric) 
field which exists at the site of that nucleus; although mainly determined 
by the external magnetic field, this local field shows small deviations 
which depend on the nuclear environment. 
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Comparatively little work has been done on nuclear resonance in 
semiconductors. This is mainly because the most widely studied 
semiconductors, silicon and germanium, do not readily lend themselves 
to nuclear resonance investigations since the isotopes of these elements 
which possess nuclear magnetic moments occur in low abundance. On 
the other hand, the semiconducting intermetallic compounds formed 
from Group II and Group V elements are particularly suitable for nuclear 
resonance investigations since the nuclei of all these elements are mag- 
netically active. The present paper describes a study of indium antimo- 
nide using nuclear resonance techniques. 

After the work was begun, a paper was published by Shulman et al. (1955) 
on nuclear magnetic resonance in pure indium antimonide and gallium 
antimonide. The research described here does not overlap their work 
since they confined their studies to crystals as near to perfection as they 
could obtain, whereas the present paper is mainly concerned with the 
effect of deliberately introduced imperfections. Preliminary accounts 
of these experiments have already appeared (Rhoderick 1956, 1957). 


§ 2. EXPERIMENTAL METHOD 


The nuclear magnetic resonance was observed with a radio-frequency 
spectrometer of the type developed by Pound, Knight and Watkins (see, 
for example, Pound 1952), used in conjunction with a permanent magnet 
having a field of 3850 gauss. The normal technique of using a lock-in 
amplifier was adopted, so that the derivative of the absorption (dy”"/dv) was 
recorded, the line shape being plotted by varying the frequency at constant 
magnetic field. Measurements at liquid oxygen temperature were made 
using a cryostat of the heat-leak variety similar to that described by 
Gutowsky et al. (1953). 

One experimental difficulty may be mentioned here. The resistivity of 
highly extrinsic indium antimonide at room temperature may be as low as 
10-4 ohm cm, for which the skin depth at a frequency of 5 Mc/s is about 
0-2mm. Not only does the intensity of the absorption fall if the radio- 
frequency field does not completely penetrate the specimen, but distortions 
of the line shape may occur (Bloembergen 1952, Chapman et al. 1957). 
The difficulty may be overcome by using powdered specimens such that 
the particle size is much smaller than the skin depth. Initially these 
were embedded in an insulating matrix of paraffin wax, but this was 
subsequently found to be unnecessary, probably because of the high contact 
resistance. The resistivity of intrinsic indium antimonide at room tem- 
perature is sufficiently high for single crystals in the form of plates about 
1 mm thick to be used. 

The indium antimonide used in these experiments was zone-refined 
material which was intrinsic at room temperature. Hall effect measure- 
ments at liquid air temperature gave the excess acceptor concentration 
as 10 cm-%, the estimated total impurity concentration being less than 
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10%°cm~*. ‘ Doped’ specimens were prepared by. using tellurium and 
zine as donor and acceptor elements respectively. In all cases the carrier 
concentration obtained from Hall effect measurements agreed within 10° 
with the donor or acceptor concentration estimated from the weight of 
tellurium or zine added. 


§ 3. RESULTS FoR PurE INprtum ANTIMONIDE 


A brief investigation of nuclear magnetic resonance in undoped indium 
antimonide was made, the results being in substantial agreement with those 
reported by Shulman ef al. (1955). The shape of the "In absorption line 
in a single crystal is independent of the orientation of the crystal axes 
relative to the magnetic field, and its width is much greater than that to be 
expected from dipolar broadening. The width between points of maximum 
slope was found to be 9-9 + 0-4 ke/s, which corresponds to 10-6 + 0-4 gauss 
for an experiment at fixed frequency, compared with the value 9-0 + 0-2 
gauss given by Shulman etal. (1955). The line shape is almost symmetrical, 
the slight departure from symmetry being attributable to the presence 
of the 1°JIn isotope, which is 4° abundant, and which has a resonant 
frequency 7-5 ke/s lower than that of "*In in a field of 3850 gauss. The 
high frequency half of the line, which is practically free from the distorting 
effect of the %In resonance, is very closely Gaussian in shape. Shulman 
et al. (1955) have shown that this can be explained in terms of the indirect 
exchange coupling between nuclei first postulated by Ruderman and 
Kittel (1954) to explain anomalously broad lines in metals. This coupling 
mechanism, which was subsequently extended to the case of insulators by 
Bloembergen and Rowland (1955), leads to an interaction of the form 
A,,|,.1, between pairs of nuclei, | being the nuclear spin vector. The 
121Sb resonance is much broader (about 17 ke/s between points of maximum 
slope) and weaker, and no precise measurements were made on this line. 

Since the !In nucleus possesses a large quadrupole moment, it is 
possible that the shape of the resonance may be affected by quadrupolar 
interactions. If a nucleus having a quadrupole moment is placed in a 
non-uniform electric field, the effect of the field gradient is to split the 
nuclear resonance line into 27 equally spaced components. If 2/ is an odd 
integer, as in the case of !*In (J = 9/2), the central component is unchanged 
in frequency, the remaining 27 — 1 components, known as satellites, being 
grouped symmetrically about this central component. In crystals having 
the zinc-blende structure, such as indium antimonide, the lattice symmetry 
is such that the gradient of the crystalline field, grad-E, vanishes at the 
position of any nucleus, so that there should be no nuclear quadrupole 
interaction in a perfect crystal. However, as Watkins and Pound (1953) 
have pointed out, the symmetry of the field may be altered by the in- 
homogeneous strains which accompany dislocations. The result of these 
strains is to produce field gradients which vary in magnitude and orientation 
from one nucleus to another, so that the satellites become smeared into a 
continuous band, leaving the central component unaffected. According 
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to Watkins and Pound, the dislocation density in the best crystals of 
potassium bromide and iodide obtainable is such that the satellites of the 
bromine and iodine nuclear resonances are entirely obliterated, the observed 
intensity corresponding to that of the central component only. Since 
it is of importance to know whether similar effects are present in the indium 
antimonide crystal, a measurement was made of the total integrated 
intensity of the In resonance in the single crystal and the value obtained 
compared with the total intensity of th 1H resonance in an equal volume 
of a concentrated aqueous solution of ferric chloride, the latter being used 
to broaden the proton resonance line and to shorten the relaxation time. 
In this experiment an electromagnet was used, the magnetic field being 
adjusted so that the !!*In and +H resonances occurred at the same frequency, 
and a calibrator was used to monitor the spectrometer sensitivity. Since 
the resonant frequency is the same in each case, the total integrated 
intensities should be proportional to y2VJ(£+1), where y and J are res- 
pectively the gyromagnetic ratio and spin of the appropriate nucleus, 
N is the total number of nuclei in the specimen, and quadrupolar effects 
are neglected. In this way the intensity of the °In resonance was found 
to be 0:8 + 0-1 times the intensity to be expected on the assumption that 
all the satellites contribute to the observed absorption. (There can be 
no quadrupolar effects in the case of the proton resonance, since nuclei 
with spin 4 have no quadrupole moment.) Since for /=9/2 the central 
component contributes only 25/165 of the total intensity, almost the entire 
contribution from the satellites was observed. It follows, therefore, 
that any quadrupolar broadening effects arising from dislocations are 
far weaker in indium antimonide than in the alkali halides studied by 
Watkins and Pound. Observations of etch pits in indium antimonide 
single crystals similar to the one used in this experiment show that the 
dislocation density is of the order of 10? em~? (Allen 1957), which is of the 
same order as that found in good crystals of silicon and germanium, whereas 
Watkins and Pound assumed densities of 108 cm~® to explain the alkali 
halide results. Since one of the causes of dislocations is the stress asso- 
ciated with thermal gradients after solidification, the low dislocation 
density in the former group of materials as compared with the alkali 
halides may be partly due to their higher thermal conductivity ; another 
contributing factor may be the higher energy necessary to form a dislocation, 
arising from the directional nature of the bonds in the predominantly 
covalently bound semiconductors. 

The results of the present experiments differed from those reported by 
Shulman eé al. (1955) in two respects. Firstly, the In resonance was 
found to be weaker and broader when the material was crushed to a powder. 
This effect is attributed to quadrupolar broadening resulting from the 
introduction of dislocations during the crushing process, and will be dis- 
cussed more fully in a future paper. Secondly, the resonance in the single 
crystal appears to be about 10% wider at liquid oxygen temperature than 
at room temperature, as can be seen from fig. 1, which shows the high 
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frequency halves of the derivative curves plotted directly from the recorder 
chart. The reason for this is not clear. There can be no ‘ motional 
narrowing ° of the line at room temperature, since self diffusion in indium 
antimonide is known to be completely negligible at 290°K (Eisen and 
Birchenall 1957), and in any case an A,,I,.1, type of interaction does not 
average to zero. The theory of indirect exchange coupling as applied to 
insulators or semiconductors (Bloembergen and Rowland 1955, Shulman 
et al. 1955) shows that the constant 4,, depends on the energy gap, H 700: 
the distance between neighbouring nuclei, Rj, and on the quantity 


<[Ars/ ae = (1678 yB/3)"ggo Ug (try) p* (ry )ep (Pe) eye * (Fe) ee 


where f and fy are the Bohr and nuclear magnetons respectively, g; and r; 
the g-factor and position vector of the ith nucleus, and w;,(r) and w,/(r) 
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High frequency halves of 1°In resonance in intrinsic and Te-doped InSb. 
The curves show the derivative of the absorption at a frequency 
vot Av, where vy is the centre frequency. The high and low frequency 
halves were identical for the doped sample. 


the periodic parts of the Bloch wavefunctions for electrons with momentum 
vectors k and k’. All three quantities depend on temperature since they 
are functions of the lattice parameter. H, is known to be about 30% 
larger at 90°K than at room temperature for indium antimonide (Roberts 
and Quarrington 1955), but, as Shulman ef al. (1955) have pointed out, 
Aj, is not very sensitive to H,; moreover, Aj, should decrease as u y in- 
creases, giving a narrower line at 90°K. The quantity Aj. 18 reminiscent 
of the factor <[x5(0)|?), which occurs in the expression for the Knight 
shift (Townes et al. 1950), and which is known to increase as the temperature 
decreases. McGarvey and Gutowsky (1953) found their results on the 
temperature variation of the Knight shift in gallium and caesium to be 
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consistent with <[/(0)|?),0cv,-!, where v, is the atomic volume, and if 
we tentatively suppose that A,, shows the same behaviour, then the 
product R,~*¢|A,.|2),,, to which Aj, is proportional, should vary as 
v, 98, This is in the right direction to explain the broadening of the 
line at 90°K, but is unlikely to lead to an increase in A, of more than 1% 
orso. For the present, the origin of this broadening remains obscure. 


§ 4, FREQUENCY OF RESONANCE 

It is well known that nuclear resonance frequencies in metals are higher 
by some tenths of 1°% than the resonance frequencies of the same nuclei 
in insulators. This effect, known as the Knight shift (Knight 1949), 
is due to the Pauli paramagnetism of the conduction electrons, which 
increases the effective magnetic field acting on a nucleus. Only those 
electrons near the Fermi surface contribute to the Pauli susceptibility 
(yp), and this makes y,, proportional to 71? for a degenerate electron gas, 
n being the electron density. In view of this, it is just possible that a 
small Knight shift might be detectable in a highly extrinsic semiconductor 
with about 10!° carriers cm~*. Careful measurements of the 1°In re- 
sonance frequency were therefore made in n-type indium antimonide with 
an electron concentration of 10!® cm-%. This material is highly degenerate, 
having a Fermi energy about 30/7 above the bottom of the conduction 
band at room temperature. No deviation was found from the resonance 
frequency in intrinsic material, within an experimental accuracy of 1 in 10. 
The explanation of this is probably that the Pauli susceptibility is very 
low on account of the small effective mass of the electrons (m* =0-02m,), 
since y,, is proportional to the density of states at the Fermi surface, which 
is proportional to m* for a given electron density. 

The situation is somewhat different in p-type indium antimonide. Here 
the larger effective mass of the holes (m*=0-2m,) makes y, considerably 
larger than in the n-type case for the same carrier concentration. The 
experiment was therefore repeated with p-type material with 10!* holes 
em~*, but again no Knight shift was detected. To interpret this result 
it is necessary to consider the theoretical expression for the fractional 
Knight shift, which is Av/v = (87/3)y,,0,<(|4(0)|?)p. Here v, is the atomic 
volume, and (0) the wave function of a hole at the nucleus, % being nor- 
malized to unity within an atomic volume. The expression ¢[/(0)|?),, 
which will be designated by P,, represents an average taken over the 
Fermi surface. If we calculate y, using the Pauli formula modified 
for the case of intermediate degeneracy, viz. y,=nf2F’/kT'F, where n is 
the carrier density, 8 the Bohr magneton, and F = F,).(¢/k7') is the usual 
FermiDirac integral, we obtain the value y,=9x10-*e.m.u. em-. 
From measurements of the Knight shift in metals, it is usually found that 
P,, is within about 50% of the value of |y(0)|? relating to the valence elect- 
rons in a free atom. If we assume the hole wave function to have 
s-character, and substitute for P,, the value of [(0)|? appropriate to the 
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5s electron in the In?* ion, which is known to be 1:5 x 1026 em-3 from 
hyperfine structure determinations (Campbell and Davis 1939), we find 
Av/v=4x 10-4, 

which is four times the probable error of our measurements. In this 
calculation we have taken v, to be one half the volume per In-Sb pair. 
Ifthe value adopted for y,, be assumed correct, the absence of an observable 
Knight shift implies that the value used for P, is too large. This is not 
unexpected, since in the case of germanium and silicon it seems well 
established that the hole wave functions are predominantly p-like, and, 
bearing in mind that (0) vanishes for a p wave function, we can explain 
the small value of P;, by assuming that the same is true of indium anti- 
monide. There is another factor which may depress the value of P,, 
namely, the suggestion by Slater and Koster (1954) that the valence band 
is predominantly built up of antimony 5s and 5p orbitals, so that the hole 
wave functions tend to be concentrated around Sb atoms. However, 
assuming that the valence electron distribution is such as to give appro- 
ximate charge neutrality around each atom, this effect cannot be large. 
The result of this experiment can be stated more specifically by saying that, 
assuming the value of 9 x 10-°e.m.u. cm™ for y,,, the value of ¢[4(0)|?), 
for holes in indium antimonide must be less than 4 x 1025 cm~ at the position 
of an indium nucleus. 


§5. Tue Errect oF [ontzep IMPURITIES 


In the course of the experiments mentioned in § 4, it was found that 
the intensity of the 4*In nuclear resonance in InSb, as represented by the 
maximum value of dy"/dv, was very sensitive to the presence of donor or 
acceptor atoms. This is shown in fig. 2, which shows the maximum value 
of the derivative of the room temperature absorption plotted against 
tellurium concentration. For concentrations of 1018 atoms cm-%, the 
effect of crushing mentioned in § 2 is completely swamped by the effect 
of the impurity atoms. This decrease in intensity is accompanied by an 
increase in width, as is evident from fig. 1, which shows the line shape for 
the specimen with 9 x 1017 Te atoms cm~? at 90°K; this has much more 
pronounced wings than the single crystal curves and approximates more 
closely to a Lorentzian than to a Gaussian line shape. No asymmetry 
could be detected in the absorption lines of any of the doped samples, 
since the additional broadening brought about by the impurities tends 
to reduce the asymmetry which is already very small in the case of the 
single crystal. In fig. 1 the ordinates are scaled to make each curve have 
the same maximum value, which helps to emphasize the differences in 
shape. With the exception of the single crystal, all samples showed the 
same line shape and same relative intensities at 90°K as at 290°K ; provided 
saturation effects were avoided, the increase in the absorption intensity 
at the lower temperature corresponded to the increase in the Boltzmann 
factor which determines the difference in the populations of adjacent 


Zeeman levels. 


552 E. H. Rhoderick on the 


The various ways in which the presence of impurities could conceivably 
broaden the nuclear resonance line can be divided into two categories. 
Firstly, it is possible that the spin-lattice relaxation time (7',) may be 
shortened in some way and so give rise to ‘ lifetime broadening’. Such a 
mechanism was invoked by Reif (1955) to explain the effect of Cd** ions 
in broadening the Br resonance in silver bromide. Reif assumes that 
the relaxation mechanism arises from quadrupolar coupling of the Br 
nuclei to the electric fields associated with rapidly diffusing Agt vacancies. 
Other processes by which 7’, may be shortened are the Heitler—Teller 
mechanism associated with the extrinsic conduction electrons, in which a 


Fig. 2 
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Variation with Te concentration of derivative of 1>5In resonance in InSb 
(290°x). The value relating to the pure single crystal is taken as unity. 
Also shown are results for two particle sizes of the pure powdered 
material. 


nucleus and a conduction electron mutually flip their spins, and its quadru- 
polar analogue (Mitchell 1957). All these processes should be negligible 
in indium antimonide at room temperature, even with 10! conduction 
electrons cm~*. To verify that this is so, a brief study of the relaxation 
time in pure and impure samples was made, using the saturation method. 
Figure 3 shows relative values of (dy"/dv),,,. plotted against the radio- 
frequency voltage across the coil for the In resonance in a pure sample 
at both 290°K and 90°K. The spin-lattice relaxation time for this sample, 
deduced from the saturation curves, is of the order of 10-2 see at 290°K 
and 0-16 sec at 90°K. This comparatively short value of 7’, at room tem- 
perature suggests that the relaxation mechanism probably involves 
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quadrupolar coupling of the nuclei to the lattice vibrations (Pound 1950). 
Also shown are experimental points for a doped sample containing 10!° Te 
atoms cm~* at 90°K, the signal-to-noise ratio for this sample being too 
poor for saturation measurements to be made at room temperature. 
Although the scatter of the points relating to the doped sample is con- 
siderable, there is clearly no significant difference between the low tempera- 
ture relaxation times in the pure and impure samples. Furthermore, 
arelaxation time of 10~ sec or more is much too long to explain the observed 
broadening. 

The second class of broadening mechanisms comprises those which modify 
the field at the position of the nuclei. One obvious explanation is that 
the tellurium may give rise to paramagnetic defects which produce strong 
local magnetic fields. However, the tellurium atoms enter the InSb 


Fig. 3 
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Saturation curves of °In resonance in pure and Te-doped InSb. 


lattice as Te* ions, since there is one conduction electron per donor, and 
there is considerable evidence that they take the place of Sb atoms, the 
surplus Sb atoms segregating at grain boundaries. Since a Te* ion is 
isoelectronic with the Sb atom it replaces, one would expect no unpaired 
electrons and therefore no paramagnetic centres. To check this, a search 
for paramagnetic resonance in a sample with 10'® Te atoms cm? was kindly 
made by Dr. E. E. Schneider, of King’s College, Newcastle, using a high 
sensitivity spectrometer. No resonance was observed, and from the 
known sensitivity of the spectrometer it is possible to say that, if paramag- 
netic defects were present, the paramagnetic resonance line width was 
more than 500 gauss. Another possible explanation is that the conduction 
electrons freed from the Te atoms may increase the constant A,, of § 3 
and so broaden the resonance. In the case of indium antimonide with 101 
electrons cm-* the degeneracy is so high that we may apply the theory 
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developed by Ruderman and Kittel (1954) for the case of metals. Accord- 
ing to their calculation the contribution of the conduction electrons 
should be negligible in this case; moreover the pronounced departure 
from Gaussian shape suggests that this mechanism is not the one responsible 
for the observed broadening. A third possible broadening mechanism 
which belongs to this category differs from the preceding ones in that it 
involves an electric rather than a magnetic interaction, namely, the 
interaction of the nuclear quadrupole moment with the Coulomb field of 
the Te+ ion.. Since both the magnitude of grad-E and the direction of its 
symmetry axis vary from nucleus to nucleus, the satellites become smeared 
into a continuous band, as mentioned in § 3 in connection with the possible 
role of dislocations. This is believed to be the mechanism responsible 
for the observed broadening. 

According to the explanation given above, the. maximum effect that 
quadrupole interactions can have is the complete obliteration of the 
satellites, the central component remaining unaltered. This is true as 
long as the quadrupolar perturbation is treated only to first order ; if the 
second order perturbation is considered, it is found that the central com- 
ponent is also broadened. For J = 9/2, the central component contributes 
25/165 = 0-15 of the total intensity in an unstrained crystal, and since fig. 2 
shows that with 9 x 10!8 Te atoms cm~? the maximum value of dy”/dv falls 
to 7% of the value for the pure material, it is evident that second order 
effects must be contributing. Figure 2 can now be interpreted as follows. 
The initial rapid fall in intensity with Te concentration is due to first 
order broadening of the satellites, while the subsequent slow decrease, 
which extrapolates smoothly to the value 0-15 at zero concentration, is 
due to second order broadening of the central component. According to 
theory, the second order perturbation should lead to an asymmetric line, 
but there was no clear evidence of asymmetry in the line shape of the 
most heavily doped sample, although the signal-to-noise ratio was too 
poor to allow any definite conclusion to be drawn. The maximum value 
of the derivative is apparently a more sensitive function of impurity 
content than is the asymmetry. 

This explanation in terms of quadrupole interactions implies that ionized 
acceptors should have the same effect as donors, since reversing the sign of 
grad‘E merely makes a particular satellite change places with its mirror 
image. The table gives a comparison between p-type samples obtained 
by doping with zine and the n-type samples doped with tellurium. The 
slightly larger intensity in the p-type specimens can be explained by the 
partial screening of the charged ions by the carriers. The theory of this 
screening has been given by Mott and Jones (1936) for metals and extended 
to the case of semiconductors by Mansfield (1956). According to Mansfield, 
the potential due to the screened Coulomb field of an ion is 


V(r) =(e/er) exp (—4qr), 


where g? = 167e?(2m*)3?4/(kT') F's /.(y)/eh3. Here m* is the effective mass 
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of the carriers, ¢ the dielectric constant, and F'y9(n) = F',.(C/kT) is the usual 
Fermi—Dirac integral. The field gradient at a vais r from the ion now 
becomes (see Cohen and Reif 1957) 

GV _ 2 oe irl 1 g2p2h 

Ree ary xp (—ar){1+qr+4q?r"}. ates, ore L) 
With 10" carriers cm~* in InSb, g has the value 1-5 x 106 em! for electrons 
(m*=0-02m,) and 4:2x10%cem for holes (m*=0-2m,). With 10! 
impurity atoms em~*, the maximum distance which any nucleus can be 
from an ionized atom is about 50 A. Using this value of 7, the term in 
eqn. (1) arising from the screening amounts to 0-92 for electrons and 
0:56 for holes. The effect of screening is therefore almost negligible for 
n-type material with 101° electrons em~*, but in the case of the p-type 
specimen the screening reduces the field gradient by an amount which is 
of the right order to explain the difference between the resonance intensities 
in the two cases. For 10!§ carriers cm~%, ¢ is 1-95 x 10° cm~! for holes 


Sample Relative 
intensity 
Intrinsic (single crystal) 1-00 
n-type (n ,=101§ cm-%) 0-23 + 0-01 
p-type (7, = 1018 cm-’) 0:25+0-01 
n-type (7 ,= 101° cm?) 0-07, + 0-005 
p-type Be sae er) 0-135 + 0-008 


0-07, + 0-005 


faa = 10? em 
compensated if 


ree, 18 em-3 
N g—N,=2 x 1018 cm 


and 7x 10° cm-! for electrons, so the effect of screening is very small 
in both cases. For a partially compensated specimen with 

Nat, = 10 em and ng—n,=2 x 101% cm 
the amplitude of the resonance is indistinguishable from that with 
nq=101® em-, since the screening is negligible and the number of ions is 
the same in both cases. 

While this explanation in terms of the interaction of the quadrupole 
moment of the nucleus with the screened Coulomb field of the ions gives a 
good qualitative interpretation of the data, the quantitative agreement is 
more speculative. First order perturbation theory shows that the change 
in frequency of the m<—>m— 1 satellite due to a charge e at a point (7, 0) 
relative to the nucleus is given by 
3eQ(2m—1) a?V 
41(2I1—1)h O22 
_ 3eQ(2m—1)_ Be( (3cos*@—1) _ Pah ay Sea?) 
7 4i(2F= 1h ee 


Av= 
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where J is the nuclear spin, Q its quadrupole moment, and the direction 
of the magnetic field is taken as 2 axis (Pound 1950). Within a solid, the 
field gradient will be modified due to polarization effects, and we allow for 
this by incorporating in the right-hand side of eqn. (2) an empirical multi- 
plicative constant B (strictly speaking a tensor). To calculate the shape 
of the broadened line it would be necessary to find the distribution of 
values of Av for different nuclei resulting from a random distribution of 
impurity atoms. This is an exceedingly complex problem, and it is simpler 
to calculate the second moment of the resonance linet. From eqn. (2) 
it is possible to show that, neglecting lattice structure and assuming a 
perfectly random distribution of impurities, the broadening of the resonance 
absorption due to first order quadrupole interactions with charged impuri- 
ties gives rise to a second moment 

FIER ae Mh Lee eh : 

(Avs Gla teh phe ae 


where N is the impurity concentration, and the effect on a given nucleus of 
impurities lying within a distance r, of that nucleus is neglected{. The 
use of this cut-off radius is necessary to prevent the second moment diver- 
ging ; its introduction is also physically plausible since the lattice structure 
and ionic size set a lower limit to the distance between an impurity atom 
and anucleus. In practice, however, 7, is determined by other considera- 
tions. As the line shape in the doped sample approaches a Lorentzian 
form, the second moment converges more and more slowly since the 
contribution from the wings becomes increasingly important. As the 
wings eventually become lost in noise, it is impossible to measure the 
second moment of the entire curve, and one is driven to the expedient of a 
cut-off frequency (Av). The simplest procedure is therefore to determine 
(Av), experimentally as that frequency beyond which the signal-to-noise 
ratio is too poor for accurate measurement, and then to find r, by substitu- 
ting this value of (Av), in eqn. (2) with the angular factor set equal to unity. 
(This slightly overestimates r, since, on account of the angular term, some 
impurities less than r, from a given nucleus will still give Av<(Ap),.) 


+ An exhaustive treatment of the effect on the line shape of quadrupolar 
broadening by point defects is given by Cohen and Reif (1957), who show that 
for small concentrations of defects the contribution to the line shape from a 
given pair of satellites becomes Lorentzian. Since the central component is 
unaltered, the line shape in this particular case should be a superposition of 
one Gaussian and four different Lorentzian curves. The author is indebted 
to Professor Cohen for letting him see the manuscript of this paper in advance 
of publication. 

{ This expression can be obtained using an approximate method due to 
Bloembergen (1955) if one writes 6 for Bloembergen’s (1+ |y..|) and corrects 
a numerical error in his expression for (Av)?. It can be obtained more 


rigorously by a method similar to that used by Reif (1955) to treat second 
order broadening. 
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In this way we find 


a; (A  47N e2QB(2m—1) 
| P| eis 7 “Tar—1yp| \AYe Gobet nti) 


This equation applies only to the contribution from the m——>m — 1 satellite, 
which has an asymmetric line shape. If we consider the contributions 
from the m<—+m—1 and —m<—> —(m— 1) satellites together, the resulting 
line shape is symmetric and has a second moment given by eqn. (4). We 


may use this expression for (Av)? to include all satellites by taking + 5/2 
as a mean value form. Figure 1 shows that the line shape for the sample 
with 9 x 10!” Te atoms cm~* approximates quite closely to the derivative 
of the Lorentzian curve g(Av)oc[(9-6 ke/s)?+(Av)?]1 up to Av=30 ke/s. 
Assuming this analytical form for the curve, we find| (Av)? |)’“"" = 137 (ke/s)2. 
The Gaussian line shape in the pure material has (Av)? = 29(ke/s)2, so the 


quadrupolar contribution to [ (Av)? }° hte calla (ke/s)?. Substituting 
this in eqn. (4) with (Av),= 30 ke/s and Q = 1-16 x 10-4 cm2, we find B = 350. 

To interpret this large value of 8 we write it in the form B = Bo, following 
Reif (1955), where 6. and f£, have the following significance. Let us 
imagine a spherical hole in the InSb crystal formed by removing one indium 
atom. We consider the field gradient at the centre of this sphere and 
suppose it to be given by 28.e/r?, so that B. describes the effect of the 
polarization of the surrounding medium. We then replace the indium 
atom and consider the effect of the distortion of its own electron cloud 
on the field gradient at the nucleus, supposing this to be equivalent to 
multiplication by the factor By. This is admittedly a rather artificial 
procedure, since indium antimonide is a predominantly covalently bound 
compound, and the concept of a sphere with a complete atom inside and the 
others outside lacks realism, but it should at least serve as a rough model 
sufficiently exact to assess orders of magnitude. If the polarization of 
the surrounding medium produces an internal field of the Lorentz type, so 
that H.g=H+4nP/3, then B.=(e+2)/8e~3, since e=16. On the 
other hand, if the internal field is of the Drude type, as in the case of ionic 
crystals where there is considerable overlap of the ions (Mott and Gurney 
1948), then H,,=H and B.=1/e. The correct value of f. will lie some- 
where between these two extremes, but we will take B.=4 since this 
demands a lower value of By, namely, By ~ 10007. 

It is well known that in diatomic polar molecules the interaction of the 
quadrupole moment of either nucleus with the electric field is much greater 
than that arising from a field gradient 2e/r), where rg is the internuclear 
distance. This has been attributed to the distortion of the closed electron 
shells, which under certain conditions may have an ‘ anti-shielding ’ effect. 


+ By using a more refined argument, Cohen and Reif (1957) arrive at the 
result B.=(2¢+3)/5e, which is slightly higher than the value adopted here. 
In view of the roughness of the approximation, it makes little difference which 
of the two expressions is used, 
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The anti-shielding factor, which we may indentify with f), has been 
calculated for ions having the inert gas configuration by Foley et al. (1954). 
Interpolating between their results for Rb* and Cs*, we may estimate 
By 70 for Agt, which is isoelectronic with In**. The value of By for 
In*+ will probably be lower because, having a larger nuclear charge, In** is 
more tightly bound. It seems, therefore, that the distortion of the indium 
core electrons is much too small to explain the large value of By required 
to fit the experiments. This suggests that the effect may be due to distor- 
tion of the valence electron distribution, and that this should be so is not 
surprising if one considers how the distortion arises. The electron dis- 
tribution around the In nucleus can be distorted only by mixing with the 
unperturbed wave function new wave functions corresponding to excited 
states. The extent to which these perturbing wave functions are admixed 
depends, among other things, on the energy difference between the initial 
state and the excited states involved. In the case of the In** core, the 
lowest excited state above the full N-shell is the 5s state, which is more 
than 10 ev higher. However, in the case of the In atom in InSb, the 
energy difference between the top of the valence band and the lowest 
unoccupied state is only of the order of the forbidden gap, i.e. about 
0-18 ev. We might therefore expect the valence electrons to be much more 
easily distorted than the core, and it is significant that the quadrupolar 
part of the hyperfine interaction in the (5s)?5p ?P3,. state of the free indium 
atom corresponds to a value of (3eQ0?V /dz?)/41(27 — 1)h of 37 me/s (Mann 
and Kusch 1950), so that we need a distortion equivalent to the addition 
of only 0-1°%, of a 5p orbital to give a first-order broadening of about 30 ke/s 
in solid indium antimonide. 

If these ideas are correct, we should attribute the very large value of B, 
in indium antimonide, as compared with polar molecules and ionic crystals, 
to the smallness of the forbidden gap. This implies that there should be 
some correlation between the value of 8) and the inverse of the gap width, 
rather like the relation between « and HL’, suggested by Briggs et al. (1954). 
Some preliminary experiments carried out in this laboratory on indium 
arsenide, which has an energy gap of 0-34 ev, show that the 15In resonance 
is affected by impurities in much the same way as in indium antimonide, 
but that the value of 8) is smaller, which is consistent with the above 
hypothesis. 


§ 6. THe Errecr or NeutTRAL IMpurirres 

The decrease in intensity of the ">In nuclear resonance on doping with 
tellurium or zinc has been attributed to the effect of the Coulomb field 
of the charged ions. An obvious method of verifying this hypothesis 
would be to neutralize the ions by lowering the temperature. Unfortuna- 
tely this is not practicable in indium antimonide, which has zero activation 
energy for ionization of donors, and an activation energy for acceptors 
which is zero except at liquid helium temperatures. An alternative 
approach is to dope the material with impurities belonging to the same 
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group as one of the constituents, which would be electrically neutral. 
Through the kind cooperation of Dr. J. C. Woolley, of Nottingham Univer- 
sity, the In resonance has been observed in a series of mixed crystals 
of indium antimonide diluted with gallium antimonide. These compounds 
are isomorphous, and combine to form a single phase over the entire range 
of compositions, providing the annealing time is long enough (Woolley et al. 
1956). Figure 4 shows that the effect of a small concentration of GaSb 
on the '*In resonance at room temperature is similar to that of Te atoms, 
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Ga CONCENTRATION (ATOMS cm>) 


Variation with Ga concentration of derivative of !!5In resonance in 
InSb-GaSb mixed crystal (290°K). 


except that for a given reduction in intensity the Ga concentration has 
to be about thirty times greater than that of the Te. This large difference 
confirms the view that the effect of the Te atoms arises almost entirely 
from the charge of the Te+ ion. The much smaller effect of the gallium 
can be attributed to strain in the InSb lattice in the neighbourhood of a 
Ga atom. The distance between nearest neighbours is 2-80 A in InSb 
and 2-65 & in GaSb, so if we assume that the length of the Ga-Sb bonds 
associated with an isolated Ga atom which is substituted for In in the InSb 
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lattice tends to be the same as the Ga—Sb distance in pure GaSb, then the 
InSb lattice around the Ga atom must be in a state of tension. This strain 
destroys the tetrahedral symmetry of the perfect lattice and gives rise to 
quadrupolar broadening. Although detailed studies of the line shapes 
in the GaSb-InSb mixed crystals have not been made, a rough order of 
magnitude interpretation of the results may be made as follows. If the 
material within a sphere of radius a, in an isotropic elastic medium is 
contracted to radius a, then the strain at a distance r from this defect is a 
pure radial one whose magnitude is given by ordinary elasticity theory as 


= Ar |r = [(a4— a2) [ay] x (a4/7)°. 
At some distance from a Ga atom, we may apply this continuum model 
and assume the radial strain to be given by this formula with 

4 — Ay = (2°80 — 2-65) A, 
and a, equal to the ‘ tetrahedral covalent radius’ of a Ga atom, which 
according to Pauling (1940) has the value 1-26 A. The relation between 
the field gradient and the strain at a given point in the lattice is a tensor 
one which has been studied experimentally for the case of InSb by Shulman 
et al. (1957). These authors find a relation between the diagonal com- 
ponents of the field gradient and stress tensors of the form 
0? V fox? a Cy [EXE ze x aA, oF AG) 

(and cyclic permutations), where x, y, and z lie along the cubic axes, X is 
the usual stress tensor, and Cj, is a ‘ gradient-elastic ° constant having 
the value + 2-4 x 10*e.s.u./dyne. We may write this in a form involving 
the strain rather than the stress as follows : 


OY [Ox? = Cig(Ci1— Cra) ene — (Cy gt Cee) | ae 


a 


Here c,, and c,, are the usual elastic constants which have been measured 
for InSb by Potter (1956), and e is the strain tensor. For a pure radial 
strain of the form Ar/r=k/r?, one can show e,,,=k(r?—32?)/r®, ete. If 
we suppose the magnetic field to coincide with one of the cube axes and call 
this the z-axis, then we finally obtain 


OV /dz? = 3k0,,(€1, — Cy2)(1— 3 cos? O)/2r3 . 2 . . (6) 


This has the same dependence on r and @ as the field gradient due to a point 
charge, Be(3 cos? @—1)/r, which was used in eqn. (2). We may therefore 
use eqn. (4) for the case of strain by substituting 3kC,,(c,, —¢,,)/2 for Be. 
Using Potter’s values for c,, and c,, and the above mentioned values of k 
and C,, we find 3hC};(¢,; — ¢y)/2= + 2-75 x 10-® e.s.u., which is equivalent 
to a value of B of 6. In the absence of a detailed knowledge of the line 
shape, the second moment cannot be estimated, but we may note that 
eqn. (4) involves 8 and the concentration of defects N only through the 
product BN. We have seen from the results with the Te-doped samples 
that the charge effect demands a value of 8 of 350. Since the strain effect 
is equivalent to a value of 8 of 6, we might expect that, if the line shapes 
are roughly the same in two the cases, we should require a concentration 


SS 
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of Ga atoms about 60 times that of the Te atoms for the same effect on the 
nuclear resonance intensity. This compares reasonably well with the 
factor of about 30 which is obtained by comparing figs. 2 and 4, The 
assumption that the magnetic field is parallel to one of the cubic axes implies 
the use of a single crystal; in fact, the material was polycrystalline, but it 
is unlikely that this will materially affect the broad conclusions of this 
analysis. 

The conclusion that the broadening produced by Ga atoms is attributable 
to strain has relevance to the interpretation of the tellurium experiments, 
since the Te atoms must also produce some strain in the lattice, and it is 
important to know the relative magnitude of the strain and charge effects. 
In In, Tes, which is a * deficit ’ semiconductor with the zinc-blende structure 
having one third of the In sites vacant (Welker 1954), the In—Te distance 
is 2-66 A, which is almost indistinguishable from the Ga-Sb distance of 
2-65 A in GaSb. If we assume that the In—-Te bond length in Te-doped 
InSb is approximately the same as the distance between nearest neighbours 
in In,Te, then the strain caused by the substitution of a Te+ ion for an Sb 
atom in InSb should have a broadening effect comparable with that ofa 
Ga atom, which we have seen to be about two orders of magnitude iess 
than the observed effect of the Te atoms. We can therefore say fairly 
confidently that the effect of the Te+ (and Zn) ions must be due to their 
charge. 


§ 7, CONCLUSION 

The aims of this research were twofold. The object of the search for a 
possible Knight shift was to investigate the s- or p-character of the electron 
and hole wave functions. The results are consistent with the hole wave 
functions being predominantly p-like, but in the case of n-type samples, 
the results are ambiguous because of the small effective mass of the electrons. 

The motive for studying the dependence of nuclear resonance intensity 
on impurity concentration was the hope that the results might be inter- 
pretable in the simple way which Bloembergen and Rowland (1953) 
used to explain their results on the Cu resonance in «-brass. In the case 
of a metal, the Coulomb field of an impurity is rapidly damped out by the 
screening of the conduction electrons, and the crystalline field is only modi- 
fied in the immediate vicinity of the impurity atom. Bloembergen and 
Rowland found that their results could be explained by assuming the 
existence of a critical radius around a Cu nucleus such that, if there is a 
Zn atom inside this radius, the quadrupolar interaction is so large as to 
obliterate the contribution of that nucleus, whereas, if there is no Zn atom 
inside this radius, the resonance of the Cu nucleus is unmodified. In this 
way they showed that the extra electron on a Zn atom is shared with a 
Cu atom in the 18 nearest and next-nearest neighbour positions. A 
similar result is possible in insulators if the interaction is so weak 
that it is appreciable only in the immediate vicinity of an impurity, 
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in which case the 7-3 variation (or 7~* in the case of the second order inter- 
action) may allow discrimination between nearest and next-nearest neigh - 
bour positions (Kawamura et al. 1956). The initial hope was that a. 
similar situation might arise in InSb, i.e. that the quadrupolar interaction 
would be appreciable only for nuclei in the immediate vicinity of a defect, 
and that in this way it might be possible to prove, for instance, whether 
the impurities enter the lattice substitutionally or not. Unfortunately, 
the large value of 8 in InSb makes the interaction appreciable even at points 
more than 50 A from the impurity, and at these large distances the lattice 
structure ceases to play a significant part. Paradoxically, nuclear 
resonance in InSb is too sensitive to the presence of imperfections to be. 
helpful in studying the immediate vicinity of a defect, and the usefulness 
of nuclear resonance in this respect tends to be restricted to those materials 
which have much smaller values of 8 and Q. With this in mind, it is 
proposed to make a brief investigation of the applicability of the technique 
to other intermetallic semiconductors, especially those having higher 
energy gaps. 
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ABSTRACT 


An analysis of Rhoderick’s experiments on the nuclear quadrupole broaden- 
ing of the nuclear magnetic resonance of !1°In in impure indium antimonide is 
carried out which is based on a detailed theory of the line shape. It results 
in values of 580 for the multiplication factor and 1300 for the Sternheimer 
antishielding factor of In in InSb, 


Rhoderick (1958) has shown that addition of impurities to indium anti- 
monide causes quadrupolar broadening of the '°In nuclear magnetic 
resonance. In n-type (Te-doped) and p-type (Zn-doped) samples, the 
broadening arises through interaction of the nuclear quadrupole moment 
with the gradient of the electrostatic field of the ionized impurity atoms. 
The broadening is first-order at low concentrations of impurity. In 
particular, the line is approximately Lorentzian with a half-width of 
9-6 ke/s in a sample containing 9 x 10!7 atoms per cm?, a concentration 
at which screening effects are unimportant. From this result, Rhoderick 
has estimated that the multiplication factor 8 is approximately 350 for In 
in InSb. The estimate is based on an approximate calculation of the 
second-moment of the resonance line. We give here an independent 
estimate of 8 based on a detailed theory of the line shape (Cohen and 
Reif 1957). 

If the resonance line were sharp in the pure crystal, the shape of the 
component arising from transitions m<—_>m — 1, m 44, in the impure crystal 
would be Lorentzian with half-width a 

if 20) 
Nae oi 
34/3 1(2L—1)h 
where N is the number of impurity atoms per unit volume. The resonance 
line is actually broadened by indirect exchange in the pure crystal (Shulman 
et al. 1955) so that the (normalized) Bo per unit Bp interval is 


1 
Im(v )= =| aye ry re rel . . . . . (2) 


in the impure crystal if it is f,,(v) in the pure crystal. The frequency v 
is measured from the centre of the line. If I, is larger than the line width 
in the pure crystal, eqn. (2) becomes 


J Dn 2 a 1 in 
a7 tye +4o m Tae Dace alee 4 6 6 (3) 


Saree eee 
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nee : ; 
where o,,” is the second moment of /,,(v). At frequencies such that v2 is 
greater than [’,,?, g,,(v) becomes 


feel 2 
Im(v) = Pn 1— Pgt= 80m! scles cavity Iie . *, i ‘ (4): 
7 Vv 


In the wings of the resonance line, where the central component does not 
contribute, the line is a superposition of these components, each weighted 
according to its relative intensity J(1+1)—m(m—1). The total (nor- 
malized) intensity per unit frequency interval is 

g(v) = 2 Pndnl”), Ma > J 2 by . : " (5) 
mAs 


where p,, is the weight of the mth component and g,,(v) is given by eqn. (4). 
Combining eqns. (4) and (5) we obtain 


Lena 
g(v)= ae ¥ mi > (DON RES| Re =e, 30°), niece er. (6) 
where I is obtained by replacing the |2m—1] in I, by [2m — 1], 
|2m —1]= 2 Pml2m— 1]. a ae) 
MAS 


The spin of 1°In is 9/2, for which |2m—1]is 4:28. This corresponds to an 
“average ° value of m of 2-64, which is close to the value of 5/2 assumed 
by Rhoderick. 

The wings of the line are indistinguishable from those of a single Lorent- 
zian of half-width I if we neglect the second term in (6). The value of 
the half-width of a Lorentzian giving a good overall fit to the line shapes 
for the Te-doped (V=9 x 10!”) sample is 9-6 ke/s. Inspection of Rhode- 
rick’s fig. 1 suggests that [is actually somewhat less than 9-6 ke/s. We 
therefore obtain a value of 580, or perhaps somewhat less, for 8. The 
correction term in (6) may be estimated by using the approximation 


TT 2 Ge 
av)= —(1- —*) eee ok 3) 


Vv 


where o? is the second moment in the pure crystal. The correction is 
small when eqn. (6) is valid. 

The multiplication factor can be decomposed by means of the approxi- 
mation 


ley lero) oe emer ear ge (9) 


where 1+y is the Sternheimer antishielding factor (Foley e¢ al. 1954, 
Cohen and Reif 1957), the remaining factor is the polarization correction 
(Cohen and Reif 1957), and « is the dielectric constant. The polarization 
correction is 0-437 for InSb for which «=16. Rhoderick’s estimate of f 
yields 1 + ~ 800; the present estimates yields 1+yZ1300. This difference 
in 1+y is probably within the error of eqn. (9), which is based on a highly 
idealized model more suitable to ionic crystals than to semiconductors. 
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ABSTRACT 


The giant resonances of nuclear photodisintegration show a strong correla- 
tion between their width and the distance from doubly closed shells. In 
particular the width becomes very great in the rare earth region of large 
permanent nuclear deformations. It is shown, with the aid of the Nilsson 
level schemes for a deformed harmonic oscillator potential, that the inde- 
pendent particle model of the giant resonance is capable of providing a 
quantitative explanation of this broadening. 


$ 1. IyTRoDUCTION 


Two apparently conflicting models seek to explain the phenomena asso- 
ciated with the giant resonance of nuclear photodisintegration. The 
first is the collective model (Goldhaber and Teller 1948, Jensen and Jensen 
1950, Danos 1952) in which the resonance is pictured due to a bulk 
oscillation of all the protons of the nucleus against all the neutrons. The 
second model ascribes the resonance to excitations of individual nucleons 
within the general framework of the shell and optical models (Wilkinson 
1954, 1956, Rand 1957). 

Many of the predictions of these opposed models are surprisingly similar : 
this is true of the energy of the resonance and its trend with atomic weight 
and of the strength of absorption in the resonance. The independent 
particle model has enjoyed certain specific successes, notably its satis- 
factory treatment of the anomalous emission of fast protons by heavy 
nuclei. This success stems from the fact that it is a much more detailed 
model from the point of view of nuclear structure and so may be immed- 
iately applied to any phenomenon in which the behaviour of single 
nucleons is of importance. It is not clear that the collective model could 
not be broken down so as to enable one to form a picture of the behaviour 
of the individual nucleons. Equally it is not clear how this could or 
should be done and it certainly seems at first sight that the conceptual 
basis of that model, totally different from that of the independent particle 
model, precludes our concentration of the energy of excitation on a single 
nucleon such as is necessary to explain the anomalous emission (see how- 
ever, Brown and Levinger 1958). We might therefore maintain that the 
phenomenon of proton emission argues in favour of the shell model and 
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against the collective model. Such a judgment would not necessarily be 
wrong but would be hasty. In fact Brink (1957) has shown that in one 
extreme case—the pure harmonic oscillator without damping forces—the 
models are formally identical. He argues that this surprising identity 
will persist in some form for real nuclei and so, by implication, that it is 
futile to attempt to distinguish between the models. Be this as it may 
the confrontation of the models continues and will do so until a full analysis 
of the relationship has been effected. 

An advantage enjoyed by the collective model is that it is always easy 
to make predictions with it because the motion envisaged is so simple and 
general. It is doubtful whether any of its important results were un- 
known to the third Lord Rayleigh. The shell model is relatively slow off 
the mark because a large amount of computation must precede the getting 
of any quantitative result. This is of course the price paid for our ability 
to make detailed predictions. 

This note describes a preliminary application of the shell model to an 
aspect of the experimental situation—the width of the giant resonances 
—which has recently acquired some topical interest following the first 
careful giant resonance studies among the rare earths (Petree et al. 1957). 


§ 2. THE RESONANCE WIDTH 


For spherical nuclei the collective model makes no simple prediction 
about the width of the giant resonance because the width is due to 
‘frictional’ forces which are not included in the model. The shell model 
can make some prediction because the two chief components of the width 
for doubly magic nuclei, the spread in energy between the several inde- 
pendent particle transition energies and the smearing in position of each 
due to the imaginary part of the optical model potential, are both fairly 
well understood. If the latter effect is set at 3Mev full width at half 
maximum the total predicted width for a doubly-magic nucleus is 4-5 Mev 
(see later for an example) which accords well with experiment. 

It has recently been found (loc. cit.) that in the region of the rare earths 
the resonance width becomes very much greater and may even exceed 
10 mev. This is the familiar region where nuclei are permanently 
deformed into prolate ellipsoids of appreciable eccentricity as is revealed by 
the rotational structure of their low-lying level schemes and by their 
greatly enhanced probabilities for E2 transitions (see e.g. Bohr and 
Mottelson 1953). 

The collective model is able to make an immediate and dramatic pre- 
diction about the giant resonance in this region: it should split into two. 
Classically speaking if the neutrons and protons vibrate against each other 
along a minor axis they would be expected to do so at a higher frequency 
than along the major axis. This splitting, whose magnitude depends on 
the particular version of the collective model which is employed, will be 
seen at least as an increased width to the resonance and possibly even as a 
separation into two peaks. This has been treated by Okamoto (1958) and 
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Danos (1956, 1958) and provides a semi-quantitative account of the 
experimental data although the full width of the resonance is probably 
not accounted for. 

It is clear that qualitatively the same effect is expected from the inde- 
pendent particle model since in the deformed ellipsoidal potential the 
independent particle levels are split up in the manner first discussed in 
detail by Nilsson (1955) and so what in the spherical potential is a single 
transition now becomes several. A glance at the Nilsson diagram shows 
that the levels in the ellipsoidal potential become completely interlaced 
at quite modest values of the deformation and so it may be thought that no 
resonance structure will persist. This however is wrong because in the 
strongly deformed potential new selection rules arise (Strominger and 
Rasmussen 1957) which limit the possible transitions. 


§ 3. Tor DEFORMED POTENTIAL 


The only way to see what is happening to the independent particle 
model of the giant resonance is to recompute it for the distorted potential. 
This is a most laborious and uncertain operation and it is not pretended 
that the present results have much more than an illustrative significance. 
However a completely consistent calculation would imply a complete 
solution of the problem of the distorted nucleus and to this we cannot 
presently aspire. 

For the present purpose we have used Nilsson’s results. This is un- 
satisfactory in that they are based on a harmonic oscillator potential which 
is unrealistic in such heavy nuclei and it is well known that these level 
schemes do not give as good an account of the ground state spins in the 
rare earths as do more realistic potentials. However the Nilsson 
Hamiltonian contains a spin-orbit term and also a further arbitrary term 
proportional to /? which is made to caricature, for spherical nuclei, the level 
ordering prescribed by nature. The great advantage of the Nilsson scheme 
is that it can be made to apply to nuclei of any mass by a suitable adjust - 
ment of the fundamental isotropic oscillator spacing. We have therefore 
used it for this preliminary examination of the distorted nuclei. 

An immediate difficulty lies in the choice of the parameters of the model. 
We may adopt the simple plan of using the same parameters as are sug- 
gested in Nilsson’s paper and then finally taking that value of the funda- 
mental oscillator spacing that locates the giant resonance at the right 
absolute energy. In this way the predictions would be unambiguous. 
This procedure of course involves us in the use of considerably greater 
values of the oscillator spacing than were envisaged by Nilsson in 
recommending suitable values for the other parameters. This large oscilla- 
tor spacing is not in itself objectionable, being associated with the use of 
the reduced mass (velocity dependent potential) in the calculation of the 
level positions (Wilkinson 1956). It has in fact been shown by Frahn and 
Lemmer (1957) that a velocity dependence of the interaction for nucleons 
moving in a spheroidal harmonic oscillator potential gives a level scheme 
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that is approximately that of Nilsson. Their numerical estimates of the 
effective Nilsson parameters for a reduced nucleon mass of 0-5 times the free 
mass (Brueckner 1955) are close to those originally recommended by 
Nilsson. Their effective oscillator spacing for these parameters is 82A~!/3Mev 
which is itself close to the value needed to give agreement in absolute 
position between the experimental and calculated giant resonance energies. 
This work thus constitutes the justification for our simple procedure of 
using the Nilsson scheme with a final adjustment of energy scale on account 
of the velocity dependence of the potential to fit the absolute resonance 
energy. 

It is to be hoped that calculations such as those of Gottfried (1956) for 
a deformed finite square well will become available over a range of atomic 
weight and will use both a velocity dependent potential and a radius 
adjusted to reproduce the experimental density of nuclear matter. The 
final refinement of a diffuse edge to the nucleus would also be welcome. 


$4. THe NucLEAR DEFORMATION 


The deformation 6 used by Nilsson is related to the intrinsic quadrupole 
moment Qy by: 
Qy = 0-8Z.R75(1+ 38). 


Fig. 1 
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The nuclear deformation 5 as defined by Nilsson deduced from the quadrupole 
moments of even-even nuclei in the rare earth region. The full line 


aL the trend of 6 used for the present computation of the resonance 
widths. 
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The values of § deduced from this expression, and using R= 1-2A1/3 fermis, 
are shown in fig. 1. The Q) we have taken from a variety of sources for 
even—even nuclei only (see e.g. Alder et al. 1956). Values of 5 derived from 
nuclei of odd mass are in accord with those shown on the figure. The 
full line of the figure shows the trend of 8 with A that we have assumed in 
the subsequent computations. 


§ 5. THE CoMPUTATIONS 


A number of values of 4, Z throughout the range shown in fig. 1 were 
chosen, with their appropriate values of 5, for detailed computation. For 
each A, Z, 5 value the Nilsson diagram was filled with the appropriate 
number of neutrons and protons and the energies and strengths of the El 
transitions were calculated. In order to limit the complexity of the work 
only transitions of the parent type 1/—11/+1 and 2/’-+ 2l’+1 were con- 
sidered since they provide the overwhelming bulk of the total E1 absorption 
strength (Wilkinson 1956). 

The energies and strengths of all possible transitions, several dozens for 
each nucleus, were tabulated and then each transition was smeared into a 
Gaussian with a full width at half maximum of 3Mev} to represent the 
effect of the imaginary part of the optical model potential. The resulting 
transition strengths were then grouped at equal energy intervals and 
plotted as shown on figs. 2 to 5. These figures represent deformations 
5=0, 0-1, 0-2, and 0-3 for elements working backwards in mass from Pb 
to about Ho. The computed points are represented by the bold dots and 
the line has been put smoothly through them by eye. In order to fit the 
absolute energies of the experimental resonances we must set one of the 
arbitrary energy units of these figures equal to about 2 Mev. 

The widths are seen to increase with 6 in the way qualitatively expected. 
For 6=0 the assembly of transitions is remarkably symmetrical and has 
a width of 4:4mev for Pb in complete accord with experiment. As 6 
increases the resonance broadens and develops an asymmetry on the high 
energy side which, for 6=0-3, becomes an actual splitting of the resonance 
into two. This behaviour runs closely parallel to that predicted by the 
collective model and gives the impression that perhaps here again we have 
a phenomenon that cannot be used to distinguish between the models. 
There is a tendency for the transitions in the higher energy peak of fig. 5 
to come from the larger values of Q, the projection on to the nuclear sym- 
metry axis of the individual j of the nucleon in question, while the smaller 
Q values are found in the lower peak. This result is qualitatively to be 
expected and is confirmed by a glance at the Nilsson diagram. High 
values of Q correspond classically to orbits lying in the plane of the minor 


+ Doubtless this should depend on 6 and on the levels concerned but we have 
no data on which to assess such a dependence. It is not necessarily to be 
expected that the smearing will increase with increasing transition energy 
since a high transition energy does not necessarily imply that the excited 
nucleon is placed high in the continuum. 
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axes of the ellipsoid and low values of Q to nucleon motion in planes at 
right angles to this. This makes complete the qualitative correspondence 
with the collective model and suggests that, for example, experiments 


Fig. 2 
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The calculated absorption curve of a nucleus with 6=0 in the region of A =210. 
The computed transitions each smeared into a Gaussian of full width 
3 Mev, are grouped in equal energy intervals and are represented by the 
bold dots. The energy scale is in arbitrary units and we must set one 
of these units equal to about 2 Mev to fit the absolute experimental 
energy of the resonance correctly. 


ABSORPTION CURVE FOR A NUCLEUS OF S=0-1 AT A™195 


u > wu 


CROSS SECTION 
n 


GAMMA RAY ENERGY (ARBITRARY UNITS) 


The calculated absorption curve of a nucleus with §=0-1 in the region of A =195.. 
See also caption to fig. 2. 
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using oriented nuclei and, or, polarized gamma-rays would not distinguish 


between the models any more than simple measurements of the widths 
seems to do, 


Fig. 4 
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‘The calculated absorption curve of a nucleus with 5=0-2 in the region of 
A=185. See also caption to fig. 2. 


Fig. 5 
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"The calculated absorption curve of a nucleus with 5=0-3 in the region of 
A=165. See also caption to fig. 2. 
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Figure 6 shows the dependence of the resonance width on 6 for a hypo- 
thetical nucleus of A ~160 and variable deformation. The resonance 
maximum is taken as 15-5 Mev for this figure. 


Fig. 6 
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The calculated width of the giant resonance for a hypothetical nucleus of 
A~160 and variable deformation whose resonance maximum occurs at 
15-5 mev. The computed points are shown by the circles. 


Finally fig. 7 shows the comparison between the calculated resonance 
widths evaluated over the whole range of interest and the experimental 
results. In order to avoid systematic errors the experimental work from 
a single laboratory (the National Bureau of Standards, Washington) has 
been used in this comparisony. 


§ 6. Discussion 


Figure 7 shows that the independent particle model gives an excellent 
account of the great increase in giant resonance width in the region of 
strongly deformed nuclei. We must however stress that the computations, 


{+ The experimental data are quoted from the work of Okamoto and are 
private communications to him from the NBS group (Petree, Weiss and Fuller). 
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lengthy as they are, are not wholly satisfactory and that they probably 
somewhat over-estimate the resonance widths for the more deformed 
nuclei. We must also recognize that the Nilsson diagram is least reliable 
for the cases of large deformation and large energy difference with which 
we are chiefly concerned here. A full computational programme which. 


Fig. 7 
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The comparison between the computed and experimental widths of the giant. 
resonance through the region of the rare earths. Mono-isotopic elements 
or those with a strongly predominating isotope are shown by a small 
cross. For the others the horizontal bar spans the mass values of the 
abundant isotopes. All these experimental data come from the National 
Bureau of Standards, Washington. The vertical bars are not intended 
to indicate the accuracy of the experimental results. 


aims at a simultaneous account of the single-particle level spacings and 
the ground state properties in this region of the periodic table would be 
extremely complex and lengthy. Until it is carried through we may 
adopt the present results as a semi-quantitative indication that the inde- 
pendent particle model of the giant resonance has it well within its power 
to give an account of the experimental situation. 
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ABSTRACT 


The energy spectrum and angular distribution of the protons emitted from 
®°Ni on bombardment with 13-5 Mey neutrons have been studied using 
photographic emulsions. About 30% of the protons emitted can be attributed 
to direct interaction processes. 


§ 1. LyTRopUCTION 


IN a previous communication (March and Morton 1958) henceforth referred 
to as I, the authors reported their results on the protons emitted from *4Fe 
and *¢Fe; for each of these isotopes the angular distributions of the protons 
emitted with energies greater than 7 Mev showed a marked peak at forward 
angles. 

It was considered desirable to make a similar study of a different medium 
weight element. The isotope ®°Ni was chosen for study because its Q 
value for the n, p) reaction lies between those of **Fe and **Fe so the 
cross section should also be intermediate between those of the **Fe and 
56He, The effects of the Coulomb barrier in °°Ni should be similar to those 
in the isotopes of iron. 


§ 2. EXPERIMENTAL PROCEDURE 


The nickel isotope (99-2 % ®°Ni) consisted of a layer 7 mg/cm? in thickness 
and 1 em? in area electroplated on to a gold foil. ‘The exposure arrange- 
ments and the procedure for the examination of the photographic emulsions 
were those described in I. Only protons making angles between 0° and 
140° with the direction of the incident neutrons were measured. The 
total number of protons observed in the area of emulsion examined was 
2498 of which 792 were background protons. The neutron flux at the 
isotope layer was determined by scanning a known volume of emulsion 
for recoil protons as described in I. The neutron energy was 


13-5 +0:1 Mev. 


+ Communicated by the Authors. 


578 P. V. March and W. 'T. Morton on the Distributions of 


§ 3. RESULTS 
The angular distributions of the protons emitted from °Ni for three 


different: ranges of proton energy are given in fig. 1. The total cross 
section for the emission of protons was obtained by extrapolating the 
angular distribution to 180°. The value obtained was 233 + 25 mb. 
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Angular distributions of the protons emitted from Ni. 


The energy spectra of the protons emitted in three different angular 


intervals are given in fig. 2. 
The solid line curves in fig. 2 are the proton energy distributions expected 


for the (n, p) reaction on the basis of the statistical theory of Weisskopf 
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and Ewing (1940). The calculated energy distribution was normalized 
to the observed distribution for proton energies greater than 4 Mev in the 
angular range 90°-140°. The majority of the protons emitted with 
energies below 4 Mev are due to the (n, np) reaction and the broken line 
was drawn through the observed distribution of these protons for the 
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Energy spectra of the protons emitted from °°Ni. 


angular range 90°-140°. The curves drawn on the other two angular 
intervals were obtained by assuming that the protons from these two 
compounds nucleus processes were emitted isotropically. 
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§ 4, DiscussION 


The isotropic angular distribution of the protons emitted at backward 
angles and the good fit of the calculated energy spectrum to the distri- 
bution of protons with energy greater than 4 Mev emitted in the angular 
interval 90°-140° suggest that these protons result from the decay of a 
compound nucleus. 

The Q values for the reactions (n, p), (n, np) and (n, 2n) are —2-0, 
—9-6 and —11-5 Mev respectively, hence at a bombarding energy of 
13-5 Mev the maximum energy of a proton from the (n, np) reaction is 
3-9 Mev and at proton energies of about 2 Mev the (n, np) reaction is 
favoured relative to the (n, 2n) reaction. Thus the large number of low 
energy protons observed is probably due to the (n, np) reaction. Figures 1 
and 2 show that these protons are emitted isotropically in agreement with 
a compound nucleus process. 

The broken line on the angular distributions was drawn through the mean 
of the experimental points at backward angles and illustrates any departure 
from isotropy at forward angles. It is seen that the high energy protons. 
are emitted preferentially in the forward direction. ‘These protons are 
probably due to the direct interaction of the incident neutron with a 
proton of the target nucleus (e.g. Austern ef al. 1953, Brown and Muirhead 
1957). The energy spectra show that the energy spectrum of the protons 
from direct interaction is approximately flat above 4 Mev and there is no 
evidence of marked structure in their spectrum similar to that observed 
by Gugelot (1954) and Cohen (1957) in the inelastic scattering of protons 
by nuclei. 

If the isotropic part of the proton distributions is attributed to the 
decay of a compound nucleus and the excess at forward angles to a direct 
interaction, the cross sections obtained are, compound nucleus: 155 mb 
(87 mb (n, p) and 68 mb (n, np)) and direct interaction: 68 mb. 

The results for °°Ni bear a strong resemblance to those obtained for 
*He and **Fe described in I. Similar results for medium weight elements 
have been reported by Allen (1957) and by Brown et al. (1957) using 
photographic emulsions and by Colli and Facchini (1957) using a counter 
technique. Colli ef al. (private communication) have measured the 
angular distribution of the protons emitted from natural nickel (67-4 % ®8Ni, 
26-7% ®°Ni) and find that the proton spectrum is much more intense at 
forward than backward angles. 

An estimate of the nuclear temperature for the compound nucleus part 
of the (n,p) reaction may be obtained by plotting log, N (L/S, nlp 
against #,, where N(£,,) is the number of protons observed with an energy 
#, in the angular interval 90°-140° and o,,, is the cross section for the 
formation of the compound nucleus by a proton of energy £,. This plot 
is shown in fig. 3. The slope of the line gives the nuclear temperature at 
the average excitation energy. The change of slope at 4 mev illustrates 
the onset of the (n, np) reaction. The value obtained for the nuclear 
temperature for the (n, p) reaction is 1-0 Mev. 
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Fig. 3 


Lh 


Fs Sy GS RT Oe 
Proton Energy in MeV 


Proton spectra of ®°Ni corrected for Coulomb barrier penetrabilities. 
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ABSTRACT 


The age of minerals determined by radioactive methods agrees roughly with 
the age of the galaxy. This suggests that radioactive alpha-decay constants 
have probably changed by less than a factor of 3 or 4 during the last 3 or 
4x10 years. It is shown that this implies that many ‘Constants of 
Nature ’, particularly e, h and c probably change more slowly than 10~™ parts 
per year. The meson coupling constant also appears to change more slowly 
than this and the beta-coupling constants by less than 10-° parts per year. 


Tue question of a possible time-dependence of the ‘Constants of Nature’ 
is occasionally raised (e.g. Dirac 1938). The evidence for their constancy 
over great periods of time is not very good although Teller (1948) has given 
reasons for believing that the gravitational constant G probably changes 
by less than 10-® parts per year. The usual interpretation of the red-shift 
of the spectral lines from distant galaxies probably implies a comparable 
constancy of the relevant atomic constants. 

It seems that we can sharpen these estimates by considering the pheno- 
menon of alpha-particle decay. This is because the rate of this decay is 
governed chiefly by the penetration of the charged particles through the 
Coulomb barrier which surrounds the nucleus and so tends to depend 
exponentially on the constants. Small changes in the constants would 
bring large changes in the decay rates. 

The ‘age of the earth’ is estimated through the relative concentrations 
of the various radiogenic leads and their parent bodies in unaltered 
minerals (Nier e¢ al. 1941, Holmes 1946, 1947). The age so deduced, 
about 3-4 x 10° years, is just about the same as is estimated for the galaxy 
(see e.g. Bondi 1952), the best current estimates running at about 7-8 x 10° 
years. We may take it that this coincidence is not due to chance. The 
implication is then that the decay constants A of the various radioactive 
bodies concerned in these estimates (?°°U, ?88U, 282Th) have not changed 
by more than a factor of 3 or 4 over the last 3 or 4x 109 years. Supple- 
mentary evidence of an essentially similar nature comes from a study of 
pleochroic haloes (Henderson 1934). 


The decay constant A is dominated by the penetrability factor so that 
we may write, to a first approximation ; 


A= constant x exp(— A) 


{+ Communicated by the Author. 
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where A=87*Ze?/hv. Here Z is the atomic number of the radioactive 
body and v is the speed with which the alpha-particle of energy H is 
emitted. This approximation, which neglects the effect of the nuclear 
radius on the decay rate, would be very poor for calculating the actual 
barrier penetrability but is quite adequate for our present crude estimate 
of the way in which A depends on e. 

If now e were to change slightly there would be no first order effect. on 
the nuclear wave functions since these are determined chiefly by the much 
stronger nuclear forces which we leave unchanged for the present argument. 
We may for example estimate that a 1°% increase in e would dilate the 
nucleus by less than 0-1°%. The intrinsic probability of the alpha-particle 
emission therefore changes very little and any effect on \ is due chiefly to 
changes in A through e and v. The change in v comes about because the 
Coulomb energy changes with e. Writing the Coulomb energy classically 
and ignoring the effect of the change in R, which is unimportant, we have: 


dv fde_ 12€2Z 


v/ e 5 RE’ 


With the earlier equation this gives: 
2 2 
s/s = -A(2- 25 z) =1-5 x 103 


for the bodies of interest. We have used R=1-2A" fermis. 

This estimate is in error because we have neglected the effect of the 
radius in diminishing A and have neglected the effect of the change in 
radius in reducing the change in v. When these omissions are corrected 
the above simple estimate of the sensitivity of A to e is lowered by about 
20%. 

We see that the increase in decay rate due to the increase in v on in- 
creasing e more than offsets the extra difficulty in escape due to the 
increased height of the Coulomb barrier. This we can understand 
qualitatively. As we increase e to the first approximation the height of 
the Coulomb barrier at the nuclear surface increases but the depth of the 
nuclear potential below the top of the barrier remains the same so the 
whole nuclear level system is simply lifted with the barrier, the alpha- 
particle escaping with greater energy but at the same energy below the 
top of the barrier. The barrier however has become thinner at a given 
energy below its top because it has been raised so the penetrability 
increases. When we take into account also the increased Coulomb effects 
within the nucleus itself and not just up to the nuclear surface we see that 
in fact the energy at which the alpha-particle escapes below the top of the 
barrier decreases as the barrier is raised so on this account also the pene- 
trability increases. 

We noted earlier that a change of A by a factor of 3 or 4 in 3 or 4x 10° 
years is unlikely so the present result implies that a rate of change of e by 
10-12 parts per year is unlikely. 
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A similar analysis of the sensitivity of A to changes in h (through A) 
suggests that h changes by less than 107" parts per year but here the 
situation is complicated because of the involvement of / in the nuclear 
wave functions. However it is unlikely that the effect of a change in h, 
through the nuclear wave function, if of the opposite sign, would nearly 
cancel that due to the change in penetrability. We may therefore let this 
estimate stand as probably conservative. 

We may likewise discuss the stability of the nuclear forces. It is difficult 
to do this from a fundamental viewpoint because we are not yet able to 
relate the real part of the effective nuclear potential to the internucleonic 
forces in any simple way. We may however discuss possible changes in 
the effective nuclear potential itself. If we consider the alpha-particle to 
be formed from nucleons which are near the top of a Fermi-sea in a potential 
of depth V and of fixed radius then : 

da [dV _ pom 2e, * si =5 x 10*, 

ah ae hv E 
We have set V=55Mev. The sensitivity of A to changes in F for a potential 
of fixed depth is roughly twice as great as this. It therefore seems likely 
that the parameters of the effective nuclear potential change by less than 
10-1? parts per year. Since the nuclear forces certainly depend on h/me 
where m is a mesonic mass it seems likely that h, c and the masses of 
elementary particles themselves change by less than this order per year. 
We may similarly include the (7) meson—nucleon coupling constants in 
this list. 

We can also make some remarks about the constancy of the coupling 
constants for beta-decay. From the pleochroic haloes in mica we get a 
picture of the natural radioactive series integrated over geological time 
(Henderson 1934, Henderson and Turnbull 1934). These series contain 
certain bodies which at the present epoch emit alpha-particles and beta- 
rays in competition with each other. One such point in the 2°5U series is 
at Ac which branches 1-2° by alpha emission to AcK, the rest by beta 
emission to RdAc. Another point is at AcC which branches 99-7% by 
alpha emission to AcC” and the rest by beta emission to AcC’. A parti- 
cularly well-measured specimen is the Renfrew biotite of Henderson and 
Turnbull (1934) which is about 1x 10° years old. It shows good haloes 
called F and H of radii 15-4 and 19-4 due respectively to Ra+ UII + Io 
and to Rn+Rak. The halo due to alpha-particles from Ac is expected 
at 17-5 4 radius which is almost exactly in the middle of a rather clean and 
halo-free gap between haloes Fand H. This is consistent with the present 
day branching of a factor of 80 in favour of beta-decay of Ac. If however 
10° years ago the beta-coupling constant had been 10 times less the rate of 
decay of Ac by alpha-emission would have been equal to the rate of decay 
by beta-emission and a visible halo would have formed in the present gap. 
Similarly had the constant been greater by this same order a halo due to 
the alpha-particles of AcC’ should be found at a radius of 32:6 u in another 
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rather clean gap between halo C at 27:41 (due to AcC + An) and halo A at 
34-6 (due to RaC’). This suggests that the beta-coupling constants do 
not change in either direction by 10-8 parts per year. 

We must stress that all these estimates are crude but probably good to 
orders of magnitude and, where open to considerable error, as in the 
discussion of nuclear forces, are probably conservative. We must also 
admit that we cannot detect a conspiracy among the Constants of Nature 
by which some or all might change in such a way as to preserve the con- 
staney of A, This, however, would surely be the ultimate luwsus natwrae. 
Nor could we detect rapid fluctuations of the constants (say with a period 
of 107 years) about their mean values. We have effectively assumed that 
their change is monotonic in time. 


Note added in proof.—t am grateful to Dr. K. I. Mayne for pointing out 
to me that the ages of minerals determined by the A—K and Sr-Rb methods, 
which depend on the beta-couplings, are in agreement with the ages 
determined by alpha-decay methods to a degree which implies constancy 
of the beta-coupling constants to the order of 10-° parts per year. See, 
e.g. J. T. Wilson, R. D. Russell and R. M. Farquhar, 1956, Handbuch 
der Physik, 47, 288, for modern references. 
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ABSTRACT 


The theory of Landau and Khalatnikov for the viscosity of liquid He IT is 
examined in the neighbourhood of the lambda point. Arguments are 
advanced against the validity of the ideal gas approximation for rotons, 
and a qualitative roton liquid model is proposed, The temperature depend- 
ence of viscosity of the liquid model is found to correspond to the observed 
viscosity above 1-8°x, and with the discontinuity in the temperature deriva- 
tive of viscosity at the lambda point. 


THE viscosity of liquid He II, as measured by the technique of oscillating 
discs or by rotation viscometer, rises monotonically as the temperature 
is decreased from 1:8° to 1°K. This trend is successfully described by 
Landau and Khalatnikov (1949, hereafter referred to as LK) in terms of 
the phonon and roton excitations of liquid “He below its lambda point 
temperature 7,. At about 1-8°K, however, the observed temperature 
derivative of viscosity changes sign, the viscosity increasing by about a 
factor of two between 1-8°K and 7’,, and a theoretical explanation of this 
behaviour has not yet been proposed. The success of the LK model at 
lower temperatures can be taken as evidence for the validity of their 
theory ; we here suggest a means for extending it to temperatures near 7’,. 

The detailed calculations of LK show that the important processes 
giving rise to viscosity above 1°K are the scattering of phonons by rotons 
and of rotons by rotons. The two types of processes are each characterized 
by separate coefficients, respectively called ‘phoncn viscosity’ 7, and 
‘roton viscosity ’ 7,, and the total viscosity » of the liquid is given by the 
sum of these two contributions. 

The temperature dependence of »,, is given by LK as 


» = const./N, 


where NV, is the number of rotons per unit volume of liquid. Since the 
density p, of normal fluid at temperatures 7’>1°K is almost wholly due 
to the rotons, and N,, is proportional to p,7', the phonon viscosity can be 
written in terms of the measured fraction 7=p,/p of normal fluid, 


Np) =a/2l, a=constant, ~. |. . 4. . 4 


+ Communicated by the Author. Guggenheim Fellow 1957-8, on leave 
from Los Alamos Scientific Laboratory. 
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for the region 7>1°K. A re-examination by Khalatnikov (1956) of the 
characteristic times for achieving equilibrium among the scattered phonons 
indicates that the coefficient a of eqn. (1) varies with the temperature, 
but much more weakly than x. The coefficient is known in terms of certain 
empirical parameters of the liquid, and in that sense, the value of a is not 
completely arbitrary. There are, however, large discrepancies among the 
several experimental determinations of the fundamental parameters, 
leading to a considerable uncertainty in the magnitude of the phonon 
viscosity. We shall therefore treat the coefficient a itself as an empirical 
parameter. In addition, we neglect its variation with temperature. 
Errors that may result from this simplification may lead to a poorer 
description of the data near 1°K, where the phonon viscosity is dominant. 
Such errors will not affect the principal object of this paper, which is 
concerned with the behaviour of roton viscosity at higher temperatures. 

The roton viscosity is found by LK to be temperature independent. 
This result stems from the high self-energy A compared to the kinetic 
energy ~k7T’, the assumption of a hard-sphere roton—roton interaction 
potential, and the assumption of relatively low density N,. Under these 
conditions, the rotons behave as quasi-particles, having a viscosity inde- 
pendent of W,, as characteristic of an ideal gas of real massive molecules. 
This prediction appears to be confirmed by experiment, for the value 
of » (exp.) —7, (theor.) is about 10 micropoise between 1-2° and 1:8°K 
although « changes by a factor of 15 in this interval. We should, however, 
accept the result with some reservation, since 7, becomes increasingly 
dominant at low temperatures and the many measurements of viscosity 
do not agree well below 1-6°K. 

The high temperature region of He II, from 1-8°K to 7',, is the region 
in which the behaviour of 7, is dominant, and it is in this range that the 
LK theory does not correspond to the experiments. The failure, as 
Landau and Khalatnikov indicate, arises from the use of the ideal gas model. 
Such an approximation is justified in the case of low density, when each 
‘ particle ’ spends only a small fraction of its time under the influence of 
others. In the low density limit, the roton—roton interactions can be 
described as elastic two-body collisions, and the empirical collision cross 
section contains in integral form whatever moderate range interaction 
potential may exist. As in the case of real molecules, details of the roton— 
roton potential become increasingly important as the density increases. 
The existence of such a potential is indicated by the theoretical examination 
of phonon exchange processes given by Cohen and Feynman (1957), and 
by Khalatnikov’s (1956) semi-empirical description of the ‘ non-ideality ’ 
of the roton gas. Furthermore, the LK ideal gas approximation suggests 
a high probability for many-body collisions, even at 1:8°K; on the basis 
of their calculations, the roton—roton collision-radius r (related to the 
collision cross section « by o=7r?) at this temperature is of the order of 
the mean free path, and the time between collisions is ~ 107 sec, com- 
parable with the period of atomic vibrations. A gas model seems 
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completely unsuitable in these circumstances, and it should be more 
profitable to compare the assembly of rotons to the molecules of a liquid. 
Such a comparison is not too badly handicapped by our ignorance of the 
detailed interaction as long as we are primarily interested in the qualitative 
dependence of viscosity on temperature and density, and can forgo the 
much more difficult problem of calculation of magnitudes. For our 
purposes, the qualitative theory of viscosity of liquids (Frenkel 1946) 
is adequate. We therefore accept, for the description of the motion of 
rotons in the high density region, the microscopic model of liquids, called 
at times the ‘ hole theory ’, which has achieved moderate success in des- 
cribing the dependence of viscosity of ordinary simple liquids upon tem- 
perature and pressure. According to this model, the motion of a molecule 
in the liquid is considered to be quite constrained by the close presence of 
its immediate neighbours; in contrast to the trajectories in a gas of low 
density, it is here taken to correspond to a sporadic migration through an 
equilibrium configuration of other molecules. Details of the interaction 
need not be known if we assume a density-dependent activation energy W 
opposing the dislocation of an individual molecule from its ‘ equilibrium ’ 
position. The viscosity then behaves according to the exponential 
law 
Thig= A exp (V7), 


where A is approximately constant. The variation of 7 with density 
is given by expanding W about its value W, corresponding to a standard 
state of numerical density 7, ; 


W=W,+a(n—M)+PB(n—n)?+.... 


Using the proportionality of V, and normal fluid fraction x above 1°K, 
and neglecting terms higher than the first, we obtain for the viscosity of 
the ‘ roton liquid ’ ; 
"=O exp (Cry 2) ee eee 
where 6 and ¢ are approximately constant : that is, neglecting the variation 
of T’ in comparison with the much more strongly varying function (7). 
Equation (2) corresponds approximately to ‘ Andrade’s Law’ (Andrade 
1934) relating the viscosity and density of ordinary liquids. The equation 
for 7, is also of the correct form in the region of low roton density, for it 
approaches the constant value 6 at low temperatures as for a gas; this isa 
fortuitous result of the rapid variation of x with temperature. We there- 
fore obtain, for the viscosity of normal fluid in He II in the range 1°K-T",, 
the sum of roton and phonon contributions . ( 


ee b a 
Uracil exp (cx). ERS et ee ee 


Comparison of eqn. (3) with the observed viscosity (Dash and Taylor 1957) 
in the range 1-2°K~7', yields the empirical constants 


a=0°675 uP deg; b=10 uP; c=0-95, 


——_— = 
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which may be compared to the values 
a~0-8 uP deg; b=11-5yP; c=0, 


given by Khalatnikov (1956). Equation (3) appears to correspond 
closely to the observed viscosity above 1-2°K, as shown in the figure; 


mean deviation between our equation and the measured values is of the 
order of 5%. 


Viscosity 7) in micropoises 


Tek 


Viscosity of He II above 1-2°K. 
—---- Np, Landau and Khalatnikov. 
— + — =: »,, liquid model. 


> n +7,. lod 
+ , smoothed experimental values, Dash and Taylor (1957), 


This liquid model of roton viscosity appears to present a simple explana- 
tion for the observed singularity in dy/dT at the lambda point. Experi- 
ments (Dash and Taylor 1957, Taylor and Dash, 1957, Welber, private 
communication) have shown that both dy/d7' and dp,/dT’ have strong 
discontinuities at 7',, but that the kinematic viscosity 7/p,, changes quite 
smoothly through the transition temperature. From the latter observa- 


tion we have 
dT’ Pn T,—aT d7 Pn T,+aT 


dn = (3) ~ 0, 
an T,+aT dT T,+aT 


dy - (#4) - (2) (FF) aah (4) 
dT [,—dT aT’ T,+aT P/T) dT T,—dT 


and since 
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From eqn. (2) we obtain the temperature derivatives of the roton liquid 
viscosity on each side of the lambda point : 


(Br (Basar AP) eBoy ar 
aT’ T,+aT  NaT T,—aT P/T) aT’ T,—aT 


Since the empirical value c= 0-95, this is seen very nearly to satisfy the 
observed condition, eqn. (4). Thus, although we do not consider the 
mechanisms of viscosity in He I, the discontinuity in dy/dT' at the lambda 
point appears to arise as a direct consequence of the abrupt halt in the 
increase of roton density with temperature. 
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ABSTRACT 


The critical magnetic field of thorium has been determined between 
0-1°K and 1-37°xK, the transition temperature, where the initial slope of the 
critical field curve is 190 gauss/deg. The data indicate a value of 162 gauss 
for the critical field at absolute zero. The critical field curve departs from 
a parabola near the transition temperature. The electronic specific heat, 
11-1x10~*7 cal/mole-deg, deduced from the magnetic measurements is in 
substantial agreement with the previously determined calorimetric value. 


§ 1. INTRODUCTION 


As the result of electrical resistivity measurements, Meissner (1929) 
reported the superconductivity of thorium at about 1:4°K. Later 
Shoenberg (1939) investigated the critical field of thorium down to 1°K 
using a moment measuring technique. His measurements indicated that 
thorium had superconducting properties similar to those found in ‘ soft ’ 
superconductors; in particular, the superconducting-normal transitions 
in a magnetic field were very sharp. Since extremely pure thorium is now 
available, we thought it would be of interest to extend the critical field 
measurements below 1°K to determine the shape of the critical field curve 
in this region and also to obtain a magnetic value for the electronic 
specific heat coefficient, y. 


§ 2. EXPERIMENTAL DETAILS 


A polycrystalline rod of thorium prepared by the van Arkel process was 
obtained from Metal Hydrides, Inc., Beverly, Massachusetts. The 
chemical analysis provided with the material indicated the following 
impurities in parts per million: 

NG 3018 Cumi2-6 0Al* 2°50 28 

ergon Nie 4-08 Ca 2°5° Sn, Co, Zn,Zr, Ge, Ps As; Cd,iSb, 

Si 33 Pb 2-5 Mg 2-5 Ba, Be, V, Mo, Bi, Ta, Hg not detected 

Grace Mie 2-5 6030-0 


The specimen was in the form of a rod 6 mm in diameter and 25 mm 
long. It was cemented (GE 7031 adhesive) into two split copper cylinders; 
one of these was in thermal contact with a pill of potassium chrome alum, 
while a carbon resistor used as a resistance thermometer was cemented to 
the other. In order to reduce heat leak from the surrounding liquid 


+ Communicated by the Authors. 


592 N. M. Wolcott and R. A. Hein on the 


helium bath to the sample assembly, the suspension threads and leads 
to the carbon resistor were passed through a ballast pill of manganous 
ammonium sulphate. 

Temperatures below 1°K were obtained by adiabatic demagnetization 
of the chrome alum from fields up to 30 000 gauss at an initial temperature 
of 1:3°x. After demagnetization, the procedure was to observe the 
temperature (differential susceptibility) of the chrome alum pill by a 
ballistic mutual inductance technique and to follow the resistance of the 
carbon thermometer. Comparison of the two temperatures indicated 
that the time for thermal equilibrium between the salt and the thorium 
was about two minutes. A steady external field was then applied and 
the differential susceptibility of the thorium was observed as the system 
slowly warmed up. As the chrome alum and the thorium were centred 
in different measuring coils, the differential susceptibility of each could 
be measured independently. When the transition temperature was 
reached for the particular applied field, a sudden change in the differential 
susceptibility of the thorium was observed. The applied field was then 
reduced until the metal again became superconducting. In this way, 
two to three points on the critical field curve could be obtained from a 
single demagnetization. 


§ 3. RESULTS 


A typical warm up curve is shown in fig. 1. The solid curve is the 
galvanometer deflection (proportional to the differential susceptibility 
of the thorium) obtained with the ballistic measuring field parallel to the 
applied external field, while the dotted curve represents the galvano- 
meter deflection obtained with measuring field opposite to the applied 
field. The deflections obtained with the metal in the normal and super- 
conducting states respectively are shown on the margin and denoted by 
NandSC. The magnitude of the change observed between the completely 
superconducting and completely normal deflections indicates that 95 to 
100% of the thorium became superconducting. 

The qualitative interpretation of the curve is as follows: Following 
demagnetization a galvanometer deflection corresponding to complete 
superconductivity is observed. An external field is then applied until 
the galvanometer deflection indicates a differential paramagnetism 
(¢=1 hour in fig. 1). The excess differential paramagnetism above the 
normal is due to the inverse Meissner effect, for near the transition 
temperature in a magnetic field the small applied measuring field (about 
one gauss) causes a large amount of expelled flux to re-enter the super- 
conductor, thus giving rise to an apparent paramagnetic differential 
moment which can be quite large. For about one hour after the measure- 
ments began, the applied external field dropped as the solenoid coils 
warmed up. This caused the galvanometer deflection to decrease during 
this period. The deflection then remained very near the normal value 
for about two hours. In fact during the first experiment, it was suspected 


EEE 
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that the metal might actually be normal; and only the fact that identical 
deflections were not obtained with the measuring field reversed indicated 
that superconducting material was still present. As the transition 
temperature is approached, the inverse Meissner effect becomes more 
pronounced and large paramagnetic deflections are obtained. This 
‘ differential > paramagnetic effect has been discussed by Steele (1952). 
In this case the maximum paramagnetic deflection was more than ten 
times the diamagnetic deflection corresponding to complete super- 
conductivity. As the transition temperature is passed, the deflection 
suddenly returns to its normal value. When the measuring field is 
reversed (dotted curve in fig. 1), the paramagnetic deflections are shifted 
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Galvanometer deflection as a function of temperature and time after 
demagnetization. 


) Measuring field parallel to applied field. 
---A--- Measuring field anti-parallel to applied field. 

N Galvanometer deflection in normal state. 

SC Galvanometer deflection in superconducting state. 


to higher temperatures since the external field is reduced by the measuring 
field; and the metal becomes normal at a slightly higher temperature. 
The transition temperature in the present experiments was defined as the 
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temperature at which the differential susceptibility returned to its normal 
value. The sharpness of the differential paramagnetic effect allowed this 
temperature to be obtained very accurately; this is in contrast to the 
results obtained with uranium, another actinide element, for which 
Hein e¢ al. (1957) found spread-out transitions with no sharp Meissner 
effect. 
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Critical magnetic field curve of thorium versus temperature. 


With no applied field, the transition temperature was 1-37°K(7',) in 
agreement with the earlier value of Shoenberg. The transition tem- 
peratures in various applied fields are shown in fig. 2 along with the 
earlier results of Shoenberg. It appears that the critical field departs 
from a parabola (shown as a solid line) near 7',. The measured points 
indicate that the slope of the critical field curve at the transition tem- 
perature, (dH,/dT')r—7,, is 190 gauss/deg, also in agreement with Shoenberg. 


§ 4. Discussion 


In order to obtain an estimate of H), the critical field at absolute zero 
H, was plotted vs. T? as shown in fig. 3. An extrapolation of this curve 
(which is linear over most of its range) to 7’=0 yields a value of H, of 
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162 gauss. Over the region where H/, is linear in 7, the critical field is 
given by H,=H,—AT® with H,=162 and \}=90-7; this is the parabola 
which is indicated by the solid line in fig. 2. 

From the critical field curve a value of the electronic specific heat 
coefficient y can be obtained, provided that the critical field curve 
represents the boundary for reversible superconducting—normal tran- 
sitions. Application of thermodynamics to the transition indicates that 
(Daunt 1955) in general y= —(V/87)(d?H,/dT?),_ where V is the molar 
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Critical magnetic field versus 7? for thorium. 


volume. If the critical field curve is parabolic in 7, then this formula 
can be rewritten y=(V/87)(dH,/dT)?2=7,. When information about the 
critical field curve is only known near 7’,, then this formula must be used. 
Using Shoenberg’s value for (dH,/dT)r-r,, Daunt obtained a value of 
y = 7-1 x 10-4 cal/mole-deg?. Since the present measurements show that 
the critical field curve is not parabolic, it is not surprising that this value 
of y disagrees with the calorimetric one. 

In order to determine y from a measured critical field curve near 7’=0, 
it is necessary to evaluate the general relationship for y. For this purpose. 


232 
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it is convenient to change the variable of differentiation from 7’ to T?; 
when this is done, it becomes 


a AVLEL 
Y a 2Q7 aT? T2—=(0 : 


Both H, and (dH,/d7T?),,_,) can be obtained from fig. 3, the former being 
the intercept and the latter the slope of the line at 7?=0. Using the 
values H)=162 and (dH,/dT?)m_)= — 90-7, the value y=11-1 x 10 cal/ 
mole-deg? is obtained. This is in good agreement with the calorimetric 
value of y= 11-2 x 10-4 cal/mole-deg? reported by Smith and Wolcott (1955). 

In summary, we may say that thorium behaves in many ways like an 
ideal superconductor, that the critical field curve though not a parabola 
is parabolic over much of its range, and that values for the electronic 
specific heat coefficient, y, obtained from magnetic and calorimetric data 
are in substantial agreement. 
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ABSTRACT 


A model of cleavage applicable to metals cleaving on the slip plane is 
developed in which a crack is initiated at the end of a low angle tilt boundary 
terminating inside the crystal. The strength of the metal is determined by 
the difficulty of growth of the crack. Satisfactory agreement is obtained 
with the experimental results for zinc. 


§ 1. IyrTRopUCTION 


RECENTLY attempts have been made to account for the brittle fracture of 
metals in terms of the stress concentrations due to piled-up groups of 
dislocations at the ends of slip lines. Though this treatment has had some 
success with polycrystals (cf., for example, Stroh 1957) two difficulties are 
met when we try to apply it to single crystals. First, it is not easy to 
think of suitable obstacles in the slip plane able to withstand the large 
piled-up groups needed for fracture, since the grain boundaries which 
fulfill this role in polycrystals are no longer available. Secondly, a number 
of metals, of which zine and bismuth are examples, have a single pre- 
dominant slip plane which is also the cleavage plane; and since piled-up 
groups produce no tensile stress normal to their own slip plane they would 
not produce in these crystals the kind of stress magnification needed for 
fracture even if they could be formed. It seems therefore worthwhile to 
consider other possible mechanisms for fracture; in this paper a mech- 
anism which is particularly adapted for the case in which slip and cleavage 
planes coincide is investigated. 


§ 2. InrTraTION OF CRACKS 


Orowan (1954) has pointed out that if a sub-grain boundary, i.e. a wall 
of dislocations, terminates within a crystal there will be a high stress con- 
centration at its end which may be great enough to produce a crack. Such 
a terminating wall may be formed from a complete wall extending right 
across the crystal if part of this is held up against an obstacle and the 
other part is pulled on by the applied stress (fig. 1). Friedel (1956) has 
estimated the stresses which will occur at the end of such a wall and con- 
cludes that a crack should form if the wall corresponds to a misorientation 
of 5° or more; this estimate is necessarily somewhat rough as it reflects 
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our ignorance of the exact value of the theoretical strength of the ideal 
lattice. This uncertainty, however, is of no great importance for the 
application to be made here, for our treatment will be based on the idea 
that the strength of the crystal which is observed is determined by the 
growth of the crack and not by its initiation ; then we need only know that, 
for some sufficiently great misorientation, a crack can be initiated at the 
end of a wall of dislocations. 

Cracks formed in this way have been directly observed by Gilman (1954) ; 
fig. 2, Pl. 25 shows two such cracks formed in zine under a compressive load. 
A second example of crack formation, also due to Gilman (1957), is shown 
in fig. 3, Pl. 25 and we suggest that the same mechanism is operating here ; 
we consider this case in some detail. The figure shows a symmetric 


Fig. | 


a | 
4 


Cleavage at a divided dislocation wall. 


bicrystal of zinc which has been tested under tension, the tensile axis 
being parallel to the grain boundary; since such bicrystals were found to 
have the same yield stress, rate of strain hardening and ductility as single 
crystals it is reasonable to assume that the same processes occur in these 
bicrystals and in single crystals. If we compare slip lines such as A and 
B lying on either side of the crack we see that while the slip lines on the 
side A are very nearly straight, those on side B are curved: thus the slip 
lines B must contain an excess of dislocations of one sign, but not the slip 
lines A. Since the slip lines are curved and not polygonized, the dis- 
locations cannot all lie in a single boundary plane as in fig. | and fig. 2, Pl. 25; 
nevertheless the curvature occurs only in a fairly narrow region and so it will 
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be a reasonable approximation to suppose that here also the dislocations 
form a plane sub-boundary. The direction of bending of the slip planes 
shows that the extra half atomic planes terminating on the dislocations 
are on the side further from the crack; under the applied stress (a tension 
parallel to the grain boundary) such dislocations would move from left 
to right in the figure. 

We therefore suggest that the following sequence of events has occurred. 
A dislocation wall formed on the left of the grain boundary and under 
the applied stress moved towards the grain boundary. Though the 
orientations of the two crystals are such that the boundary is transparent 
to dislocations and slip can continue through it from one crystal to the 
other, the grain boundary will nevertheless offer some resistance to the 
passage of dislocations. In the present case part of the wall has been held 
up at the grain boundary while the remaining part has been pulled through 
it into the second grain. Thus the wall has been split into two, and we 
have the situation shown in fig. 1 where a crack has been formed by the 
stress concentration at the end of the wall. 

The misorientations in figs. 2 and 3, Pl. 25 are 14° and 8° respectively ; 
Friedel’s estimate of 5° as minimum misorientation to produce a crack is 
consistent with these results. 

Some sort of barrier, for example the grain boundary in the case shown 
in fig. 3, Pl. 25, is still needed to hold back part of the wall while the other 
part is pulled on by the applied stress. But now the barrier may be relat- 
ively weak since a large number of dislocations can become blocked simul- 
taneously, whereas in a piled-up group only the leading dislocation is 
blocked and the other dislocations all push against it. A possible barrier 
is now a low angle boundary formed from dislocations intersecting the 
slip plane; such a boundary corresponding to a misorientation of 1° or 
even less should be quite adequate as a barrier. 

It is interesting to note that the present mechanism for initiating a 
crack accounts for the tendency for fracture (in single crystals of zinc, for 
example) to occur near the grips in a tensile test if the greatest possible 
care is not taken in mounting the specimen. For then inhomogeneous 
strains will develop near the grips and low angle boundaries will form 
most readily here. 


§ 3. SPREAD OF THE CRACK: APPROXIMATE TREATMENT 


We shall now accept that cracks can be initiated in the manner dis- 
cussed above and consider the conditions under which they will grow; 
the stress at which this occurs will determine the strength of the crystal. 

When a wall of dislocations divides into two, the dislocations on either 
side of the division will attract each other and so oppose the separation of 
the two parts of the wall. If the distance between these parts is x the 
number of dislocations in each which interact is of order a/h, where h is 
the spacing of the dislocations in the wall. There are thus (a/h)? pairs of 
dislocations each interacting with a force of about (b?/27x, and so the 
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resulting force tending to reunite the wall is Gb2x/ 2h? = GO?x/27, with 
§=b/h the misorientation across the wall. If the applied load produces a 
resolved shear stress o, and so a force o,b on each dislocation, the total 
force on a wall of length L is o,b(L/h)=o,L@. Hence the two parts of the 
wall will attain an equilibrium separation « given by 


OillnaoD. oe er 


Now consideration of the geometry of fig. 1 shows that, as far as its 
effect in causing the crack to spread is concerned, the divided wall is 
equivalent to a wedge of thickness Ox driven into the cleavage plane ; such 
a wedge may be shown (Stroh 1954) to produce a crack of length 


= (OxG [Say 5 2 ee ee 

where y is the energy per unit area of the new surfaces formed. On 
substituting for 0x from (1), this becomes 

Cera ALA 2yG) 0 ae: Cane ee 

Once it has formed the crack will be subjected to the normal component 

oc, of the applied stress which will extend it in the manner considered by 


Griffith (1920); the condition that a crack of length ¢ should spread 
catastrophically is then 


077°C = 8yG'/ 7. 40 a 9.31. 2 a Pe 
From (3) and (4) we obtain the condition for fracture of the material as 
6 O,l = 4yG4/r. on ae sat Gal a 


This equation which constitutes our main result will be discussed in § 5. 
The derivation which has just been given, however, involves a number of 
simplifying approximations which have been made so that the physical 
ideas involved will not be obscured by purely mathematical complexities. 
We therefore reconsider the argument in the following section under less 
restrictive assumptions. 


§ 4. SPREAD OF CRACK: Morr Exact TREATMENT 


Our previous treatment was‘deficient in two respects which will be 
remedied here. First, various approximations were introduced which 
result in the constant on the right-hand side of eqn. (5) being incorrect by 
a numerical factor. Secondly, the material was assumed to be elasti- 
cally isotropic which will certainly not be the case for single crystals. As 
we shall be particularly interested in the variation of the strength with 
crystal orientation, we wish to be sure that no errors are introduced by 
neglecting the anisotropy of the crystal. Though in this paper we shall 
confine our discussion of experimental results to those for zinc, we shall 
not assume hexagonal symmetry, but consider the general triclinic case. 
In this way our equations will be directly applicable to other crystals when 
experimental results for these become available. This will not complicate 
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the treatment in any way. The main approximation which remains is 
that we assume that, after the dislocation wall has divided, its two parts 
are plane; in fact they are likely to curve slightly under their mutual 
attraction. 

A general discussion of some problems in anisotropic elasticity has been 
given by Stroh (1958) and we follow here the treatment and notation of that 
paper, referring to it for fuller details. We take wv, axis parallel to the 
dislocation lines so that the stresses depend only on 2, and 2, and intro- 
duce three complex variables z,=2, + p,2, («=1, 2, 3) and their complex 
conjugates 2,=2,+p,%.; p, and p, are defined as the roots of a certain 
sextic equation whose coefficients are the elastic constants. We shall 
also need a set of complex constants L;, (¢, «=1, 2,3) which depend on 
the elastic constants, and we denote the reciprocal of the matrix [L;,] by 
[M,;|. d; are the components of a real vector which is essentially the 
product of an elastic modulus and the Burgers vector. All these constants 
will be eliminated from our final equations and we do not need here to know 
their precise values. It may be shown that three of the components of 
the stress field due to a dislocation lying along the x, axis are 

Galan (L200 2M.) 0, . 2 - ©(6) 


ia~a o) 
a 


where we have used the convention that a repeated Latin suffix is to be 
summed. (Summation over Greek suffices will always be indicated 
explicitly as in eqn. (6).) The stress components not included in (6) will 
not be needed here. 

We have now to sum the stresses (6) over all the dislocations in the 
divided wall, the position of which we suppose is given by 7,=0, %,<0 
and 2,=x,2%,>0. In finding the stress at a point P(2,, x.) at a distance 
from the wall large compared with the spacing 4 between the dislocations 
we may replace the summation by an integration. Then considering the 
part of the wall x,=0, x,<0, we see that the factor z,~1 in (6) must be 
replaced by 

[es eer 
racy na RA ERE Ta 
if the length L’ of this section of the wall is large. In the same way we 
obtain from the section of the wall of length Z in the plane x,=¥«, a con- 


tribution of 


1 1 


7, 08 (2, —% + py%o) — p 5, OB Pals 


and taking both sections of the wall together we see that z,~' in (6) must be 
replaced by . 
(ph) log (2, —x)/z, — (ph) log (L’/L) + 7t/p,h . . (7) 


to give the stress of the whole wall. Now we specify that the wall is to be 
a simple tilt boundary. Such a wall, when undivided, produces no remote 
stress field; hence the stresses we obtain must vanish when «=0, and so 
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the terms independent of x in (7) can give no contribution to the stress. 
Dropping these terms we obtain from (6) and (7) the stresses due to the 
divided wall 
Ci = (1/47h) 2 NE a log (Z, i 2) [ey pe LM jf. log (Z, = a) [Z\ d;. 
& ee LS) 
The force on each dislocation of the wall due to the remainder is boy», 
and the total force on one section of the wall is obtained by summing the 
forces on all the dislocations composing it. Replacing as before the 
summation by an integration, we find the attractive force between the 
two parts of the wall is 


F = (ba/2ah?){u,(1 + log (L/x)) +p}, 9. - » «ss = (9) 
where : = 

beet > (Ley Mt Lea ew) d,/b, ‘ " a ¢ 2 (10) 
and ; a 

bo= 3 > (UaaDe log p,M 4; ‘Phas Dy log p,M ,;) 4,6 ; * (11) 


in obtaining this, we have used the result that the constants are related by 
Ly,+Pelis,=9. Since d,; is proportional to 6, ~; and «, depend only on 
the elastic constants of the crystal. The force due to the applied stress 
separating the two parts of the wall is o,bL/h and on equating this to F, 
eqn. (9) shows that the equilibrium separation x is given by 


a{u,(1 + log L/a) + pg} = 2arhLo,/b. 7: Ares ire See 


Now suppose that a crack has been formed in the plane «,=0 by the 
stresses of the dislocations as in fig. 1. The normal stress o,, across this 
plane is from eqn. (8) and (10) 


Oo9 = (bpt4/2h) log |(x, — x) /2r,|. soto) vo Rene bo 
We choose as unit of length half the length of the crack, so that the crack 
extends from 7,=0 to x,=2. Stroh (1958) has shown that the stresses 
near the tip of the crack x,=2 under a variable stress o,, are the same as 
if a uniform tension 

T,=n7! | Cont U2(2—a,)H2da, . . . . . (14) 


0 
acted on the crack. Substituting (13) in (14) we find that in the present 
case the mean tension is 
or if the crack has length 2c the effective tension on it is 
T=, oi, 7/27he, ee (15) 


where the normal component o,, of the applied stress has now been added. 

We assume that the spread of the crack depends only on the tensile 
component of the stress at its tip, and this is completely determined by 
To ascertain the value of 7’ for which the crack will spread we follow 
Griffith (1920), and consider the energy of a crack of length 2¢ under a 


Cleavage of Metal Single Crystals 603 


uniform tension 7’. The elastic energy is 47B,57c2, where By. depends 
on the elastic constants of the material, and so, on adding the surface 
energy, the total energy of the crack is 

W = —4nB,.T?c? + 4ye. 
The condition dW /de=0 for the critical size of the crack, gives 


Toll = (4y/mBy)¥®, 2 . 2 1... (16) 
On substituting (15), we see that our cracks will just be able to extend if 
oC? + bu,x/2rhell? = (4y/7B,,)"? ; fe eM am CLIT) 


but if the left-hand side of this equation is less than the right-hand side 
growth of the crack will increase the energy of the system and so can- 
not occur. The left-hand side is a minimum when c=by,%/2zho,, 
and on substituting this value in (17) the condition that the crack can 
grow at every stage becomes 


Birnie yl Dame cerns nie (08s 1 eae (LS) 
Finally eliminating x between eqns. (12) and (18), we obtain 
0 0,L TF (y/7By){log (1449/05) 43 3}, : : : -) s (19) 


where 6=1+ ,/,; this equation may be regarded as the more precise 
form of eqn. (5). 

The constant 6 may be shown to be zero for an isotropic material ; 
generally it will be small compared with the logarithmic term and so may 
be neglected. Then yp, appears only in the argument of a logarithm and 
it is sufficient to know it in order of magnitude, which will be that of an 
elastic modulus. Thus the elastic constant of the material enter eqn. (19) 
essential only through B,,. The general relation between 6,, and the 
elastic constants is fairly complicated, though 6,, may always be 
calculated numerically by a straightforward procedure (Stroh 1958). For 
hexagonal metals 


Bog = — (811833 — 813")811_ | Y (Py + + Po *); (20) 
where p, and p, are the two roots with positive imaginary part of 
P*(8y1? — 8497) + "(811844 + 2511513 — 2513”) + S11533 — 513° = 9. (21) 


In these expressions s,,y are the elastic coefficients referred to the con- 
ventional axes, i.e. with the x, axis along the hexagonal axis, and not to 
the axes used elsewhere in this paper. 


§ 5. COMPARISON WITH EXPERIMENT 

In order to put eqn. (19) in a form more suitable for comparison with the 
experimental results we make a few minor modifications. If the crystal 
has diameter D the greatest value of L which we can obtain is D/cos x 
where y is the angle the specimen axis makes with the slip plane; since 
large values of L give the weakest crystals it is this value of L which 
will determine the strength. Also, in the previous section it was assumed 
that the dislocations were free to move under any stress no matter how 
small; in fact there will be a resistance to the motion, due to intersecting 
dislocations, impurities, the Peierls-Nabarro force etc. We may allow 
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for this by taking the effective shear stress on a dislocation to be not o, 
but o,—o, where o, is a constant stress which measures the resistance to 
the motion. 
Then eqn. (19) may be written 
OO = ROOK Vidg, toa ee ee (22) 
where 
k=(y/aBogD) log (u,9/o,). . - - + + + (23) 
If we neglect a small variation due to the logarithm we may regard k as a 
constant and independent of the crystal orientation. For a given crystal 
orientation o, and a, are readily calculated from the applied tensile stress 
and so (22) expresses the variation of the crystal strength with orien- 
tation. The comparison with experiment can most conveniently be made 


Fig. 4 
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Variation of strength with orientation. Experimental values for zine due to 
Deruyttere and Greenough (1956). 


by plotting o, against cos y/c,. Deruyttere and Greenough’s results (1956) 
for zinc shown in fig. 4 agree with the linear relation predicted by eqn. (22). 
Similar results for bismuth have been obtained by Vickers (1958) and will 
be published elsewhere. From the slope of the line in fig. 4 we obtain a 
value of k equal to 6-2 x 103 (g/mm+?)? which is in reasonable agreement with 
the value 8-1 x 103 (g/mm?)? calculated from eqn. (23). 

The intercept in fig. 4 gives a value of o, of 55g/mm2?. The axis corre- 
sponds to crystals with their basal planes normal to the stress axis and for 
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such crystals the ductility tends to zero. Crystals in other orientations 
have an observable ductility and for these the flow stress was found to be 
independent of the orientation; the initial flow stress of Deruyttere and 
Greenough’s crystals was 55 g/mm?, agreeing with the value of oy. 

When the basal plane is nearly perpendicular to the stress axis ( y= 90°) 
the resolved shear stress o, is small compared with the tensile stress; or, 
since o, must be at least equal to co, before the dislocations can move and 
fracture can occur, the strength of the crystal will be large in these orien- 
tations. This is in accord with the observation of Gilman (1957) who 
obtained a fracture stress as high as 4500¢/mm2 for a zine crystal with 
x=89°. In this orientation oc, will be large and eqn. (22) reduces to the 
simple condition that o,= 09, or fracture occurs as soon as the dislocations 
can move. Using the value of o)=55g/mm obtained above we estimate 
the strength when y=89° to be 3000g/mm?. Since this value depends 
very sensitively on the orientation this is as good agreement with Gilman’s 
result as we can expect (y=89°15’ would give exact agreement); it 
shows, however, that our model is able to account for the high stresses he 
observes. 

We are also able to understand from our model the variation of ductility 
with orientation observed by Deruyttere and Greenough (1956). We 
have seen that when y~90°, cleavage occurs as soon as the dislocations 
can move, so the ductility tends to zero in this orientation. For other 
values of y, plastic flow will occur as soon as the applied stress is such that 
o,=0 , but according to eqn. (22), cleavage will not occur until o, and the 
applied tensile load have increased; hence plastic flow must take place to 
work-harden the metal before fracture can occur, and a finite ductility is 
observed. 


§ 6. PoLYCRYSTALS 


If the model which has just been developed is accepted for single 
crystals, it is natural to ask whether the same mechanism can occur in 
polycrystals. However, this mechanism becomes more difficult in a poly- 
crystal since the constraint of the surrounding grains will oppose the 
movement of a subgrain boundary across a grain; while, on the other 
hand, it is the very presence of the grain boundaries that provides the 
obstacles which make a piled-up group mechanism possible. Thus it may 
well be that different mechanisms operate in single and polycrystals. It 
might be possible to decide between these two mechanisms for polycrystals 
by considering the variation of fracture strength o with grain diameter d. 
On the model of this paper we obtain (taking the length L of eqn. (19) 


proportional to d) 
‘eae 08 <3 ee C9) 


this law has also been proposed recently by Cottreil (1958) from a different 
model, The model based on piled-up groups of dislocations gives 


(Gees tag ts 82 ue ye ne ee (25) 
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Equations (24) and (25) agree provided the constant stress oy is small; 
this is the case for zinc (Greenwood and Quarrel 1954). However in the 
case of magnesium the experimental values (Hauser et al. 1956) seem to fit 
eqn. (25) better than eqn. (24) suggesting that piled-up groups are 
responsible for fracture here. 
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ABSTRACT 


The ligand field theory is applied to the hydrated manganous salts. The 
splitting of the degenerate levels of *A,g and ‘Hg is first investigated, and 
then the intensities of the transitions from the ground sextet ®Aig to these 
levels are calculated. The properties of the odd vibrations responsible for 
these transitions are also discussed. It is found that the calculated oscillator 
strengths are in good agreement with experiment if an appropriate model is 
used for the vibrations, 


§ 1. LyTRODUCTION 


Ir isnow generally admitted that the optical absorption spectra of inorganic 
complex ions are interpreted as transitions between the various d” or f” 
electronic configurations of the metallic ions involved, and, so far as the 
so-called ionic complexes are concerned, general features of these spectra 
are explained fairly well by the ligand field theory (Moffit and Ballhausen 
1956, Jorgensen 1957). The probabilities of these transitions are very 
small because initial and final states have the same parity in the case of 
free ions. Though it is well established qualitatively that these tran- 
sitions are electric dipole ones caused by the odd vibrations of the com- 
plexes Van (Vleck 1937, Broer et al. 1945, Satten 1957), as to the quanti- 
tative calculation of the intensities there is only the recent work of Liehr 
and Ballhausen (1957) for Ti (III) and Cu (II) complexes which can be 
treated as one-electron systems. 

In various complexes of the first transition group elements, Mn?* (3d°) 
complexes are rather special because all absorption lines correspond. to 
spin forbidden intersystem transitions and the intensities are quite small. 
This is due to the fact that the ground state ®A,,¢ is the only sextet obtained 
from the 3d® configuration. To explain the spectra, therefore, it is 
necessary to take the spin-orbit interaction into account. There are two 
approaches in the calculation based on the ligand field theory. In the 
weak field scheme, energy levels of the free ions are first determined by taking 
account of the Coulomb interactions between electrons, and the ligand 
field is regarded as giving rise to the mixing of these free ion states and to 
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cause the change of level structure. Orgel’s calculation (Orgel 1952) is 
based on this point of view. In most part of the present paper, how- 
ever, we shall adopt the so-called strong field scheme extensively developed 
by Tanabe and Sugano (1954). According to this model, the field acting 
on each electron is considered to be of cubic symmetry. The fivefold 
orbital degeneracy of d-electrons is lifted in such a field, and we get the 
lower triplet de orbits of T,, symmetry and the upper doublet dy orbits 
of Eg symmetry. After accommodating the electrons in these orbits, we 
introduce the Coulomb interactions. The energy differences between de 
and dy orbits, usually denoted by A or 10g, is considered to express the 
measure of the strength of the ligand field. Both approaches eventually 
lead to the same result. Figure 1 shows the positions of the lowest quartet 
levels as functions of the parameter A in the case of cubic ligand field. 


8s 


It is easily seen that the transitions in which electrons jump from de to 
dy or vice versa correspond to the transitions between curves with 
different slopes in this Orgel diagram. Such transitions are usually 
observed as broad bands or complex patterns consisting of numerous 
lines. ‘This may be interpreted as the result of excitation of even vibrations 
due to the change of ionic radius, because the charge cloud makes appreci- 
able change in such electronic transitions. 

To avoid inessential complexity, therefore, we shall restrict our attention 
to the transitions from ground state ®A,, to the levels of 4A,. and 4K, 
which are degenerate in the cubic field and henceforth denoted by 
{*A,,; *Eg}. These are transitions between the dédy* configuration, and 
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are observed as a narrow band in the spectra obtained at not very low 
temperatures (Jorgensen 1954, Holmes and McClure 1957). Recent 
experiment carried out by Pappalardo (1957) at very low temperatures 
has revealed the existence of fine structure of this band, which gives us 
some information about the detailed mechanism of the transitions. The 
aim of the present calculation is to clarify the mechanism by analysing 
this fine structure quantitatively. 


§ 2. Wave FuNcTIONS IN THE Cask or CuBIC FreLp 

We shall start from the Mn?* ion placed in a cubic field. Explicit forms 
of the wave functions can be easily written out by following Tanabe and 
Sugano’s (1954) scheme. Let the three de and the two dy functions be 
denoted by €, , € and u, v. Angular variations of these functions are of 
the forms yz2/r?, za/r?, xy/r?, (32% —r?)/r? and (a? — y?)/r? respectively. Then 
the wave functions of the states under consideration are, apart from the 
closed shell configuration, expressed as follows. 

SBA ig: de'(*Ayg)dy*(3Agg) ground states 


eA OF2)) = encwos i, hE ee Pe (251) 
PAyg: dF (*Agg)dy*(PAgg) 
[*A12(3/2)) = {3([Enlue] + [Enlav]) 
— 2([Enluv] + [Enluv]+ [EnZuv])}/4/380  . . |. (2.2) 
"He: deb(?Hg)dy*(PAgg) and de?(*Az¢)dy*(1Kg) 


*H_0(3/2)) = 3  (2LEndue] — [Enbww] — [EHouv]} 

+ 5 (enkun] - [ena PE DC 
[He 3/2)) =  (LEnduo] — [éilun]} 

+ 5 {lene — [éntva. ae Bare @3'5) 


For simplicity, here we mention only the states with maximum J, 
values which are shown in the parentheses, and each of the square bracket 
notations stands for a normalized Slater determinant obtained by accom- 
modating electrons into the orbitals mentioned in such a way that the 
spin of electrons in the orbitals without upper bar is up and with bar is 
down. The coefficients « and f are 


eae onde B= (3/7). wees ja, ee -7_(2.4) 

If we choose « = (3/7)#2 and B= —(4/7)"”, we get the wave functions of the 

upper ‘Eg states from 4D of the free ion which we will denote by 4He* in 
the following. 

The states 4A,, and *Eg and 4H,® have the same energy 10B+5C 

measured from the ground state ®A,., where B=F,—5F, and C=3oF, 


P.M. 2T 
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are Racah’s parameters and are considered to be measures of the Coulomb 
repulsion between electrons. If we adjust these parameters so that the 
calculated energy values of 4A,g, 4T,¢ and 4T,¢ agree with the experimental 
ones 25000, 18800 and 23000 cm respectively (Holmes and McClure 
1957), it turns out to be necessary to put C/B=4 and A/10~ B~800cm™. 

For the sake of calculation in the following sections, we shall here 
determine the wave functions of the 4T,, states. By solving the three- 
- dimensional secular equation of the energy matrix for 4T,g states given by 
Tanabe and Sugano, with the values B=800cm+, C=3200cm"}, 
A=8000cm-!, we can determine the coefficients A,,, ,, and v, of the wave 
functions 


[T, o"2(3/2)) =A, [EnoZu] + a {[EELuv] + [nAlur]}+v,[Enuvd] (2.5) 


as follows. 
Table 1. 


Energy in cem-t 
measured Ay Len 
from ®A,g 


18 950 0-977 0-144 
34 200 0-040 0-597 
42 700 — 0-208 0-788 


Two other components |*T,g"”) and |#T, 2”) of each triply degenerate 
orbital state are easily obtained by appropriate transformation of the 
coordinate axes. 

In some calculations, it is more convenient to express the wave functions 
in terms of the weak field language. Since *A,¢, and the lower ‘Kg states 
come from the 4G states of the free ion and they do not couple with other 
multiplets through the cubic field, we can express the wave functions of 
these states by linear combinations of the M, eigenfunctions |M,) of the 
4G states. The explicit forms are as follows :— 


[?Ai2) = (7/12)12| O>+ (5/24)14] to): | —4)} = 1b) 
[AH gM) = (5/12)"2] Oy —(7/24)"24| 4) + ]—4y} | \(2ahay 
[He = (2p ys [oak en 


§ 3. SPLITTING OF THE LEVELS 


The level structure shown in fig. 1 is obtained by neglecting the non- 
cubic components of the ligand field and the spin-dependent forces. If 
we take the effects of these forces into account, the degeneracy of the 
energy levels may be partly removed, and it gives rise to fine structure in 
the absorption spectra. In this section, we shall consider the effects of 
these perturbations on the energy levels under consideration. 
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Recently it has been shown that an appreciable amount of covalencies 
exist even in the so-called ionic complexes (Owen 1955). This effect is 
especially noticeable for dy electrons because the overlap of dy orbitals 
with the ligands is larger than that of the de orbitals. We shall, therefore, 
examine the influence of the difference of covalencies between dy and de 
orbitals on the splitting of the energy level {4A,,; 4Eg}. 


3.1. Spin—Orbit and Spin—Spin Interactions 
From eqns. (2.6) and (2.7) it is easily seen that the spin-orbit inter- 
action, which is expressed by AL.S for the (20+1)(28+1) states with 
given LZ and S, has no matrix element within each of the manifolds °A, ¢(8S) 
and {*A,g; 4Eg}(4G). Therefore, the changes of energy by this spin-orbit 
interaction are given by the second order perturbation terms : 


Kale d1.55)P 
AE.» ‘ 


Since AZ,,~10+cem- and ¢€~300cm~, splittings 6#, have the order of 
magnitude 10'cm~!, which can hardly be detected by the experiments. 
Usually the first order effect of the spin-spin interaction has the same 
order of magnitude as the second order one of the spin-orbit interaction. 
In an (LZ, 8) term, the spin-spin interaction is given by the expression 
(Pryce 1950 a, b) 
—p{(L.S)?+4(L.S)—$L(L4+1)S(S4+1)}. 2. (8.2) 


8H, ~> (3.1) 
b 


From the spectra of free Mn?* ion (Moore 1952), the coefficient is found to 
be p=0-334cem-!. Though the operator p(L.S)? has matrix elements 
between the twelve states of {*A,,; *Eg}, splittings caused by this pertur- 
bation are shown to be at most of order 25p~8cem. 


3.2. Effects of Covalency 


We shall consider the effects of covalency only for the dy electrons. In 
the simple molecular orbital treatment, this effect is taken into account 
by describing each of the electron orbitals as a proper linear combination 
of the dy and the ligand orbitals: 


cos 6|3dy) +sin O|c) , Le eek 


where |c) stands for an appropriate linear combination of the o-orbitals of 
the ligands. 

Since the greater part of the contribution of the interelectronic Coulomb 
repulsion comes from the charge clouds in the proximity of the nucleus of 
the central atom, it may be a reasonable approximation in the evaluation 
of the Coulomb and exchange integrals to replace the radial part of the dy 
orbitals by cos@ times of the original one and neglect the contribution 
from the second term of (3.3). In this approximation, Tanabe and 


2T2 
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Sugano’s energy matrices of the relevant levels (Tanabe and Sugano 1954, 
I, p. 764) are modified as follows: 


*Be(D, 4G) 


de3(?Kg)dy?(FAgg) 9Byt+3C )+4B' +20 OV GNsy 
de®(*Ayg)dy?("Kg) — 24/3B' 6B +30. 8B 220" 
Gea Me) tL) ee Ye 
4A o(4G): de® (*Agg)dy*(FAgg) 10B’ +50", 


(3.4) 
where energy of the ground state ®A,g is chosen as zero, and B’ = B, cos? 8, 
B" = B, cos! 0, C’ =C, cos? 6 and OC" = Cy cos* 6. 
Fig. 2 


ie) 005 O10 015 


ate 

Variation of the energy levels as functions of the covalency parameter 
«= 1 —cos?@ is shown in fig. 2. Even in the case of the free ion, By and Cy 
cannot be determined so that all calculated spectral lines agree with experi- 
ment. Hence, we have assumed C,/B,=4 and used B, = 900 cm obtained 
from the 4G<—+ ®S separation 26 900 cm of the free ion in drawing fig. 2. 
From these curves, the covalency parameter ¢« is found to be about 
0-08 ~0-10, and, if our assumptions are not wrong, the 4A,¢ level is lower 
than the 4H, doubly degenerate one by about 100 em. 

Since B~ 800 cmt obtained in § 2 is regarded as some average of By, B’ 
and 6", B’~800cm-! is a quite reasonable value. However, it is to be 
noted that this argument is valid only for the choice of the same ratio 
C) DiC aD 

3.3. Distorted Field 

Since the crystal structure is known for only a few of the manganous 
salts, it is impossible to predict the level structure for each salt quanti- 
tatively from the theoretical point of view. However, it can be shown 
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group-theoretically that the trigonal field does not remove the two-fold 
degeneracy of the E levels. In this sense the trigonal distortion is 
rather special. We shall, therefore, treat here the effects of the rhombic 
distortion and see how it affects each of the three orbital states of our 
{*A,¢; *Hg} level. 


Let us assume the potential of the form 


Av= > NC A) ae eee eH, 4: 


where 


il 
Av(x,y,2)=A f2- . (x24 ”} +A, = — (224 — xt — yt + 12x7y? — 6y222 — 622?) 


+By = (2° — 92) + By > (— 622n8 + bye? + at 9), - (3.64) 
By Esse apa Arve+4yt{ /Zye - Ja maryo} 
+ Byr(Y¥q2+ Vo %)/4/24+ By(V2+V.y2 . . (3.68) 


Then we get the following non-vanishing matrix elements related to the 
states {* Aves! Eg}: 


(HE PIAV Age) = — 243 Bair (3.7) 
es Tah7 at 
«*E g|AV [tAag) = iA aie wa = (S38) 


where +A,¢ is the state obtained by the configuration de®(4Agg)dy?(tA,¢g) 
and comes from the 4F states of the free ion. In the strong field repre- 
sentation, the wave function of this state is expressed as follows: 


*Ave(3/2)) ={[énlua] + [Enlvo//2, . . . (3.9) 


and its energy measured from the ground ®A,g is 226+ 70 irrespective of 
the cubic field strength. So, if we apply the second order perturbation 
theory, the energy denominator common to both *Hg® and 4Hg® is 
12B+2C ~ 20B~ 16000cm-1. 

This result shows that the 4A,¢ state is affected neither by tetragonal 
nor rhombic distortion and that if the distortion is tetragonal (Bb, = 56,=0) 
only the 4E, state is depressed. For the rhombic field, if we carry out 
a partial diagonalization using first order perturbed wave functions, we 
can eliminate the matrix elements between the *Ay¢, and the {*H,®, 4K ,} 
states and get the following Hamiltonian matrix for the perturbed *K,g 


states : W W.Ww.* 
| uf? ab aby 
A Fi Nee! ail 
Soe Gi We hep) aa 
where W,= (A rel AV {Ee ) 


W.= (AAsg|AV FEZ > 
and AE , = E(*Agg) — (4H g®) = E(*Agg) — H(*Eg®). 
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The eigenstates are, therefore, linear combinations of perturbed 4E,® and 
4,2 states and their eigenvalues are 0 and —(|W,|/?+|W,|?)/A#4. Hence 
only one of them is depressed by the rhombic field. 

Although it is quite difficult to know the magnitudes of the parameters 
A and B even when the crystal structure is known, it may safely be con- 
cluded from work on the cobalt salts (Abragam and Pryce 1951) that the 
splittings due to the non-cubic components of the crystal field are much 
larger than those due to the effects mentioned in (1) and (2). 


§4. Opp VIBRATIONS 


In those hydrated salts where the crystal structure is known, the 
manganous ion is at a centre of symmetry. Under these conditions 
electric dipole transitions between the levels we are considering can only 
take place with the excitation of one quantum (or an odd number) of an 
odd vibration. 


Fig. 3 


Q; ( v3) 


For the present purposes, it is a good approximation to regard each of 
the ligand molecules as a mass point because its internal molecular 
vibrations are in general of much higher frequencies. Then the problem 
is reduced to the vibrations of an octahedral XY, molecule. It is well 
known that the normal vibrations of this molecule can be classified into 
the modes of the symmetries A,g(1), E¢(2), Tya(3), Tyu(3), Tag(3) and 
T.u(3), where figures in the parentheses denote the degeneracies of the 
vibrations. The frequencies of these vibrations are usually denoted by 
V1, Vg, Vg, Vg, Vs ANA vg(v3< v4) according to the above order. 

Normal coordinates of the even vibrations A,g, Eg and T,¢ are given, 
for instance, in Van Vleck’s paper (1939). For the present calculations, 
however, the normal coordinates of odd vibrations are necessary. Those of 
the T,, vibrations can be readily obtained by symmetry consideration as 


Qi3 = (X_+ X5— X3— X,)/2, Qu=(¥3+ Y.—Y,—Y,)/2 
and Qis=(2,+2,—Z,—Z;)/2. . . . . . . (4.1) 


It is easily seen that the x-components of the two T,, vibrations are 
expressed by proper linear combinations of (M,/M)"2X , (X,+X,)/2 and 
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(X,+X3+X;+X,)/2, where M and M, are the mass of a ligand and the 
central atom respectively. These must be orthogonal to each other and to 
the coordinate 


Qie= {(Mo/M)Xo+ (X1+X.+X5+X,+X,4+ X,)}/(6+ M,/M)¥2 (4.2) 


of the translational motion in x-direction. The following set of the co- 
ordinates 


Q7' ={4M?/(My+2M)(M,+ 6M)} PX, + {4.M?/(My+2M)(M,+ 6M)}2 
i x (Xq + Xq)— (My + 2M) /4(My + 6M)}17(X,4+ X,+X5+X_), (4.34) 
an 
Qio ={2My/(My + 2M) PX, — {Mo/2(My+2M)}2(X,+X,) (4.36) 


satisfies the conditions. In terms of these coordinates, the part of the 
vibrational potential energy which depends upon the x-components of the 
T,y coordinates may be expressed as 


PaO. 0. i) el ee ee (Led 
w= Or Q0)(4 4) (Gi, (4.4) 
Actual normal modes Q, and Q,) are given by two such ortho-normal 
linear combinations of Q,’ and Q,, that diagonalize this 2x 2 matrix: 


Q,=Q,' cos0+Q,, sin@ and Q,)=Q,'sin?—Q,,) cos@. (4.5) 


To get the explicit expression of @, and Qo, therefore, we must make some 
assumptions as to the forces between constituents of the XY, molecules. 

The simplest assumption which Liehr and Ballhausen (1957) adopted in 
their calculation is the so-called valence force model proposed by Wilson 
(Yost et al. 1934). In the following, for simplicity, we shall assume 
M,/M =3, which is the ratio of the mass of Mn atom to that of a water 
molecule. Then in this model, our matrix is given by 


f 4,60) 


where K is the X—Y force constant, K, the interaction constant for the 
effect of extension of the X-Y bond, and K,, denotes the constant for the 
bending of the valence bonds. 

Another typical assumption is the central field model. Let the mutual 
potential of the central X atom and a ligand Y be f(r), and that between 
two ligands be g(r). Then, for M,/M=3 as before, we get 


Oy eye Se DE Oe Pare 
39 + 30R 203K 58! | 
2 : (4.6) 
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where R is the equilibrium X—Y distance, and 
f" =(@f/dr?),_»,f' =(dfldr),_p and g” =(d’9/dr?),_ jor: 


The interaction between three pairs of opposite ligands, 39(2h), has been 
neglected in deriving the matrix (4-65). 

In the case of aquo-complexes, covalency of the forces is expected to be 
very small, so in the former model K —2K, is much larger than K,. If 
we adopt the latter model and regard the six water molecules as dipoles of 
moment , we may put g(r) = (3/2)u2r-*. Since the distance 1/2 & between 
two adjacent waters is larger than twice of the ionic radius of oxygen atom, 
exchange repulsion between water molecules is neglected in this expression. 
Then, using the equilibrium condition f’(R)+ 4/2%g'(,/2R)=0, we get 
f’ =a’ and the matrix (4-66) is reduced as 


63 ld 6 ” 

eae eee 

| 20 20 | 4.60’) 
Ta 8 pete 8 — 
20 4 SW CeTEAT: 


Thus, in both models, the off-diagonal elements are found to be quite 
small compared with the differences of diagonal elements. It is, there- 
fore, reasonable to assume that Q,’ and Q,, express the actual normal 
modes to a good approximation. In fact, the choice of normal modes by 
Liehr and Ballhausen corresponds to the very small 6 of the expression 
(4.5). 

In this approximation, the frequencies v, and v, are given by 


ise 
= Ka (valence force) (4.7a) ° 
63, 
307 (central force) (4.76) 
(5 e 16e= 
3 (K—2K,)+ TR (valence force) (4.8a) 
vo 
Dees aon 
Eu T 399 > (central force) (4.8) 


where A;=47*»,?. On the other hand, frequency of the To vibration is 
given by : 


2K, (valence force) (4.9a) 


3 
qe (central force) (4.96) 
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So the ratio ve/vs= (Ag/Ag)"? is 
Ve ye (valence force) (4.10a) 
Vg | 0°49. (central force) (4.106) 


For the sake of comparison, we shall give here the expressions for the 
frequencies ,, vp, v; of the even vibrations as 


M,- K+10K, (valence force) 
f'+4g" (central force) 
(K-—2K, (valence force) 
Md, = , 3 
Rae aro. (central force) (4.11) 
|) — aR 
(4K, (valence force) 
MA, = ’ 
D> n" aS , 
7 + oF f (central force) | 


The perturbation responsible for the mixing of odd states is given by 


15 
V0, — VS 0Vi0Q)y.5. = 2 412) 
si 


where V is the ligand field acting on the 3d electrons. This may be 
expressed in terms of the cartesian coordinates of the ligands as 


= Vi,(Xo+X3+X;+X,—4Xo)+ Ve Xa 2X 6) VaQ13+ ies, 
(4.13) 


where Vije= OP (x;, Yi, 2;) etc., . : ° . . (4.14) 
a 


and 2;, y; and z, are the coordinates of 7-th electron. V,’s in (4.12) are 
easily derived from J, etc. by using the relation between Q’s and X’s (4.3), 
or 


6 + 6 
10 6 
X,+X,—2X)= = ie: . . 5 . (4.15 b) 
) 8 : 5 
y 7 = Ve and ope (; Vigt 2V ss) J ° 5 (4.16) 


We shall assume here that the field V is due to the six point dipoles of 
moment p vibrating about the equilibrium positions (+ R, 0, 0), (0, + R, 0) 
and (0,0, +R) without changing their directions given by a(e 1, 0,0); 
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(0, +1, 0) and (0, 0, +1) respectively. Then v,,, Ya, and v3, are written 
explicitly as follows: 


M1g= 8% q+ — ( 4/6 Yu) — /10 Y55) 
- a io 7 
Oa eg 
C5 i 6 Ys + 5(4/6 Y = 4/10 Y 93) 


| 
2 y JF (4.17) 
ad Ae JF Yat : 5 ss) 
10 2 2 
= - 2 (5 Yo + Ji Si 
1B Pgs [Bene Jaa) 


where 4g (—Y" 4+ Y,-) epr! e 
Yin = Jaa (iL, m= Iss, 5, ae F (4.18) 


§ 5. CALCULATION OF THE INTENSITIES 


For the sake of brevity, here we ignore the effects of the perturbations 
discussed in § 3, and use the wave functions given in § 2 as the unperturbed 
ones. Since electric dipole transitions between ground states and various 
3d° excited states are highly forbidden because of the same parity and 
different multiplicities of the initial and final state, to compute the tran- 
sition probabilities, we must take the change of wave functions up to the 
second order taking account of H,, the interaction with odd vibrations, 
and H,,=¢€>1,.s, the spin-orbit interaction, as perturbations. Then 
transition probability between a pair of such perturbed states |®g :7) and 
|4f:n’) is proportional to the square modulus of the matrix element 
<*f ;n'|Pl'g ;n) of the electric dipole moment vector P= Se r,, where 
each of n and n’ stands for a set of the vibrational quantum numbers. 
Bearing in mind that the ground state ®A,¢g is the only sextet obtained 
from 3d° configuration and that H,, has no matrix element between 
SA, ¢(=%S) and {#A,¢; 4Eg} originated from 4G, we get 


(fs n'IPlgs m) 


Bk (L, HE ed i’) 
oy Chm Phes wu; n'|H tk; n) 4h; nlH.,|°93 2) 
ee me (EL, deh rH, dy) 


Mv Ee, By 
i > Gir (mW ad bale nu; n' |H,|Sv ; n'y <p: n'|H,|69 ; n) 
ie (H,—E \(#,-—£ ) : 
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The intermediate states denoted by Greek letters are of odd parity, and 
lie far beyond the 3d° states. The even intermediate states 44 which 
couple with the ground sextet through the spin-orbit interaction are the 
*T,g states mentioned in $2. So the energy denominators of the third 
and fourth terms are much larger than those of the first and second terms, 
We may, therefore, neglect the third and fourth terms in the following 
calculations. 

The remaining terms are still unmanageable because of the many un- 
known intermediate states with odd parity. Main contribution will, 
however, come from the states in which one of the five 3d electrons is 
promoted to the 4p orbits owing to the large overlap of the wave functions. 
It will, therefore, be a good approximation to replace the energy differences 
H,—H,~H,—E, by their appropriate mean value AZ. Then the sum 
over the odd intermediate states can be carried out according to the closure 
properties of the eigenfunctions 

d< flAlm)<m|Bli) = (f|ABIi). 
Thus we get _ 
Cf n'PDas my = sy PP te 
k 
Let us first consider ie patie element <*k|H,,|°7>. The °S states 
couple only with 4P states through the spin-orbit interaction H,, and the 


$0? 
4P states are mixed in the *T,, states mentioned in §2 in such a way as 


(5.1) 


(4P_(M,)|*Tie”(M,) = re {y/2(rn Tt Vn) a ioe oe wre (5.2) 


where t=2, y, z. Using this relation and the Racah coefficients, we can 
determine the matrix elements of H,,, which are shown in table 2, where 


On = AV Ant Yn) — Bn} /5. 
Table 2. Matrix elements of H,, 


$0? 


IT, o(nex)(3) | AT, glre)(4) | 4Tg(ra)(8) | Pya(mvd(4) | #P,o(r2(8) | #T e204) 
*Ag(8) — /$oni | When 
®Ae(8) —~V/ sont Vion V Zant 
SA j2($) —oni|2 on/|2 a/ dant 
SAig(—#) — v¥ 30nt/2 V 30n/2 
°Aje(—$) 
SA 1g(—$) 


Matrix elements of PH, can be reduced to 
n'|PH,|*k; n)= x (n'1Q;|\r) AfIP Vth > 
= Son IQ; fn) CF 2 Pee. 3),(r,)|*k>, +. (5.3) 


where use has been made of the fact that p aod v, are both odd functions 
of the one-electron coordinates r,. Since our interest lies in the absorption 
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spectra at very low temperatures, we may assume that |n) expresses the 
state of no vibration. Then |n’) must be the state in which only j-th 
vibration is excited; n;=1. For such vibrational states 

(n'|Q;\n) = /h[2Me), Sia salen 
where w, is the angular frequency of the vibration associated with the 
normal coordinate Q;. Thus one gets 


(ifs wIPHBei ny=%, for AZ plenleth)-- (6) 


Since p and v, are independent of the spin variables, non-vanishing matrix 
elements of (5.5) are obtained only between the states with same M,, and 
are independent of the value M,. Using the wave functions (2.2), (2.3) 
and (2.5) for M,=3/2, one can reduce the matrix elements into the sum 
of one-electron integrals. Table 3 shows those associated with the 
z-component of the *T,¢ states. Matrix elements between {*A,g; 4H} and 
Te", 47, ."} can be easily derived by appropriate transformations of 
the coordinate axes. 


Table 3. (4f|[¥p(r,)o,(r,) PT) 


of *Tie("*) 
4A — (5/6)1 (An + rp) <o]p-v;| > 
4H) (An te Vn) <v|p-v;| > / 4/42 
3K — (3/14)? (A, + yp) <u p-r;|6> —4 pn <Elp-v; [> //7 


Ultimately, it is found that all matrix elements can be expressed in 
terms of the following nine one-electron integrals: 


Cele ey=3, 15 (5 (Fa) + Taam) + 3g Gm) 


(Sleroylu) = V10(~ 358) + a1 Bae Uae 5)) oe 


R2/ * 91\R4/ * T1\Re 
élrisuley = 73 S 7m 


Zle> ye) = sy (2 Se (i (i | 
Cawley = IP (5 a vy pa) ) : (5.6) 
CClesavlr) =a (F («ae 

yah) = 75 — 5) — 7 Ge) Th) 


= DES 26 ,r4 9,78 
chal) =4, (5 CFs) ~ at Ge) +) ) 


(flersyln) =. | 
J 


In deriving these, use has been made of the expressions (4.16) and (4.18). 
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§ 6. Resunrs AnD ComPARISON witH EXPERIMENTS 


To work out the intensities or oscillator strengths numerically, it is 
necessary to evaluate the values of ¢(7?/R?), ¢r4/R*) and {r§/R8), Analy- 
tical functions of Slater type (Slater 1930) and numerical Hartree-Fock 


Table 4. 


ere R2»y <n| > PAS 


Hartree—Fock 0-075 0-012 0-004. 
Slater 0-21 0-07 0-04. 


functions (Hartree 1955) give entirely different results. Though the 
latter are of course more reliable in the case of the free ion, if we use these 
values and R=2-35x10-8cm and w=4~x10-¥e.8.u.cm, it is found 
impossible to obtain correct order of magnitude for the intensities. In 
fact, if one calculates A by using the same model, the relationship 

. 25 eu 74 

AS 1D¢= = oa \ pi)? eid ee re Gi) 
gives 


S = 0-054 


for the same values of R, » and A=8000cm~-!. Since the wave functions 
may be considerably distorted in the crystal, it is quite difficult to estimate 
correct <r”) values. In effect, this difficulty is related to the validity of 
the ligand field theory itself. Several attempts of non-empirical calcu- 
lation of A reveal this situation (Kleiner 1952, Tanabe and Sugano 1956). 
We shall, therefore, treat these as empirical parameters, and tentatively 
assume the following two sets of values: 


Case I: 

(77/R?)=0-17, <7r*/R*)=0-05 and <r°/R®)=0-025, (6.2) 
Case IT: 

(r2/R2) =0-20, r4/R*)=0-05 and (r°/R8) =0-01, (6.26) 


The spin-orbit interaction parameter is found from the free ion spectra 
as C=300cm~!. The mean energy difference between odd parity states 
AE is assumed 120000cm- from the spectral data of the 3d*4p con- 
figuration states (Moore 1952). It is to be noted that if we consider only 
3d44p configuration in intermediate states the terms containing ¢7°/R®) 
never appear in (5.6). For this reason we have assumed rather small 
values for (r*/R®) instead of using different A# value for these terms. 

For the vibration of the Mn2+ aquo-complexes, there are no experi- 
mental data of the frequencies. However, for several aquo-complexes of 
other iron group elements, frequency v, of the A,g vibration has been 
determined from the structure of the optical spectra (Schultz 1942). 
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From these values, one can presume v,=250~300cm- for hydrated 
manganous complexes. As is seen in §4, vy=(5/3)"?»,, so that one gets 
v4=320~400cm-1, Since the Jahn-Teller effect can give rise to the 
E,(v,) vibrations which appear as sharper peaks than those of the 
A,e (v,) vibration because of its weaker coupling with the environment, 
the frequency attributed to the A,g vibration by Schultz might be v, of 
the Eg vibration. However, we can see, from eqn. (4.11) that the two 
frequencies are very nearly equal in that K>K, and f’>g”. Hence this 
estimate of v, need not be altered. 

As to the frequencies vz and vg, no information has been obtained experi- 
mentally. Though Liehr and Ballhausen (1957) adopted the valence 
force model, the central force model would be more appropriate for the 
aquo-complexes. If we use the simple model discussed in § 4 and assume 
g(r) = (3/2) 2r-3 with w= 4 x 107!8e.s.u. em (this value makes allowance for 
the induced dipole moment on the water molecules), eqns. (4.7b) and 
(4.96) give 

¥,=1700en and %=8iem-. +) 3) aa VE 


The calculated oscillator strengths thus obtained are shown in table 5. 
These are the values averaged over the polarization direction of the 
incident light. The intensities of the transition accompanying excitation 
of the high frequency T,, vibrations are in general small because of the 
factor (5.4). 

Table 5. Calculated oscillator strengths 


Calculated oscillator strengths 


Electronic Associated 

transition vibration Coney aT 
SAie>tAie Ae) eel (me 3 x 10-8 
SAig>4 1g Tiu(v3) 6 4 
PAye>*Aye Tou(v 6) 7 a 
SA g> tHe) ata) 0:04 0-10 
8A g> tHe? Tul) 4 5 
SA g> Kg Tou(v.) 0:8 0-8 
®A go *Eg® Tul) 0-04 0-10 
SA, g> “Hy” Tyu(v3) 2 5 
SA p> 4H, one) 0-8 0:8 


Magnetic dipole transition, if it exists, would give an intensity of the 
same order of magnitude. However, in our case, we need not consider it, 
because magnetic dipole transition takes place only between the states 
with the same Z and S and we cannot expect any perturbation which mixes 
the °S and 4G states appreciably with each other. 

Experimental spectra have been obtained for various manganous salts 
(Gielessen 1935, Pappalardo 1957). For most of these salts, however, it 
cannot be allowed to treat [Mn(H,0O),|?+ as an isolated molecule. We 
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can, therefore, compare our theoretical results only with the spectra of 
MnSiF,.6H,O obtained at 20°K by Pappalardo. 

He observed three comparatively intense absorption peaks in the 
region 24500~26000cm- of the wave number. The middle one at 
25 330cm™ is considered to be a double peak of separation 14em-!. So. 
one can regard these as four lines centred at 25420, 25337, 25323 and 
25075 cm respectively. We may ascribe the first and second of these 
lines to the transitions *A,g—>4A,¢ accompanying the excitation of lower 
Tyu(vs) and T,, vibration respectively, because the separation of these 
two peaks, 83cm™, is quite consistent with our estimate of vz and v,. 
Then the remaining two may be attributed to the transition °A,,—>4E,®) 
and ®A,g—>4E, with the excitation of the lower frequency T,, vibrations. 
Then some of the other lines can also be assigned as shown in table 6. It 
is, however, impossible to discover transitions corresponding to the faint 
lines found at 24946 and 24896cem~!. Furthermore, lines corresponding 
to the transitions °A,z,—>*A,¢+hv, and ®A,g> 4H, + hv, are not observed 
in spite of their comparatively large intensities predicted by our theory. 
The wave number corresponding to the transition ®A,g+hv,>‘*A,¢ is 
expected near 25170cm~!. The intensity of this transition is, however, 
quite small at 20°K because of the Boltzmann factor exp (—/hv,/kT’) ~ 0-003. 


Table 6 

Wave Observed ah ate eae: 
number oscillator Assignment ees 
in em-t strength Case | Cen, 

25 420 Gx 104 SA got Aig this 6 x 10-8 4x 10-8 

25 337 7 SA g>tAig thr, ff 4 

25 323 SA go> tHe” +hvg 4 5 

25 170 SA ge thve*Aye(?} 0-02 0-01 

25 075 3 SA g>*He™ + hv 4 5 

25 002 SA g>tHe +hv, 0-8 0-8 

24 946 a 

24 896 a 


§ 7. CONCLUSION 


Though some of the details are still unexplained, the general features 
of the observed spectra are explained and the calculated intensities give 
correct orders of magnitude. To make a decisive assignment, some more 
detailed experimental information, for instance crystal structure, direction 
of the polarization, etc., are awaited. 

It may be stated that the forces responsible for the vibrations can be 
described by central force model rather than by valence forces. 
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ABSTRACT 


The solution of the elastic equations is considered for the case in which 
the state of the solid is independent of one of the three Cartesian coordinates. 
The stresses due to a dislocation, a wall of parallel dislocations, and a crack in 
an arbitrary non-uniform stress field are obtained. The results hold for the 
most general anisotropy in which no symmetry elements of the crystal are 
assumed. 


§ 1. IyTRODUCTION 


KsHELBy et al. (1953) have developed the theory of anisotropic elasticity 
for a three dimensional state of stress in which the stress is independent 
of one of the Cartesian coordinates, and have applied this to find the stress 
field of a dislocation. In the present paper, which follows their treatment, 
the stresses due to a dislocation are treated more fully, and the interactions 
of dislocations considered; also, the stresses round a crack subjected 
to an arbitrary non-uniform applied stress are obtained. The object 
will be to present the results in a form which is, analytically, as simple 
as possible. It is hoped that, in applications of the theory, this will 
often allow of the properties of the system studied to be deduced without 
the need for numerical computation, and that when such computation 
is unavoidable, as when definite numerical values are required, the labour 
involved will be reduced to a minimum. For this purpose, the properties 
of a number of constants introduced in the theory and which are related 
to the elastic constants are investigated in some detail (§3). In §2 some 
general relations are considered ; most of these are given in Kshelby et al. 
but are included here, both for ease of reference, and so that the whole 
theory may be presented in a uniform notation. 


§ 2. GENERAL EQUATIONS 


The stresses o,; are related to the elastic displacements u;, by the 
equations 

O55 = Ci OU, OX, A hid pe ee e8) 

where i,j, k, J=1, 2, 3 and the convention of summing over a repeated 

Latin suffix is used. The elastic moduli ¢,;,,; have the symmetry properties 


ate! at pare Cele 
Cit = gina = Cite = Cee (2) 


+ Communicated by the Author. 


P.M. 2U 
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It is often convenient to replace the pairs of suffices (i, 7) and (h, 1) by 
single suffices M and N according to the scheme that 11 corresponds to 
1, 22 to 2, 33 to 3, 23 to 4, 31 to 5, and 12 to 6. Besides securing some 
brevity in writing, this has the advantage that it enables the elastic 
constants to be considered as the elements of a matrix, but has the 
disadvantage that it does not make the tensor character of the constants 
apparent. In the sequal either four suffix notation ¢;,,; or two suffix 
notation ¢,,y will be used according as to which is the more convenient 
in each cage. In addition to the symmetry relations (2) the elastic moduli 
are subjected to further restrictions which arise because the elastic energy 
density must be everywhere positive. If the elastic strains are 
COU [OC ered ean , €g = 0U,/Ox_ + Ouy/0x,, the energy density is 

Ben Carey ous 5 tins iene nee eee 
provided not all the e,, are zero. The condition for (3) is that the deter- 
minant |c¢,,, | and its principal minors of all orders should be positive. 

On substituting (1) in the equilibrium equations 

do;;/0x; =, Perey er 
we obtain 
C45 j40°U | Ox Ox, = 0. spent St 
Now we suppose that w, is independent of x, and, following Eshelby 
et al. take 
U,=A,f (41+ Xo), 45 tye Se area Om 
where f (z) is an analytic function of the complex variable z; (6) is a solution 
of eqns. (5) provided the constant vector A,, satisfies the equations 


(Cina + PCa po + Pons + DCion2)A y=. AD) 
Values of A;,,, not identically zero, can be found to satisfy these equations 
if p is a root of the sextic equation 


Cinna + PC Ke + PCat + Coke | = 9. ee 
Eshelby ef al, have proved that eqn. (8) has no real root, so that the roots 
occur in complex conjugate pairs. The three roots with positive imaginary 
part will be denoted by p,(«=1, 2, 3) with complex conjugates p,; the 
corresponding values of A), obtained from eqns. (7) are A,, and A,,. 
Summation over «, which of course is not a tensor suffix, and generally 
over Greek suffices will always be indicated explicitly. It will be assumed 
that the roots p, are all distinct; equal roots may be regarded as the 
limiting case of distinct roots. A general expression for the displacement 


may then be written 
Ui) A ned glen) o> Aneta ee) ee 
where 2, =2%,+ ),%o. : 
It is convenient to express the stresses in terms of a vector, the com- 


ponents of which will be denoted by ¢,; with the stresses independent 


of x3, eqns. (4) are satisfied identically if the stresses are derived from the 
¢; by the relations 


Oy = —09;/0%2, Fi. = 04/0. Perm arse. (KY) 
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These equations determine in terms of the ¢, all the components of the 
stress except o33. But since in the present case o, is linearly dependent 
on the other stress components (the equation showing this is just the 
condition, expressed in terms of the stresses, that the strain e; = u3/0x4 be 
zero), we may regard the problem as solved when the remaining five 
components of stress, or equivalently the functions ¢,, are determined. 
Since 


0, /0X, = 043 = — 0¢4/0X2, 5 Pier, been ated aL) 
¢, and ¢, may be expressed in terms of a single function y by 
$1 = —Ox/0x, $,=0x/0%,; sptopmeren ee MeL) 


however the theory takes on a more symmetrical form if this is not done, 
and we shall generally prefer not to introduce y. Eshelby et al. have 
shown that the resultant force acting across a curve C is A¢;, where Ad, is 
the change in ¢; on going along the curve. It may also be shown that the 
moment about the x axis of the forces acting across C is 


A (261 — 2% 142+ x). 
From eqns. (1), (9) and (10) we obtain 


$i Liat al@e) + SLinf lx)» i ees) 
where the three vectors L,, “aes defined oy 
Lig = (Cor + Palirea) A kas a ee UES, 
or alternatively, 
Tere Pa. Cura bonne) oh + ee 5) 


these expressions being equivalent by eqns. (7). For (11) to hold, we 
must have 

Fag t Patten =, Ce gh) Bea rs ered WILD 
a relation which also follows directly from (14) and (15). Use of eqn. (16) 
will simplify the expressions for the L;, and so we shall prefer to express 
our results in terms of the Z;, rather than the A;,,; as will be seen in §7 
the L,,, can be obtained in terms of the p, without first determining the 
A,. Nevertheless, in the general case, these expressions still remain 
rather cumbersome, and we proceed first to develop the theory as far as 
possible without introducing the explicit values of the A;,, or Lj,. 


§ 3. Propertigs or A,,, L;, AND SOME RELATED CONSTANTS 


From the vectors A,, and L;, we may form the 3x3 matrix whose 
elements are A,,L,,; this matrix we now prove has skew-Hermitian 


symmetry. For from (14) we have 

PA iplin a PpAigCinirA ke 1 Bp AspCisnePrA kee 
and from (15) : 

A pling Pa ae AiCnraA ko A iC n2P oA ka} 
subtracting these two equations, i 
A LialBp —B2) = AipCn tr} tat Bp AipCink2PcA nat (PpAigCiarrA ka 
+ Aisi n2PoA ta): ee nil T) 
202 
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Now each of the three terms on the right of (17) has Hermitian symmetry, 
and so the left side must also; that is 
A yp Lig( By — Pa) = A, whip Da P Dp) 
or since (p,—p,) has negative imaginary part te so is not zero 
Aplin = — Aialag, emg beater Gk) 
which is the required relation. 

Since A,,, 4,, form a set of six vectors in three dimensional space, 
not more than three of them can be linearly independent; the three 
vectors A,,, corresponding to the three roots p, with positive imaginary 
part, do however form such a linearly independent set. For suppose 
that this were not so and that a relation 


eA, =0 Wet so% Se eet. 


existed with not all the €, zero. Multiply eqn. (17) by E,£p and sum over 
x and B. The left hand side becomes 


DEA ipi(Pp =D, joa a5, DipPpepXAsab r= DAs MjoPaba =0 


on using (18) and (19); also after using (19) and its complex conjugate 
the only term that remains on the right-hand side is 


(38 pPoAip)einnal Yé,P 1A kx) = 0- 


But [C2], being a principal minor of the 6x6 matrix [c,,,], is by (3) 
positive definite; hence 


Ena ha = 0. bok clos goa oe 
From eqn. (7) we obtain 


Cai daPa A +(c ike + Sj271 26a Ata + Cana2baPa Aj.=9, 


or using (1 9) and (20) 
Ca ala A ka — 0 > 
again by (3) |¢1;, | #0, so that 
DIS eee => 0. . . . . . * (21) 


Now for a fixed value of k, (19), (20) and (21) may be regarded as a set of 
three equations in the three unknowns €,A4,,; the determinant of these 
equations is equal to (p,;—p.)(P2—ps3)(P3—P1)/P1P2P3, Which is not zero 
if all the p, are distinct ; hence the only solution is €,4,,=0. But A ty 18 
not identically zero cane so we must have €,=0, or no relation of the form 
(19) with non-zero €, is possible. 

Next we prove that the linear independence of the L,, follows from that 
of the A;,,; in the course of establishing this result some constants which 
will also be needed later are introduced. Again we suppose that the 
L,,, are not linearly independent so that we can write 


Dealiy = Mal fine Ae Rie 4 5) 
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with not all the €, zero. From (14) we have 
22E pAipligDy ME a. 2é pAipCinn >. Pa i€ ake + DE pAipCionad& athar (28) 
(od e4 B % 


and the left-hand side is seen to be zero on using (18) and the conjugate 
of (22). Also, from (15) we have 


Cin) Pa tA rate t Cine>Apaba =—>1,,6,=0.°. . .» (24) 


Now the determinant |¢;,;,| is by (3) positive and not zero, and hence 
the matrix [s’;,,] reciprocal to [¢,,;,,] exists. Then from (24) 


a ee kaba = — 8" iCute 2 faArw, 
and substituting this in (23), we obtain ‘ 
(dF e4 ip) Bul D> 0Aix) = 9, 2 et ae AS) 
where ; : 


Bin = Ciny2 — CiogrS’ njCju2- eect poe tec 0) 

The expression 
Cia 1121 C131 C12 
ree lesan | x Coi11 ©2121 2131 += Cati2 

C3111 3121 «©3131 Cs1l2 

Ciorr  Cjn01 ig31 Gale 
is easily verified by expanding the determinant on the right by its 4th 
row and 4th column. With 


Ay= |e] = Cy5 C55 C56 | > de eee) 


we have 


A,B o2= > 
Ci5 Con C55 C56 
Cig Cog C56 66 
Cyr Cia C15 = C16 
Es Cc @ @ Cc p 
A, B33 = OS Soa y gel Cae Peer (25) 


A 1Po3 = Ai P32 = 
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and the remaining elements of B, are zero. It follows from (3) that 


Boo > 0, and Pa,> 0. ee dr te 4) 
Further, consider the determinant 


BoA, is the cofactor of cy, in A, BsgA, of Cy2, and —Bo3A, of Cy4. Hence 
it follows from Jacobi’s theorem that 


Bo2B33 — Bos? = A/A, > 0. oe Ber ant 1) 
On account of (29) and (31) eqn. (25) can be true only if 
»¢,A,,=0 and »é74.5=0. PA AA Bark ly) 
The third equation, ~% a 
dé.A1.=9, co eh pene eee an oa} 


does not follow from (25); it may however be deduced in a similar way 
but starting from the equation formed from (15) in a manner analogous to 
that by which (23) was formed from (14). Then (32) and (33) show that 
the linear dependance of the L,, implies that of the A,, which has already 
been proved untrue ; hence the linear independence of the JL,, is established. 
In deriving (33) we would have to introduce constants y, analogous to 
8, and which will also be needed later. The non-zero components of y, are 
C1, Cig Cra C46 ] 


Agyu = a > 


As/33= 2k be UeW ges (OA) 


Asyi3 = Asysi = i ? 


where 


Ay= |Cjo,9| = Cog Cyq Cygg | - (SD) 
Cog 4g gg 
Also Vu > 0, Y33 = 0, and V11Y33 — Maas => AV aN. > 0. 5 (36) 
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Now having established that the L;, are linearly independent, we are 
able to introduce three vectors M,,, reciprocal to the L,,, and which are 
defined by 

M yi Lip=Sap- Oe Neel, gn eet ee PENT) 
From (37) the relation 


Lin M 4; = 8; . . . . ° . (38) 
follows. Also eqn. (18) may now be written ; 
>4,.M,i — = >4;.M,;- . . . . . (39) 


Since, for each «, eqn. (7) determine only the ratios of the components, 
A;, (and hence L;, and M,,;) contains an arbitrary factor independent 
of 7 but which may depend on «. However this factor does not appear in 
the product A;,M,; so that }A;,M,, is unambiguously defined. A related 


tensor which will play an important role in the theory is 
By =4i>(4,,M,;—4;,M,;) 5 . . . . . (40) 
B;; is clearly real, and being equal to the imaginary part of a skew-Hermitian 


matrix (cf. eqn. (39)) is itself symmetric. 


Finally we seek a relation between 5;; and L;,, M,, not involving the 


A,, explicitly. From eqns. (15) and (38) we have 
8; = >L,,M ef a 24 ie Dre aj Canz2A ka aaj? 
and from (14), 
eiaPa Maj = C21 DA aPa Mag + Girne XA bo M ay. 


From these two equations we obtain, on using (26) and (40), 


1 i 
2j (Qi sy Qi;) es —Pi7,Byj, ee re ee aL) 
where 
is = LiPo May. Siders Pee (Se) 
In the same way with 
Py a DS ye ea Catt. aire @er a ies (43) 
we obtain ; 
il 7 
5 (Piy — Pay) = Vir B y- ween 2 ieee (44) 


The matrices [P,,] and [Q,;] are reciprocals of one another, as is easily 
seen on forming their product and using (37) and (38). 

Since [f,,] and [y,,] are singular matrices, neither (41) nor (44) alone 
can be solved for B,;. However taken together (41) and (44) are sufficient 
to determine B,; and when numerical values are required it will usually 
be easier to solve them than first to determine the A,, and to use the 


definition (40). 
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§ 4. DISLOCATIONS 
Consider the following displacements, which are of the form (9): 


1 |e we 
Uy, = Darj  tae log pa Dor 2A rae log zZ,. . . . (45) 

Along a closed path encircling the x, axis, w;, changes by an amount 
b,= Zt pont A has) . . . . . (46) 


lo 


Also from (13) the stress functions corresponding to the displacements 
(45) are 


eee 1p. ea ee 
ee ; gag EE 3 ; : Re a YE 
d; ort dHiPa log ey 9 1 2a log Ry ( ) 
the change in 4; along a closed path about the x, axis is 
Ad; ee DLiDPa ats LiDa); toe es (48) 


and it has been noted (§2) that this change in ¢; just represents the 
resultant force Ff; acting. 

Thus in general (45) represents the displacements due to a dislocation 
with Burgers vector 6;, together with a line of body force F;, along the 3 
axis. Since the D, and D, constitute a set of six independent constants 
we may expect that it will always be possible to choose them to correspond 
to any given values of 6; and F;. We now determine these constants so 
that we have a pure dislocation without a line of force along its axis, that 
is we take F';=0. 

It is always possible to express D, in terms of new constants d; by the 
relation 

D,=HM 4 d;; 
then substituting this in (48) we see that F;=0 if d; is real.. Equations 
(46), (47) and (40) now show that the stress functions 


al ae 
$;= Foy a Pi aj log 25, + Lig Mas log 2.) d; eee PAC 8 
represent a dislocation with Burgers vector 
6,=B;; d; . . . . ° . . (50) 


for all real values of the vector d;. It will be seen later that the deter- 
minant |5;;| is not zero; hence given the value of 6, it is always possible 
to solve eqns. (50) for d;, so that (49) can represent a dislocation with 
arbitrary Burgers vector. However the stresses are more closely related 
to the vector d; than to 6; and so we shall usually prefer to express our 
resuits in terms of the former. (The nature of d; may perhaps be made 
clearer by reference to the isotropic case where d ; are the components of 
Gb for a screw dislocation and of Gb/(1—v) for an edge, with @ the rigidity 
and v Poisson’s ratio). The stresses derived from (49) by eqns. (10) are 


OF (1/427) > (LinD 2% 1M ,; “rh Lip fy MG) d,, ° (51) 
and C2 = (1/4ar) > (L,,2,7° as Ee eae) d,. 2 C s (52) 
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The elastic energy per unit length of the material lying inside a cylinder 
with generators parallel to the x, axis, and which meets the plane a, =0 
in a curve C, can quite generally be written 


U=3| U;o,;n; ds, eee ee ras (OD) 
/ C 


Taal 


where 7; are the direction cosines of the outward normal to C. On using 
(10) this becomes 


U=—3| u(ddy/ds)ds. ee ee ee oa) 


We apply (53) and (54) to calculate contribution to the energy of a 
dislocation from the material between radii R and ry (R>1r,) from the 23 
axis (fig. 1). This region may be made simply connected by making a 
cut in the plane x,=0 from x,=7,) to z,=R; in the cut region we take the 
displacement wu; to be single valued but discontinuous across the cut. On 
the surface of the cylinder radius R, (54) gives a contribution to the energy 


of 2ar 
as { u,(d¢,[d0) dd, 
0 


where we have taken z, =2,+ p,%_,=R(cos#+p,sin@). From (49) we see 
that d¢,/d@ is independent of R, and from (45) that u; is of the form 


u;=klog R+4 (8), 
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where k is a constant and g() is independent of R. Thus the part of U, 
involving £ is 


—iklogR [¢ (d¢,/d0) dd =0 


since ¢, is single valued; U, is thus independent of the radius k. From 
the surface ai the inner cylinder of radius rg we obtain a similar contri- 
bution which will be of the same magnitude but, since the direction of 
the normal n, is now reversed, opposite in sign. Thus the net contribution 
to the energy from the two curved surfaces is zero. 

The displacements on the two sides of the cut differ by b;, and the stress 
across the plane of the cut 7,=0 is, from (52), 


O49 = (1/4721) > (L,,M,;+ L; 1 ,;)d,; = d,| 2rr0e. 
Hence (53) gives the energy as 


U =(b,d,/47) IE dxe,|x, =(b,d,/47) log (R/ro); - - - (55) 
using (50) we have the alternate expression for the energy of the dis- 
location 

Uc (By d, dia) loo (Tita): ee tae 


Ui pnaet thea} 

Since U must be positive for all non-zero values of d;, we see that [B;,;] 
is positive definite. In particular the determinant |B;,| is positive and 
not zero; this establishes the statement made previously that eqn. (50) 
can always be solved to give d; in terms of 8;. 


§ 5. INTERACTION OF Two DISLOCATIONS 

The energy of interaction of two dislocations can most easily be found 
by the method due to Cottrell (1949) of forming one dislocation in the 
presence of the other. Consider a dislocation with Burgers vector b, 
lying along the x3; axis and suppose a second dislocation with Burgers 
vector 6;° is formed parallel to the first and through the point (X,, X5, 0). 
This dislocation may be formed by making a cut along the plane 7,= X, 
from x,=X, to x,= 0, and giving the surfaces of the cut a relative dis- 


placement 6°. The work done in this process against the stresses of the 
first dislocation is 


V =), | = oi) dry ; 
X, 
using (52) this becomes 
V = — (62/47) »> (L,,M,;logZ,+L,,M,;logZ,)d;%, . . (57) 


where now Z,=X,+p,X». Introducing polar coordinates (r, 4) (fig. 2) 
we have 


log Z, =logr + log (cos 6+ p, sin @). 
Then from (57) and (38) we obtain the radial force 
P,= —0V [er =b2d/2nr ; 2 eee 8) 
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by using (50) this may also be written in the symmetrical form 

diye = B; AMA, /2ar. . . . . . . (59) 
The expression for the tangential component of the force is rather more 
complicated and is given by 
Lev Oa s ( Mf p,cosd—sind  - ir. 22808 6—sin 5) 


a 00 ar “cos 0+p,sind § “4 cos0-+ 9, sind 
(60) 


a 


Fig. 2 


An important special case of the forces between dislocations is that of 
two like parallel edge dislocations in different slip planes; the well-known 
result of isotropic theory that the dislocations have an equilibrium position 
with the line joining them normal to the slip will not be true in general. 
For suppose the dislocations have Burgers vectors 


b,=(b, 0, 0) es er ete, Oey (61) 
so that the slip planes are parallel to the plane x,=0, and let the dis- 
locations pass through the origin and the point (0, x, 0), the component 
of force in the slip plane may be found from (58) and (60), or more simply 
by noting that this force is 


By = doin = F042 ; 
on using (51) and (38) this becomes 
P= —(b/4ra,) > (Lo,M yj + LoyM 43) dj = — b dy/27p. . (62) 


a 
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The value of d, used here is to be obtained by solving eqns. (50) with 6; 
given by (61); it is shown in §7, by consideration of a special case, that in 
general dy will not be zero, and so the dislocations will not be in equilibrium. 

It is easily seen from (52) that o,, changes sign as x, varies from — 
to + 00, so that it will vanish at an odd number of points; that is two 
parallel edge dislocations have an odd number of equilibrium positions. 
It is convenient to characterize these by the angle 6 between the line joining 
the dislocations and the slip plane, since, from the form of the forces (58) 
and (60), the equilibrium positions depend on @ but not on r. If, then, 
the dislocations at O and P in fig. 2 are in equilibrium, we can add any 
number of similar dislocations at different points along the line OP; every 
pair of dislocations, and hence the system as a whole, will be in equi- 
librium. Thus the values of 6 corresponding to the equilibrium of a 
single pair of dislocations define possible orientations in which dislocation 
walls may be formed; moreover it is easily seen that these are the only 
orientations in which the dislocations forming a wall will be in equilibrium. 
In particular we are led to the conclusion that a wall of edge dislocations 
normal to the slip plane (a simple tilt boundary) will not be possible 
unless 6=47 is an equilibrium position for two dislocations, which, we 
have seen, is not always the case. 

This last conclusion appears inconsistent with that obtained from the 
treatment of Frank (1950); in this the dislocations walls of least energy 
(and which are therefore the most stable) are obtained as those for which 
the material is unstrained at a distance from the walls. This gives the 
simple tilt boundary normal to the slip plane as the most stable arrange- 
ment which can be formed from edge dislocations all of the same kind; 
Frank’s arguments are very general and should apply quite independently 
of the crystal symmetry. We therefore consider, in the next section, the 
stresses due to a wall of dislocations and verify that the remote stress 
field leads to results in agreement with Frank. 


§ 6. DisLtocaTion WALLS 


Consider a wall of like parallel dislocations, with arbitrary Burgers 
vector 6;; the dislocations are to be equally spaced in the wall, the distance 
between neighbouring dislocations being h. We take 2, axis normal to 
the wall, and, as usual, 2, axis parallel to the dislocation lines. The 
stresses due to the wall can be found by summing the stresses (51) and (52) 
over all the dislocations in the wall. Thus the stresses will be given by 
two expressions of the form (51) and (52) but with z,-!=(#,+p,a.)7 
replaced by 


N 
lim t+ p,(t,—vh)\ = —— Got —— 
N>0oo 3 a s Pad 2 A 5 ph cot Pal (ay + PX). e (63) 


When «, is large (63) reduces to im/p,h; then from (51) and (38) we 
find 


O77 = 0, 
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and from (52), (41), (42) and (50) 
OD Ve (7/4h) > (Li.Py 1M 4; — £6, M,;) d; 
= BiB ry d;/2h = B;,b;,/2h. 


Since the only elements of 8; which are not zero are those given in (28), 
the non-vanishing stresses are 


090 = (Byob. 5 Bogb3)/2h, 


093 = (Ba3D2 + Bg3b3)/2h. 
These stresses clearly vanish if 6; is given by (61), and by (31) this is the 
only value of 6; for which they vanish. Thus the simple tilt boundary 
normal to the slip plane is the only type of boundary which can be formed 
from dislocations all of the same kind, and for which the remote stress 
field is zero. This agrees with Frank’s results. 


Fig. 3 


and 


ee B 


a b 


It remains to reconcile the treatment here with that of the previous 
section. Consider a wall built up of a finite number of dislocations. 
First suppose that the wall occurs in an infinite crystal (fig. 3(qa)); the 
treatment of § 5 will now be valid; but the present section fails to consider 
the regions of high stress concentration M and N at the ends of the wall, 
and so cannot give any information of the total energy and the equi- 
librium orientation of the wall in this case. Secondly, consider a finite 
erystal with the wall extending right across it (fig. 3(b)); now the surface 
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of the crystal must be stress free and, to secure this, additional stresses 
(which possibly alter the equilibrium positions of the dislocations) must 
be added to the stresses of §5. On the other hand suppose the stresses 
in the crystal are those obtained in this section, that is they are given by 
eqns. (51) and (52) with z,-1 replaced by the expression (63); with the 
Burgers vector given by eqn. (61) the stress on the surface of the crystal 
is then negligibly small except within a distance of order h from the 
points A and B where the surface is met by the wall, and the stresses 
here are easily shown to give zero resultant force and moment; hence 
if the stresses on the surface are relaxed the change in the stress field will, 
by Saint-Venant’s principle be confined to regions at the ends of the 
wall with dimensions of order h. Also the stress o,, on each dislocation 
due to the others may be shown to be zero, and so the dislocations are 
in equilibrium. Thus the stress field has essentially been found. 

We see then that there are two distinct physical situations: if the wall 
terminates inside the crystal, its orientation must be determined by the 
considerations of §5; while if the wall extends right across the grain, the 
treatment of the present section will apply, and the wall will be normal to 
the slip plane. 

However if the angle 0 (fig. 2) defining the equilibrium between a pair 
of edge dislocations is a right angle, then incomplete walls will also be 
formed normal to the slip plane and the distinction between two cases 
largely disappears. Owing to the high symmetry of the lattices com- 
monly occurring in metals, this will be so in a number of important 
instances. For if either the x, or the x, axis, i.e. either the slip direction 
or the normal to the slip plane, is a two fold axis, then the equilibrium 
positions of the dislocation at P must be arranged symmetrically on 
either side of the plane 2,=0; also the total number of equilibrium 
positions was seen in §5 to be odd; hence an equilibrium position must 
occur on zv,=0. ‘The only important case not covered by these symmetry 
requirements seems to be that of the body-centred cubic lattice with slip 
direction [111] and slip plane (112) or (123) (but not 110)). Here dis- 
location walls may be expected in different orientations according as they 
are complete or incomplete. 

On the other hand dislocation walls terminating inside the crystal may 
not occur with any great frequency as, owing to the high stress concen- 
trations at their ends, they represent a state of high energy. They are un- 
likely to survive any prolonged heat treatment in which dislocation climb 
can occur; while even if only glide occurs, the nature of the stresses at the 
ends of the wall is such that dislocations will tend to be captured, and in 
this way the wall can extend across the crystal. 


§ 7. EVALUATION OF THE CONSTANTS 


In this section we consider how the various constants we have intro- 
duced, and particularly the L;,, may be evaluated. It will be convenient 
to use instead of the elastic moduli c,,,, the elastic coefficients Sy 
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This also has the advantages for numerical work that compilations of 

experimental values, for example those of Schmid and Boas (1950), and 

of Hearmon (1946), usually give the values of s,,, rather than of c Pies: 
The relation between stress and strain is then 


Car = Sypwon- ° . ° . ° ° ° e (64) 


In the present case the strain e; = 0u3/0x, is identically zero, and we may 
use this to eliminate the stress component o;=0,; from eqns. (64); we 
obtain 
Qi 
eur — SanOy> . . . . . . . . (65) 
where 
. £ 
Suw =Sun — 81738 3/ Sag: . . . . . . (66). 


Equation (66) gives S,,3;=S3,.=0, and omitting these elements from the 
matrix [S,,,], we may show that the remaining 5 x 5 matrix is reciprocal 
to (30). 

The constants §;; and y;; are now easily expressed in terms of the S,,,; 
from Jacobi’s theorem we obtain the non-zero elements in eqns. (28) and 
(34) 


20 Saa/(So2Sa4— So”), | 
Bs3 = So9/(Sa9Sa4— Soy”), (67): 
Bos = B32 = — Syo4/(So2Sq4— So4”), | 
and 
¥11 = S55/(S11S55— S35”): | 
¥33 = S1a/(Su855— S15”), : (68) 
¥13 = Ya1= — S45/(SuS55— S,5”). | 


We now consider solutions of the elastic equations which involve only 
one of the roots p,, that is we take the displacement to be of the simple 
form (6) rather than the more general form (9). The stress functions ¢; 
may then be written, on using (16), as 


d= — Palioaf (L1 + Pala): ] 
$2 = Ly, f (1+ Poo); ! od pce tek oe GEN 
b3 = Li f(t + Pa%2)- | 

The stresses deduced from these functions must satisfy the compati- 


bility conditions. Expressed in terms of the strains, the compatibility 
equations reduce, in the present case, to 


As ee 670) 
Ao GREP OS shGE RC 
and 


Oe, 0s _ ae eee ne 
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Inserting (65) in (70) and (71) and using (10) and (69) we obtain 

Log lDotS 41 — 2942516 + Pa(25 12 + See) — 29,526 + So] 

— Lay Pa? S15 — PalS1a+ S56) + Pa(So5 + S46) —Soal=9, . (72) 
and 

Log { Da? S15 — Pa(S14 + S56) + Pa(So5 + Sag) — Soa] 

— Dey Po? S55 — 2PaSa5+Sag]=9. - (73) 
Eliminating L,, and L,, from (72) and (73) we find that p, must be a root 
of the sextic equation 

[p*Syy — 293515 + p?(2S yo + Seg) — 2PS26 + So9][D7S55— 29S 45 + Sas] 

— [p>Sy5— p?(S14+ S56) +: P(S25+S46)—Soa]P?=0. . (74) 
Now eqns. (8) and (74) both just determine the values of p for which a 
solution with the displacements of the form (6) exists. Thus these two 
equations must be identical and we have simply expressed the coefficients 
in terms of the S,,, instead of the c,,y. 

Now we recall that the vectors L,, are undefined to the extent of an 
arbitrary constant factor (for each «); if L,,40 we choose this factor so 
that L,,=1. Then L,,=(—>p,, 1,1,), where J, is a constant which has 
still to be determined. Substituting this value of L,, in (73) we find 


dy = [Pa S15 — Pa®(S14 + S56) + Pa(So6 + Sag) —So4]/[Pa?555— 202845 +Sa4], (75) 


and 1, is determined provided p, does not satisfy the equation 


PS 20S ce PaO ee ee eT 
Then if none of the roots of (74) also satisfy (76), we may write 
= Pr = Pages 
Le = il ] } > . . < ; . (77) 
L be be 


and 
lz3—l, I3pe—Inp3 P3— Ps 
[M,J=L"| 4-1, lp3—lyp, pi—Ps |, (78) 
l,—l, 1py-lp, P2- Pi 
where L denotes the value of the determinant |Z;,|. To determine B,,, 
we must first evaluate P;; and @,, defined in eqns. (42) and (43) and then 
solve (41) and (44). However, the expressions become increasingly 
complicated, and it is probably best to procede numerically at this point. 
If eqns. (74) and (76) have a root in common we choose the labelling « 
s that this root is p;. Comparing (74) and (76) we see that we must 
have 


D3°S15— P3"(Sy4+ S56) + Pa(So5+Sag)—Soe=0, . . . (79) 
and (72) and (73) will be satisfied if L,. is of the form (0, 0, 1). It is of 
course possible for eqn. (79) to be true without the coefficients vanishing 
identically, but this does not appear to be of much physical significance. 
A more important case is when (79) holds because the coefficients vanish 
identically; a sufficient condition for this is that the x, axis should be a 
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two-fold symmetry axis. Then (74) factorizes into a quartic and quad- 
ratic equation, and also (72) and (73) give 1,=1,=0. Thus 


ee fewer | 


ne , aie ae 


and 
(3-3)? Po(P2—P1) 4 | 
(Mal) —(@e-P)" —p(p.—p,) 4 0}, . 2 . (81) 
0 0 1| 


‘On account of the simple form of these matrices, the remaining constants 
are readily evaluated explicitly. From (43) 


Pit+Ps Pip, OY 


2S id 2 SS es ORR Pe ee uit -7)' 
0 0 Pal 
and from (42) 
0 —1 0 
[Qis]=| Upp, Vpit+lfp, 9 |. « «© « . (88) 
0 0 1/p5 


Also in the present case, the non-zero elements of f;; and y,; reduce, by 
(67) and (68) to 
Bo2= Soot, 33= Saas 
and Ppt eee os) 
yu=Sy 4; y33= S555 
since S,,=S,,=0. Then substituting (82), (83) and (84) in (41) and (44), 
we find 


Sy1(21 + P2— Pi — Pe) S11(P1P2— PrP2) 0 
[B;;] = U S1(91P2— PP) — S59(1/p,+ 1/p.—1/p, — 1/P,) 0 
0 0 S55(D3 — Ps) 


(85) 

By substituting this value of B,, in eqns. (50), d; may be expressed in 
terms of the Burgers vector, and hence the stresses of the dislocation may 
be obtained from (51) and (52) by a straight forward substitution. 

In particular we note that if b;=(b, 0, 0), then d,=0 only if p)p.= P,P, ; 
but this would imply a relation between the elastic coefficients which 
will not in general hold. The statement following eqn. (62) is thus 
verified. 

§ 8. Cracks: THE GENERAL SOLUTION 

Suppose we have an infinite crystal in which there is a crack along that 
portion of the plane «,=0 which is given by —1<a,<1. An arbitrary 
stress is applied to the crystal such that in the absence of any stress 
relaxation due to the crack there would be a traction of 7;(z,) over the 


P.M, gee 
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surface of the crack. We must therefore add a stress field which produces 
a traction of —7,(a,) over the surface of the crack and which vanishes at 
infinity. 

We consider the semi-infinite regions z,>0, and x,<0 separately. In 
a> 0 a displacement which has the form (9) and which remains finite at 
infinity is 


Uy, = > Alga [ apr a' (p) Xp (tpz,) + ¥ Are | aeF.*(p) exp(—1pz,). (86) 
The ean stress functions are | 

b= DLps | del s'(0) exp (ip2=) + SLs |, doF "exp (—iv%.). (87) 
In the ati 7m <0, we may take the ee functions to be 

$i. = d Li | dor a (p) exp (—tp2,) + Dhow | 40F) exp (7pZ,)- (88) 


We have now to relate the solutions in Re two regions. In the part of 
the plane z,=0 outside the crack the material must be joined together, 
and the boundary condition on the stress is then that the components 
o;, Should be continuous. In the crack the components take on pre- 
scribed values and so are a fortiore also continuous here. Hence o,.= 
0¢,/0x,, is continuous over the whole plane 7,=0. Now an arbitrary 
constant can be added to ¢; without affecting the stresses, and if this 
constant is chosen suitably then 4; will also be continuous across 2,=0. 
Putting v,=0 in (87) and (88) and comparing the resulting equations we 
obtained the boundary conditions 


2 Lin 7 aie = LLP) 4 


denoting these oe by ane) , we have 
ite *(p) M h;(p) ] 
. (89) 
i “(p)= MihiAp). 


On substituting from (89) in (86) we obtain the displacement in the 
region 27> 0 


U,= > A, M ly bi(p) exp (tpz,) dp + 2a hip) )exp (—%pz,) dp. (90) 


In the region 2, <0, the displacement will be 
uUyn= 2 AicMas le W(p) )exp ( — tpi, )dp+ > ApMas |v b(p) exp ( ipey )dp. (91) 


Subtracting (90) and (91) we find that the difference in displacement on 
either side of the plane 2,=0 is 


Aw;,= — 20B;,; ie {¥i(p) exp (ip2,)—y,(p) exp (—ip2,)}dp, . (92) 


where B,; is defined in eqn. (40). Outside the crack the displacement 
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must be continuous and so Au,=0. Since |B,,|40, this gives 
| 0 (ilp) exp (ipa) —Y,(p) exp (—ipx,)}dp=0, |a[>1. . (93) 


Also substituting (89) in (87) we obtain the stress functions in the region 
Xy> 9 


b= LLMs | “Ule) exp (0pz, fap > teal if “Uie) exp ( — pz, ) dp. (9 4) 


Now on the surface of the crack the stresses o;,. = 0¢;,/0x,, must be equal to 
—T;(x,). Hence from (94) we obtain 


iy {,(p) exp (ip2,) —x;(p) exp (—ipx,)}pdp=—7,(a,), —1<a,<1. (95) 


If we write ¥;=7,' + i#;” where y%,’ and %," are real, eqns. (93) and (95) 
become 


| {yj sin px, + ¥;" cos pay} dp =0 [ra|> 1, 

0 

and poo 

| {w, sin px, +o,” cos px,}p dp = 47,;(x) —l<a,<1. 
0 


These are equivalent to the two pairs of dual integral equations 


fe us; (p) Sin px,dp =0 ig Soll 
pre we\26); 
fs : ;'(p) sin px, .pdp=}{7,("1)—7(—%)}, 0<a%,<1 
and z= : 
ie is," (p) Cos px, dp = 0 LO; 
(97) 


ir b;"(p) cos px, .pdp=4{7,(%1) +7(—2)} O<ay< 1. 


y} 
Equations (96) and (97) are special cases of a general set of dual integral 
equations studied by Titchmarsh (1937) and Busbridge (1938). Their 
general solution leads to 


asi ( aoe uJ (up Jade fn 7(ue(1—2)-MEdE, . . (98) 
and 
wi," (p) = Te =| pS (up Jdu {7 7(mE)(1—&)-H2dé. =... (99) 

On substituting (98) and (99) in eqn. (90) we obtain the displacements 
ty= od AneMas+ Apallles) | du | rl ode. — 2) Ag 
on “ = 
= PE AeaMas | dale? —p2y7" |" ryled)(u +2061) de 
TT 0 
ES Apa Mas [ae e—n2)7" | aus) 2.-B) PAE. (100) 
Ty 0 jai 


2X2 
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This expression will be valid both in the region x,>0, and in x, <0 pro- 
vided we choose the sign of (z,?—?)"? so that 


arg (2,2—p?)"2>argz, as |z,|—> 2. 


The stress functions 4, will be given by an expression of the form (100) 
but in which A,, is replaced by L,,,. Once the ¢;,, are determined the 
stresses are readily obtained by differentiation according to eqns. (10); 
it is to be remembered that the stresses obtained in this way are to be 
added to the stresses which would occur in the infinite material if the 
crack were absent. 

A particular case of some importance is that in which the applied stress 
field is uniform. With 7, constant, the equation corresponding to (100) 
for the stress functions reduces to 


Pie = 4 » {LM 4; ere r= Tyne ati 2a] a Digg M a [(2,2 an | Wis ee 2] \T;. 


a 


From this we obtain the stresses 
Cy 3 >: {La M aj Pa [Za(%" ca ee oe 1] g Lig M us Ba [Z4(Z" a LY tas ot ]}7;, 
a 


on2=% > OIE aN Ce eee Lig a Bae LD at, ads Fe 
a 


aja 


§ 9. THE STRESSES NEAR THE TIP OF A CRACK 


The stresses near the tip of the crack are of special importance since 
there will be a large stress concentration here, and we shall consider now 
the nature of these stresses. : 

First we reduce the repeated integral in (100) to a single integral. This 
may be accomplished by replacing the variable € by y=€u and changing 
the order of integration; we obtain 


b,= — => {LM ia 7;(€) dé cost 2 
+ Dig M ag ee 7;(€) dé cos al Oneeks SHLD 


on dropping a constant term which can contribute nothing to the stresses. 
Then from (101) we obtain the stresses 


ep) {LiPaM a(e2— ye | (6) deep", 

v4 =a! 
+ LabeM (G2 — 1) |) “(6 de BG, 2), _ (102) 
and i. 


Y [inaMajeat 1% | (6) age), 6) 


+ Dagl E21 ayeyaec—eym@,—esb. (108) 
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Now we are interested in points near the tip of the crack where A ale 


hence we write z,=1+¢, and take \¢,,| small compared with unity. If 
r and @ are defined as in fig. 4 we have 


C, =r(cos 0+ p, sin 6). oi deere laa (Lod) 
Then 


Ale 


rE) — EM, 8) dem = | (e\(1-+ Eg), 


/—1 / =—1 


=T',, say. Pica Ke LOD) 


7’; may be regarded as a suitably weighted average value of 7; over the 
crack ; it has Peers defined so that if 7; is constant 7; =7,. Substituting 


(105) in (102) and (103) and remembering that ¢, is Snel we obtain 


25 5 as + » ig PM (26) =F: Ln, DM ,,(26,) YT ;,, 


ef ee = : 106 
Ops = FY (Lip gy 2Lg) 2 + L pelt (22) 1}7,,. | es 


Fig. 4 


In particular we note that these stresses decrease with the distance from 
the tip of the crack as r—!?, just as in the isotropic case. 

In the plane of the crack we have 0=0 and ¢,=r, and so from (106) the 
stresses acting across this plane are 


Crean) J hene . . ° . 5 . ° (107) 


it is interesting to note that these stresses depend neither on the symmetry 
nor the elastic constants of the material. 
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If the length of the crack is 2c then we have only to replace 7 in eqns. 
(104) and (107) by r/c. 


§ 10. Tum ENERGY OF A CRACK 
If the two sides of the crack undergo a relative displacement Aw I when 
the surface tractions 7, are relaxed, the elastic energy of the material will 
change by an amount 
rl 
Bey 7 \Au, dig ane LOS 
d=] 
Now Au, is given by eqn. (92) (—1<a,<1), and substituting from (98) 
and (99) in this we obtain 


2 (i dp re dé 
Au,= = By | came) ue 6) Gas . (109) 
From (108) and (109) the energy of the crack may be written as 
1 
U=IBy | WeldeP (Tw), © - - » (110) 
where 
1/2 
Pi(p)=a | 7(usiné)(1+sin6)d0. . . . (111) 
—7/2 


T,, as defined in eqn. (105) is related to the present 7'j(u) by 


In the particular case in which the applied stresses are uniform, 7',() =7; 
and the energy (110) becomes 


or if the crack has length 2c then this energy is 
U= 37B j j:7j7 pC. . ° . : . ° . (112) 


Similarly, in the general case a factor c? must be inserted in eqn. (110). 
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The Capture of Slow Neutrons by Protons} 


By A. R. Baker} and D. H. Winxrysong 


Cavendish Laboratory, Cambridge 
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ABSTRACT 


The result of a new experimental determination of the capture cross 
section of protons for slow neutrons is briefly reported. It is oy =0-327--0-004 
barns and is obtained by a comparison with the absorption of boron for 
slow neutrons. 

A survey is presented of various direct and indirect methods for determining 
the neutron—proton capture cross section and the final value resulting from 
their inter-comparison is og=0-3315+0-0017 barns. 


THE cross section o,, for the capture of slow neutrons by protons is the 
simplest process which may be expected to be markedly influenced by 
mesonic interaction effects (Austern and Sachs 1951, Austern 1953). Since 
these effects are expected to give a contribution of a few per cent to the 
cross section we must know o, with high accuracy. Other uncertainties 
in the evaluation of the interaction effect are due to lack of precise 
knowledge of other relevant two nucleon properties such as r,(—, — €), the 
D-state percentage in the deuteron ground state and the effective n—p 
potentials. However these uncertainties are presently being reduced by 
both experimental and theoretical study and this is a suitable moment 
for an accurate evaluation of o,. 

Several measurements of co, have appeared during the last decade. 
None however has used the traditional method of comparing the slow 
neutron density distributions due to a steady neutron source in water 
and a boron-rich solution. This method has the disadvantage that the 
bulky detectors usually used introduce a perturbation into the neutron 
distribution which cannot be evaluated with adequate reliability. We 
have however returned to this method and have overcome the difficulty 
by using as slow neutron detectors unmounted discs of boron loaded 
nuclear photographic emulsion of diameter 6mm and thickness 100 yp. 
The perturbation in this case is small and can be accurately calculated 
owing to the simple form of the detectors. The neutron distribution was 
determined by counting under the microscope the individual events due 
to the reaction !B(n, «)7Li in 54 detecting discs in each of two cylindrical 
tanks of height 4 ft and diameter 5 ft, one filled with distilled water and 
the other with an aqueous solution of boric acid. The neutron source was 
a Po-Be mixture. A total of about 200 000 tracks was counted. As a 
nee ene ee 
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result of this investigation we determined the ratio between o, and o, 
the cross section for slow neutron capture in our boron. The ratio between 
oy, of our boron and of the standard Harwell boron was measured by a 
pile oscillator method. o, for the standard Harwell boron is taken as 
769-4 + 3-8 barns (at a neutron velocity of 2200 m/sec) from the work of 
Egelstaff (1957) which considers the work of many laboratories. 

From this experiment we find the value o,,=0-3271 + 0:0042 barns. A 
full account of this work will be published by the first named author of 
this note. 

We now attempt a comparison of the various data from which o,, may 
be derived in order to arrive at the best current value for this constant. 
There are three chief approaches. In the first the mean life 7 of thermal 
neutrons in water is measured directly. In the second a measurement is 
made of the ratio o,/o, which is then converted into a value for o, by 
using some independently determined value of og. In the third we 
combine measurements of L (the diffusion length for slow neutrons in 
water) and D (the diffusion constant for slow neutrons in water) at the 
same temperature in order to arrive at a value for the mean life 7= L?/D. 


Table 1. Direct Determinations of the Slow Neutron Mean Lifetime +r 
in Water Reduced to 22°c 


Author 


Antonov et al. (1955) 206-5 + 4-2 
Bracci and Coceva (1956) 202 +6 

von Dardel and Sjostrand (1954) 204-4 + 2-0 
Meads et al. (1956) 203-4 + 2-6 


Mean value r= 204-3 + 1-4 psec 
(y2=0-54 ; P=0-9) 


In table 1 we present the results of recent measurements of 7. We 
have rejected the measurements by von Dardel and Waltner (1953) and . 
Scott ef al. (1954) the former owing to the unknown and probably large 
effect of the perturbation by the detectors of the slow neutron density 
(see von Dardel and Sjéstrand 1954) and the second because of the 
uncertain allowance for the effect of higher harmonics. Where necessary 
slight corrections have been applied to the original results in order to 
reduce them to an effective water temperature of 22°c. This correction 
is particularly important for one of the measurements presented by 
Antonov et al. (1955) which was carried out at 80°c. The lifetime quoted 
for the results of these workers is obtained by averaging that determined 
at 23°c with the suitably corrected figure at 80°c. It is seen that these 
several results are in excellent accord with one another. 
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In table 2 we summarize the recent accurate determinations of o,/loy 
including that reported in this paper. The two U.S. results were originally 
determined relative to the standard U.S. boron. In quoting o,/o, in 
table 2 we have corrected this ratio to correspond to the standard Harwell 
boron using o, U.S./o, Harwell = 0-986 + 0-003 suggested by the work of 
Green et al. (1954)—see also Egelstaff (1957). The mean value of op/oy 
given in table 2 therefore refers to standard Harwell boron. Again the 
accord between these various measurements is excellent. 


Table 2. Measurements of o/c, reduced to standard Harwell boron 
using the Ratio o, U.S./o, Harwell =0-986 + 0-003 


Author Cnlon 
Baker and Wilkinson (this paper) 2352 + 28 
Hamermesh ef al. (1953) 2325 + 28 
Harris et al. (1953) 2306 + 51 


Mean value og/oy = 2338 + 18 
(7=0'86 "7 P=0-7) 


Table 3 summarizes recent measurements of the diffusion length L. 
In presenting these results we have corrected the quoted values for 
buckling where necessary and have reduced them to a value appropriate 
to 22°c on the assumption that L is proportional to the square root of the 
absolute temperature. We have admitted without further correction the 


Table 3. Measurements of the Diffusion Length for Slow Neutrons in 
Water Reduced to 22°c 


Author LI (em) 
Barkhov et al. (1957) 2-684 + 0-020 
De Juren and Rosenwasser (1953) 2-735 + 0-015 
Wright and Frost (1956) 2-700 + 0-010 

Mean value L=2-707 + 0-008 cm 


+0 
(2=5:3; P=0-1) 
iar 


preliminary value due to Wright and Frost (1956) which may require 
subsequent revision. These values of L are seen to be in barely adequate 
accord with each other and for this and other reasons connected with the 
buckling correction and possible systematic errors we have preferred to 
increase the statistical error by 50%. 

In table 4 we present recent results on the diffusion constant D. These 
we have reduced to an effective temperature of 22°c under the assumption 
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that a 1° rise in temperature increases D by 130 cm?/sec. This correction 
is a very small one and is in all cases well below the stated error of the 
measurements. These determinations are seen to be in good mutual 
accord. 


Table 4. Measurements of the Diffusion Constant D for Slow Neutrons 
in Water Reduced to 22°C 


Author D (em?/sec) 
Antonov et al. (1955) 34 870 + 1000 
Bracci and Coceva (1956) 34 850 + 1100 
Campbell and Stelson (1956) 34 670 + 1000 
von Dardel and Sjostrand (1954) 36 340 + 750 
Dio and Schopper (1958) 35 440 + 600 


Mean value D=35 420 + 370 cm?/sec 
(=e. P=O7) 


When these values for ZL and D are combined we find the result 
T= 206-9 + 2°8 psec. 

In table 5 we compare the values of o, obtained by the three inde- 
pendent methods. In quoting the value derived from the ratio op/o, 
we have used the above-mentioned value of o, for the Harwell boron. 


Table 5. Summary of Determinations of o,. In quoting o, from 
op/o, the value op = 769-4 + 3-8 barns for standard Harwell boron 
has been used 


Method oy (barns) 
Lifetime 0-3335 + 0-0023 
On Oy 0-3291 + 0-0030 
Land D 0-3293 + 0-0045 


Mean value og =0-3315 + 0-0017 barns 
(2=163-P=0-5) 


We see that these results are in excellent mutual agreement and since they 
are independent of each other we quote the final mean value: 


oy, = 90-3315 + 0-0017 barns. 
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A Simple Formula for use with Carbon Thermometers 
at Low Temperatures 


By O. V. LoUNASMAAT 
Clarendon Laboratory, Oxford 


[Received March 10, 1958] 


SEVERAL empirical or semi-empirical formulae have been proposed for 
use with carbon thermometers at low temperatures (Brown et al. 1951, 
Clement and Quinnell 1952, Lacaze and Peretti 1953, Hoare et al. 1955). 
Unfortunately, the mathematical forms of these are rather too complicated 
to be convenient when the amount of experimental data to be handled 
is large. A simpler formula was therefore sought, with the particular 
requirement that it be valid in the temperature range 4 to 12°K where 
frequent calibration is usually inconvenient. 

The thermometer used was a 47 ohm, 0-25 watt radio resistor, taken 
from a ‘ Lab Minipack ’ card supplied by V.E.S. Wholesale Ltd., London. 
After removing the paint from the surface, an insulated copper wire was 
wound round it and bonded in place by Araldite 985E. The ends of the 
wire were soldered to a calorimeter whose temperature it was ultimately 
desired to measure ; for present purposes it served simply as a gas 
thermometer bulb. This assembly formed part of an apparatus used 
for determining the entropy diagram of helium (Hill and Lounasmaa 
1957). The details of the apparatus do not concern us here except that, 
as in most low temperature apparatus, thermometer calibrations against 
helium and hydrogen vapour pressures could conveniently be included 
in every experiment, whereas a gas thermometer calibration at inter- 
mediate temperatures involved a slight modification of the apparatus 
and so necessitated a separate experiment. 

A first calibration in the range 4 to 12°K showed that when 1/R 
(&=thermometer resistance) is plotted against temperature an upward 
convex curve is obtained. On the other hand, 1/R’, where R’= R— diy 
is nearly a linear function of temperature if the constant Ry is suitably 
chosen. For a thermometer which had a resistance of 176-3 ohms at the 
hydrogen triple point (13-957°K) R, was found to be 94:0 ohms. When 
the calibration was repeated some months later, the thermometer having 
been cycled repeatedly between room temperature and helium tempera- 
tures in the meantime, little change was found. The resistance at the 
hydrogen triple point had increased by 1-6 ohms, but if R, was increased 
by the same amount the two curves of 1/R’ against temperature coincided 
within the experimental accuracy. This is illustrated in the diagram, 
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where the temperature derivative of 1/R’ is plotted as a function of 
temperature for the two runs. 

As the quantity d(1/R’)/dT is sensibly constant, both with temperature 
and between experiments, the required relation may be written 


T+0=A/R' 


where @ and A are constants for a particular thermometer ; in the present 
case, 0+ —1°K and A~1000 ohm deg. It is to be emphasized that 6 and 
A do not vary from run to run, provided that in each run R, is adjusted 
so as to bring R’ to its standard value at one convenient temperature, 
preferably the hydrogen triple point. If the resistance measurements 
are to be made with a Wheatstone bridge, R’ may be measured directly 
by inserting a resistance equal to R, in the balancing arm. The reduction 
of data then becomes particularly rapid, using either a desk calculator 
or tables of reciprocals and a large scale plot of 1/R’ against 7’. 


1.00 


Ra 


£ 11% 


olexp.; ell exp. 


0.95 


This calibration formula was tested and confirmed during a large 
number of experiments. With slight modifications its use can be extended 
to other temperatures. 
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The Nucleation of Dislocation Loops during Cleavage 
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Ir has been shown recently by Gilman (1956) that dislocation loops are 
nucleated by a propagating crack if the speed of the crack front falls below 
a certain critical value. Gilman deduced the presence of the loops from 
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etch patterns. The same phenomenon has been studied also by Forty 
(1957) who used interferometry and microscopic examination, and in- 
ferred the existence of loops from the geometry of the cleavage face. 

We have developed recently a method for the decoration of dislocations 
in KCl which leaves the dislocation pattern after deformation practically 
undisturbed. This makes it possible to observe the dislocation loops in 
an even more direct way. Photograph 1 represents an example. The 
erystal was intentionally cleaved in such a way that the crack front stopped 
a number of times. The instantaneous position of the crack front can 
clearly be recognized on the photograph 1. 

It is evident that a number of small dislocation loops have been formed. 
They are lying in two families of (110) planes. A large number of them are 
single loops. A few systems of concentric loops are however also visible, 
suggesting that surface sources have been active during a short time. 

The decoration procedure is as follows. One uses KCl doped with 
0-75% of AgNO, (by weight). The dislocations are already slightly 
decorated in the crystal as grown. A much better decoration, even of 
freshly introduced dislocations, can however be obtained if such crystals 
are x-irradiated for +6 h in contact with the window of a copper target 
x-ray tube operated at 40 kv, 20 ma, and subsequently annealed in air 
during 2 h at 600°c. At high magnification and in crystals irradiated and 
annealed for a longer period than those of which the patterns are shown, 
it is observed that the decoration consists mainly of small cubic cavities, 
similar but smaller in size than those found in hydrogene treated specimens 
(Amelinckx et al. 1958). The cavities are filled with gas under a pressure. 
of a few atmospheres, and the walls are covered by silver. This can be 
deduced by dissolving the crystal in water and observing it under the 
microscope. Every time the dissolution front reaches a cavity a small gas 
bubble, a few times larger than the cavity giving rise to it, is produced. 
The nature of the gas was studied by means of mass-spectroscopy and by 
chemical means; details are given elsewhere (loc. cit.). It mainly 
consists of oxygen and nitrogen, resulting from the decomposition of the 
NO,- and NO,~ radials under the influence of the ionizing radiation. The 
presence of oxydes of nitrogen could not be established. 

The method has serious advantages, in that it leaves the pattern of freshly 
introduced dislocations largely undisturbed, probably as a consequence 
of their pinning under the influence of x-irradiation. Even the anneal 
at 600°C hardly disturbs the pattern. In particular the small half loops 
do not pop out as is invariably the case if they are not pinned. It is clear 
that this method is a promising tool for investigating the dislocation 
configuration after slip. Photograph 2, for example, shows the first 
stage of polygonization in KCl, i.e. the rearrangement of edge dislocations 
in their glideplane, so as to form vertical rows. Edge dislocations in their 
glideplane have also been observed in potassium bromide crystals decorated 
by means of gold (Barber et al. 1957). 

A more detailed account will be published elsewhere. 
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THERE has been much recent interest in the problem of the low fatigue 
properties exhibited by age hardened aluminium alloys. Current 
theories suggest that the poor properties are associated with the relative 
instability of the dispersed precipitate under conditions of cyclic stressing. 
Migration of solute atoms and over-ageing are thought to occur in heavily 
deformed localized regions. This produces zones which are relatively 
soft and in which subsequent plastic deformation tends to be concen- 
trated leading, ultimately, to accelerated fatigue cracking. 

There is some support for the suggestion that moving dislocations 
generate vacancies which increase solute diffusion rates in localized 
regions (Broom et al. 1956). Broom et al. (1957) have indicated that this 
proposal requires any structural changes to be associated with slip planes. 
Furthermore, since vacancies are generated during unidirectional extension 
as well as during fatigue, then some changes might be apparent in fractured 
tensile specimens. It was with the object of detecting such changes that 
these workers studied high purity alloys of aluminium—5-59, zinc—2-65°% 
magnesium and aluminium-—6%, zinc—3°%, magnesium—1°% copper. There 
was some evidence of slip band effects in polished and etched sections of 
the simpler ternary alloy but, because of a tendency to intercrystalline 
cracking in this material, work was concentrated on the quaternary alloy 
which had the basic composition of the commercial material DTD 683. 
Intense slip bands were observed in polished and etched sections from the 
interior of both tensile and fatigue specimens. However, modifications. 
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to the distribution of the precipitate, as shown by Hanstock (1954) in 
fatigue-tested DTD 683 specimens, were not observed in either case. 

We have also been studying metallographic changes associated with 
fatigue and tensile stressing of a series of high purity aluminium—zinc— 
magnesium alloys. Standard 0-25 in. diameter tensile specimens were 
machined to a diamond-turned surface finish from forged bar which had 
previously been heat treated. The heat treatment cycle, which consisted 
of 14 hours at 460°c, water quenching, and age hardening for 16 hours at 
125°c, was chosen so that all alloys had a common U.T.S. of approximately 
25 tons in?. The specimens were tensile tested to failure, carefully centre 
sectioned in a longitudinal direction, and mounted in a cold setting resin. 
The centre section was mechanically polished and etched in Modified 
Kellers reagent (87-59% H,O, 10% HNOs, 1:5% HCl, 1-0% HF) which 
clearly showed the distribution of the precipitate and any deformation 
markings. The structural changes observed in a tensile specimen of an 
alloy of aluminium—5-1°% zinc-3-0°% magnesium are of interest in view 
of the earlier remarks. 

The room temperature tensile test of the alloy gave the following 
results: 0-1°% P.S., 20-6 tons in?; U.T.S., 26-3 tons in?; and elongation 
on a gauge length of 1} in. of 17-7%. Significant metallographic changes 
were confined to the necked section. In some grains widely spaced 
bands were observed which, because of their crystallographic nature, were 
presumed to indicate regions of intense slip. An example of this is shown 
in fig. 1, Pl. 27 where it is the distribution of the precipitate which is of 
particular interest. Intense precipitation has occurred in the region of 
the slip bands with consequent depletion of the surrounding matrix; a 
process which must have occurred during the comparatively short test 
period of approximately three minutes. Closer examination showed that 
the lines of precipitate and the coarse slip bands did not coincide. This is 
clearly evident in fig. 2, Pl. 27 where an actual plane on which slip dis- 
placement has occurred is revealed by a step in the grain boundary. It 
thus seems that coarse slip has occurred through the softer, depleted 
matrix which is the same result as that obtained for certain age hardened 
aluminium alloys tested in fatigue (Forsyth and Stubbington 1957). 

A further example of precipitation caused by deformation in a tensile 
test is shown in fig. 3, Pl. 27. In this case the slip lines have been 
generated at the tip of an advancing internal crack and the distribution 
of the precipitate is similar to the bands of precipitate which Hanstock 
(1954) has observed in fatigue tested DTD 683 specimens. 

The metallographic changes described above were less evident in alloys 
with a higher zinc : magnesium ratio where failure became progressively 
more intercrystalline or in alloys with a lower ratio of zinc : magnesium 
where deformation was more homogeneous. It thus appears that this 
alloy is particularly susceptible to structural changes when subjected to 
plastic deformation. This conclusion is significant as most of the 
commercial alloys based on the aluminium—zinc-magnesium system have 
closely similar zinc and magnesium contents. 
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Measurements of Thermoelectricity Below 1°k 
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RECENT measurements of the absolute thermoelectric power S of the 
alkali metals down to ~2°kK (MacDonald et al. 1958) have shown that, 
even down to so low a temperature, there is very considerable variety in 
Fig. | 
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the behaviour. In particular, the Thomson heats » derived from our 
results (= 7' dS/d7') show, in the case of rubidium and caesium, a rather 
remarkable oscillation with temperature; in caesium the Thomson heat 
shows a positive maximum at ~ 3°k with a value as high as ~ + 2-5 wV/°c. 

We therefore decided to extend these measurements of thermoelectric 
power to the temperature region below 1°K; as far as we know, no 
measurements of thermoelectric power below 1°K have previously been 
reported. Before attempting to work with the alkali metals themselves, 
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we thought it best to try some preliminary experiments in this tem- 
perature region on the ‘noble’ metals, gold and silver, which are much 
easier to handle than the alkali metals. A few: experiments on gold and 
silver specimens are now complete and the results appear of considerable 
interest in themselves. 

Figures | and 2 show the observed thermoelectric force on two samples 
of gold and one of silver. The electromotive force itself can be measured 
with rather good precision (around 1%). No great accuracy, however 
is claimed for the temperature measurement in these experiments, 
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Ferric Ammonium Alum was used as the refrigerant, and the correction 
from 7* to 7' was made using the recent scale resulting from the experi- 
ments of Cooke eé al.t (1956); only a very rough form factor was used, 
however, for the salt pill which was an irregular cylinder. 

It will be seen that the gold specimen referred to as A shows an absolute 
thermoelectric power of about —4V/°c around 0-5°x. Now the usual 
theory of electron transport predicts for a free electron model: 


CS Te et eer eS 


where C, is the electron specific heat per unit volume, and N the density 
of free electrons. For gold this would give S~ —3-5 x 10-9 volts/°c at 
~ 0:5°K, i.e. less than one thousandth of the observed value. 

Of recent years there has been considerable interest in the ‘ phonon- 
drag ’ effect to which attention was first drawn by Gurevich (1945, 1946). 
This recognizes that the heat flow in the /attice (i.e. a ‘ phonon-current ’) 
can, through scattering of phonons by conduction electrons, ‘ drag ’ these 
electrons with it so causing a thermoelectric current. The magnitude of 
this effect depends primarily on the lattice specific heat, C,,,,, and a crude 
analysis indicates a maximum thermoelectric power of the order Cra 
due to the ‘ phonon-drag’. In semiconductors it appears established 
that the ‘ phonon-drag’ can play quite a major role, and it also now 
appears quite probable that in metals at fairly low temperatures the 
effect may be important. 

However, at 0-5°K, C;,,,/C., in gold is only the order of 10%, assuming 
reasonable theoretical values; consequently it appears impossible to 
look in this direction for an explanation of the relatively enormous 
thermoelectric power which we have observed at these very low tem- 
peratures. 

A more detailed formula (e.g. Wilson 1953) for the ‘ normal’ thermo- 
electric power in a metal is: 

ga ak*T (dlog n(£) it dlog oe (2) 
= ae dlog B gore ple 


(We note that C,/Ne=n?k?T'/2el, cf. eqn. (1).) For ideal free electrons 
the first factor in the brackets equals }, and the second is presumed (Wilson 
loc. cit.) to vary from } at low temperatures to 3 at high temperatures. 
One would clearly have to postulate enormous values (i.e. of the order of 
103) for one or other of the logarithmic derivatives to account for the 
observed thermoelectric power in gold below 1°K. 

These experimental results seem to us, therefore, to lie beyond the 
limits of current theoretical explanation. The interest of measurements 
in this temperature region is evident. 
a 
+ We are grateful to Dr. A. H. Cooke for letting us have a copy of his data. 
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The magnitude of thermoelectric power observed in these metals also 
clearly points to their possible use as differential thermometers, for 
example, in calorimetry below 1°K. 
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REVIEWS OF BOOKS 


Theoretical Physics (3rd edition). By Gwora Joos and Ira M. Freeman. 

(London and Glasgow: Blackie & Son Ltd.) [Pp. 885.] 70s. 

THIS new edition of the well-known textbook by Professor Joos is as useful 
and well written as earlier editions have been. It covers the whole of theoretical 
physics from vector analysis to quantum theory, including classical mechanics, 
relativity, electromagnetism, optics thermodynamics, statistical aspects of 
heat and the atomistic nature of electrical phenomena. This edition has 
been brought up to date by some short additions including a section on matrices 
and by a revision of the chapter on nuclear physics. The treatment through- 
out is a connected development of selected useful topics rather than exhaustive. 
The mathematics never becomes frighteningly formidable, while the application 
of the theory to well-known physical phenomena is continually pointed out. 
In particular ‘* modern ’’ physies, such as the vector model of the atom and 
detailed wave mechanical calculations, receives its due weight. 

Only the nuclear physics chapter is really disappointing, in that it is general 
and discursive, instead of a theoretical discussion of selected aspects in accord- 
ance with the spirit of the rest of the book. Incidentally it contains one odd 
assertion—* probably every elementary particle has spin 1/2”’. 

The book can be recommended strongly as a textbook for honours students 


and as a reference book for the shelf of any work-a-day physicist. 
Velie 


An Introduction to Fourier Analysis and Generalised Functions. By M. J. 
LicgHTHimy. (Cambridge University Press.) [Pp. 79.] 17s. 6d. 

WHEN talking to physicists about Dirac’s delta function, it is customary to say 
‘ Of course, to the pure mathematicians, this is a very improper function, but 
it would take a lot of explanation to get it right, and in practice there is usually 
no difficulty’. To the bright boys of the class, we can now add “ Professor 
Lighthill’s little book shows the simplest and most elegant way of setting up a 
consistent formulation’. The theory of generalised functions, developed by 
Schwartz and Temple, seems to solve most of the epsilontic problems that arise 
in applied mathematics, without recourse to Lebesque integration and other 
sophistications. It is strongly recommended to undergraduates and research 
students who care for rigour (if not for the rigor mortis of absolutely pathological 
functions !). J.M. Z. 


Magneto-Hydrodynamics. Edited by R. K. M. LAnpbsHoFF. (Stanford 
University Press. London: Oxford University Press.) [Pp.x+115.] 32s. 


Tuts book is essentially an edited collection of the papers by different authors 
presented at a symposium on Magnetohydrodynamics held in California at the 
end of 1956. <s 

Five theoretical papers are followed by seven describing experimental work, 
though in the latter is included one on a hydromagnetic waveguide which gives 
the theory of an experiment only being planned at the time of writing. 

The first paper is an interesting attempt to classify the problems of magneto- 
hydrodynamics according to the relative importance of the many parameters 
involved, The next two are of the descriptive type common in this field 
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(in one, the reader is asked to visualize the use of a rolling pin on a lump of 
cosmic dough). There follow analytical treatments of a shock wave problem 
and the pinch effect. The experimental section is largely concerned with shock 
waves. 

The papers are bound to be of interest to workers in the field, but surely this 
rather heterogeneous collection would have appeared more appropriately in 
the research journals than in a book. Pate e. 


Magnetohydrodynamics (Interscience Tracts on Physics and Astronomy, No. 4). 
By T. G. Cownrne. (New York: Interscience Publishers, Inc. [Pp. vill + 
115.] 


Tus is a book to be commended, short but comprehensive. Its length seems 
inversely proportional to its authority, and careful reading, supplemented by 
the references, should give a good appreciation of magnetohydrodynamics. 

Following a concise presentation of the fundamental mathematics, various 
problems are considered, grouped in respective chapters on magnetohydrostatics, 
wave motion, magnetic fields and instability, dynamo theories, and ionized 
gases. Many of the discussions are conducted in words, with quantitative 
injections from order of magnitude arguments; this reflects the impossibility of 
anything like a comprehensive mathematical treatment of the interlocked 
electromagnetic and hydrodynamic equations. On the other hand, some of 
the more analytical investigations do receive attention, results being quoted 
and derivations hinted at. Most of the theories are criticized, and reasons 
given for their rejection often outweigh those for acceptance; but some are 
commended with phrases like “the suggestion deserves attention’, and 
Bullard’s work on the dynamo theory draws the response * Thus the possibility 
of dynamo maintenance of cosmic magnetic fields can be regarded as finally 
established ”’. 

One evident omission is any discussion of magnetic field effects on shock 
waves, and it is perhaps disappointing that so little is said about the important 
case of anisotropic conductivity. Pe 


The Klements of Classical Thermodynamics. By A. B. Preparp. (Cambridge : 
University Press). [Pp. 159, Exercises and Index.] Clothbound 25s.; 
Students Edition, paper covers; 15s. 

Ir has often been pointed out that thermodynamics is a subject of which it is 
easy to acquire knowledge but difficult to acquire understanding. One can 
readily become familiar with the manipulation of partial derivatives without a 
real appreciation of the concepts on which the structure of thermodynamics is 
based. Dr. Pippard’s book, whilst paying adequate attention to technique, is 
particularly to be recommended for providing the reader with an understanding 
of thermodynamics. Many awkward points, which are glossed over in other 
treatises, are discussed clearly and comprehensively, 

The book contains nine chapters and is restricted to the study of homogeneous 
systems. In the introductory chapter the author points out that classical 
thermodynamics approaches more closely to an ideal logical development than 
any other branch of natural science. The next three chapters which are devoted 
to the zeroth, first, and second laws keep this aspect well to the foreground; 
gpa yn. 1s se to the approach of Caratheodory in developing the second 
aw, as well as to the more usual Carnot cycle method. Chapters 

six consider various consequences and HARK of he cee eos 

dynamics; these include the Maxwell relations, various properties of the equation 
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of state, adiabatic changes, specific heats, the Joule Kelvin effect, radiation, and 
surface phenomena. Chapter five also contains a brief discussion of the third law 
of thermodynamics and Chapter six ends with a detailed account of the establish- 
ment of the absolute scale of temperature. Chapter seven deals with the 
thermodynamics inequalities and conditions of equilibrium, and contains as an 
illustrative example a discussion of the equilibrium between a liquid drop and its 
vapour, 

The last two chapters which discuss respectively phase equilibrium and higher 
order transitions, constitute, in the opinion of the reviewer, the most valuable 
and original section of the book. Applications are drawn widely from modern 
research in solid state physics, and in themselves show how much help can be 
obtained from thermodynamics in elucidating these phenomena. The analysis 
of second order transitions, and the correlation of the difference between first and 
higher order phase transitions with the existence or non-existence of metastable 
states, should dispel the confusion with which this topic has long been surrounded. 

The book is not intended for the complete beginner, but anyone with a 
rudimentary knowledge of thermodynamics cannot fail to derive benefit and 
stimulus from perusing its pages. The Cambridge University Press is to be 
warmly congratulated on the innovation of a paper-backed edition which brings 


the price to a very reasonable level. It is to be hoped that this experiment will 
be continued. C.D: 


Soviet Sputniks. (London: Soviet News, 1958.) [Pp. 52.] ls. 3d. 

Tus booklet gives a popular account of the Russian satellites, based on material 
published by Soviet scientists. As popular accounts go it is a good one and 
well illustrated, but it reads, unfortunately, as if the material has been extracted 
piecemeal with a pair of scissors and assembled with a pot of glue. 

The method of launching is described after a brief historical introduction, but 
most of the book is taken up with the various phenomena that can be investigated 
by means of artificial satellites. 

There is little propaganda and only once are we told how it is “far from 
accidental that it was the Soviet Union which first succeeded in building and 
launching satellites’. This is an opinion with which all may concur. 


J.R.S. 


[ 664 ] 


BOOK NOTICES 


Pressure Measurement in Vacuum Systems. By J. H. Leck. (Published on 
behalf of the Institute of Physics by Chapman & Hall Ltd.) [Pp. 144.] 30s. 

Communication, Organization, and Science. By J. Rorustern. (Indian Hills, 
Colorado : Falcon’s Wing Press.) [Pp. Ixxxv+110.] $3.50. 

Contributions to the Theory of Games. Volume III. Edited by M. DresHEr, 
A. W. Tucker, and P. Wore. (Princeton University Press. London : 
Oxford University Press.) [Pp. vi+435.] 40s. 

Edmund Burke and the Natural Law. By Peter J. Sranuis. (Ann Arbor : 
University of Michigan Press.) [Pp. xiii+311.] $5.75. 


|The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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Symmetrical zine bicrystal under tension Gilman (1957). 
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Photograph 1 


Dislocation loops along crack front in KCl crystal. Loops are lying in (110) 
and (110) planes. Note systems of concentric loops ( x 350). 


Photograph 2 


Dislocations in glide planes. Note the tendency to form vertical rows. The 
deformation was caused by indentation somewhere in the bottom left 
corner ( x 500). 
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Precipitation in slip band regions. x 500. 


Fig. 2 


Lack of coincidence between lines of precipitate and slip bands. x 1300. 


Precipitation associated with slip bands generated at the tip of an internal 
erack. x 500. 


